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ABSTRACT 


\ 

A  general  procedure  for  analysis  of  the  response  of  concrete  gravity 
dans,  including  the  dynamic  effects  of  impounded  water  and  flexible  founda¬ 
tion  rock,  to  the  transverse  (horizontal)  and  vertical  components  of  earth¬ 
quake  ground  motion  is  presented.  The  problem  is  reduced  to  one  in  two 
dimensions,  considering  the  transverse  vibration  of  a  monolith  of  the  dam. 

The  system  is  analyzed  under  the  assumption  of  linear  behavior  for  the  con¬ 
crete,  foundation  rock  and  water. 

The  complete  system  is  considered  as  composed  of  three  substructures  — 
the  dam,  represented  as  a  finite  element  system,  the  fluid  domain,  as  a  con¬ 
tinuum  of  infinite  length  in  the  upstream  direction,  and  the  foundation  rock 
region  as  a  viscoelastic  half plane.  The  structural  displacements  of  the  dam 
are  expressed  as  a  linear  combination  of  Ritz  vectors,  chosen  as  normal  modes 
of  an  associated  undamped  dam- foundation  system.  The  effectiveness  of  this 
analytical  formulation  lies  in  its  being  able  to  produce  excellent  results  by 
considering  only  a  few  Ritz  vectors.  The  modal  displacements  due  to  earth¬ 
quake  motion  are  computed  by  synthesizing  their  complex  frequency  responses 
using  Fast  Fourier  Transform  procedures.  The  stress  responses  are  calculated 
from  the  modal  displacements. 

An  example  analysis  is  presented  to  illustrate  results  obtained  from  this 
analytical  procedure.  Computation  times  for  several  analyses  are  presented  to 
illustrate  effectiveness  of  the  procedure. 

The  response  of  idealized  dam  cross-sections  to  harmonic  horizontal  or 
vertical  ground  motion  is  presented  for  a  range  of  inportant  system  para¬ 
meters  characterizing  the  properties  of  the  dam,  foundation  rock  and  impounded 
water.  Based  on  these  results,  the  separate  effects  of  structure-water  inter¬ 
action  and  structure- foundation  interaction,  and  the  combined  effects  of  the 
two  sources  of  interaction,  on  dynamic  response  of  dams  are  investigated, 
leading  to  the  following  conclusions. 

Each  source  of  interaction  generally  has  significant  effect  on  the  com¬ 
plex  frequency  response  functions  for  the  dam.  The  fundamental  resonant  fre¬ 
quency  of  the  dam  decreases  and  its  apparent  damping  increases  because  of 
structure- foundation  interaction.  The  higher  resonant  frequencies  and  asso¬ 
ciated  damping  are  affected  similarly  but  to  a  lesser  degree.  These  effects 


are  qualitatively  similar  whether  the  reservoir  is  empty  or  full,  except  at 
the  resonant  frequencies  of  the  fluid  domain.  Because  of  hydrodynamic  effects, 
the  response  curves  are  complicated  in  the  neighborhood  of  the  natural  fre¬ 
quencies  of  water  in  the  reservoir;  the  resonant  frequencies  of  the  are 
reduced  — ■  the  fundamental  frequency  by  a  significant  amount  but  the  higher 
resonant  frequencies  by  relatively  little;  and  the  fundamental  mode  exhibits 
highly  resonant  behavior.  The  hydrodynamic  effects  in  the  dam  response  are 
qualitatively  similar  whether  the  foundation  rock  is  rigid  or  flexible.  The 
fundamental  resonant  frequency  of  the  dam  is  reduced  by  roughly  the  same  de¬ 
gree,  independent  of  the  foundation  material  properties.  However,  the  appar¬ 
ent  damping  at  the  fundamental  resonant  frequency  is  dominated  by  effects  of 
structure- foundation  interaction  and  varies  little  with  the  depth  of  water. 

The  fundamental  resonant  frequency  of  the  dam  is  reduced  by  each  of  the  two 
sources  of  interaction,  with  the  influence  of  water  usually  being  larger. 
However,  there  are  no  general  trends  regarding  the  comparative  effects  of 
water  and  foundation  on  the  higher  resonant  frequencies  or  on  the  resonant 
responses  of  the  dam.  The  response  of  the  dam,  without  water,  to  vertical 
ground  motion  is  small  relative  to  that  due  to  horizontal  ground  motion, 
but  it  becomes  relatively  significant  when  the  hydrodynamic  effects  are 
included. 

The  displacement  and  stress  responses  of  pine  Flat  Dam  to  the  S69E  com¬ 
ponent  of  the  Taft  ground  motion  only,  and  to  the  S69E  and  vertical  components 
acting  simultaneously,  are  presented.  For  each  of  these  excitations,  the 
response  of  the  dam  is  analyzed  four  times  corresponding  to  the  following 
four  sets  of  assumptions:  (1)  rigid  foundation,  hydrodynamic  effects  excluded; 
(2)  rigid  foundation,  hydrodynamic  effects  included;  (3)  flexible  foundation, 
hydrodynamic  effects  excluded;  and  (4)  flexible  foundation,  hydrodynamic  ef¬ 
fects  included.  These  results  lead  to  the  following  conclusions. 

The  displacements  and  stresses  of  Pine  Flat  Dam  due  to  the  Taft  ground 
motion  are  increased  significantly  because  of  hydrodynamic  effects.  Com¬ 
pared  to  the  response  of  the  dam  including  only  hydrodynamic  effects,  the 
stresses  in  the  upper  parts  of  the  dam  are  significantly  increased  due  to 
structure- foundation  interaction.  Stresses  at  the  heel  of  the  dam  are  in¬ 
creased  to  a  lesser  extent  because  of  the  stress  relaxation  due  to  founda¬ 
tion  flexibility.  The  influence  of  structure-foundation  interaction  and 
structure-water  interaction  on  the  response  of  a  dam  depends  in  part  on 


the  change  in  the  earthquake  response  spectrum  ordinate  associated  with 
changes  in  frequency  and  apparent  damping  due  to  these  interaction  effects. 

As  such,  these  effects  would  depend  on  the  resonant  frequencies  of  the  dam 
and  the  shape  of  the  earthquake  response  spectrum  in  a  neighborhood  of  these 
resonant  frequencies.  Although  considerable  stresses  in  Pine  Flat  Dam  are 
caused  by  the  vertical  component  of  Taft  ground  motion,  they  partially  can¬ 
cel  the  stresses  due  to  the  horizontal  component,  resulting  in  reduced  res¬ 
ponse  when  both  ground  motion  components  are  considered  simultaneously.  The 
contribution  of  the  vertical  component  of  ground  motion  to  the  total  response 
of  a  dam,  including  hydrodynamic  effects,  depends  on  the  relative  phasing  of 
the  responses  to  horizontal  and  vertical  ground  motion,  which  in  turn  depends 
on  the  phasing  of  the  ground  motion  components  and  the  vibration  properties 
of  the  dam. 
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.1 


1.  INTRODUCTION 

1.1  Objectives 

It  is  extremely  important  to  design  dams  which  store  large  quantities 
of  water  to  safely  withstand  earthquakes,  particularly  in  view  of  the  cata¬ 
strophic  consequences  of  dam  failure.  The  damage  to  Koyna  Dam,  in  India, 
which  was  designed  according  to  standard,  widely  accepted  procedures,  shows 
that  concrete  gravity  dams  are  not  as  immune  to  earthquake  damage  as  has 
commonly  been  believed  [1,2]. 

Reliable  analytical  procedures  are  necessary  to  design  earthquake 
resistant  dams,  and  to  evaluate  the  safety  of  existing  dams  during  future 
earthquakes.  These  procedures  should  provide  the  capability  of  evaluating 
the  dynamic  deformations  and  stresses  in  a  dam  subjected  to  a  given  ground 
motion.  Special  attention  should  be  given  to  the  interaction  of  the  dam 
with  the  impounded  water  and  with  the  foundation  rock  or  soil.  These  fac¬ 
tors  complicate  the  otherwise  routine  finite  element  analysis  of  concrete 
gravity  dams . 

The  objectives  of  the  present  study  ares  (a)  to  develop  reliable  and 
effective  techniques  for  analyzing  the  response  of  concrete  gravity  dams 
to  earthquake  ground  motion,  including  effects  of  dam-water  interaction  and 
of  dam- foundation  interaction,  and  (b)  to  examine  the  significance  of  these 
interaction  effects  in  earthquake  response  of  dams. 

1.2  Review  of  Past  Work 

During  the  past  25  years,  the  finite  element  method  has  become  the 
standard  procedure  for  analysis  of  all  types  of  civil  engineering  struc¬ 
tures.  Early  in  its  development,  it  became  apparent  that  this  method  had 
unique  potentialities  in  the  evaluation  of  stress  in  dams,  and  many  of 
its  earliest  civil  engineering  applications  concerned  special  problems 
associated  with  such  structures  [3,4].  The  earliest  dynamic  finite  element 
analyses  of  civil  engineering  structures  involved  the  earthquake  response 
analysis  of  earth  dams  [5],  Using  the  finite  element  method  and  exclud¬ 
ing  hydrodynamic  effects,  a  number  of  questions  were  studied  concerning 
the  response  of  concrete  gravity  clams  to  earthquakes.  These  questions 
were  prompted  by  the  structural  damage ,  caused  by  the  December  1967  earth¬ 
quake,  to  Koyna  Dam. 
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The  analysis  of  hydrodynamic  pressures,  due  to  horizontal  ground  motion, 
on  rigid  dams  started  with  Westergaard1 s  pioneering  work  of  1933  [7].  Con¬ 
tinuing  with  the  assumption  of  a  rigid  dam,  more  comprehensive  analyses  of 
hydrodynamic  pressures  on  the  dam  face,  due  to  both  horizontal  and  vertical 
components  of  ground  motion,  have  been  developed  [8,12]. 

In  the  finite  element  analyses  mentioned,  hydrodynamic  effects  were  not 
considered;  whereas  in  the  studies  of  hydrodynamic  pressures,  the  dam  was 
assumed  to  be  rigid.  Additional  hydrodynamic  pressures  will  result  from 
deformations  of  the  upstream  face  of  the  dam,  and  the  structural  deforma¬ 
tions  in  turn  will  be  affected  by  the  hydrodynamic  pressures  on  the  upstream 
face.  To  break  this  closed  cycle  of  cause  and  effect,  the  problem  formula¬ 
tion  must  recognize  the  dynamic  interaction  between  the  dam  and  water. 

Finite  element  analysis  of  the  complete  dam-water  system  is  one 
possible  approach  to  including  these  interaction  effects.  Application  of 
standard  analysis  procedures  with  nodal  point  displacements  as  the 
degrees-of-freedom  was  only  partly  successful  [13].  A  different  finite 
element  formulation  of  the  complete  system,  in  which  displacements  are 
considered  as  the  unknowns  at  the  nodal  points  for  the  dam,  and  pressures 
as  the  unknowns  at  the  nodal  points  for  the  water,  has  been  applied  to 
small  problems  [14] ?  but  this  approach  appears  to  require  prohibitive 
computational  effort  for  practical  problems. 

The  more  effective  approach  is  to  treat  the  dam-water  system  as  com¬ 
posed  of  two  substructures  —  dam  and  fluid  domaiA  —  coupled  through  the 
interaction  forces  and  appropriate  continuity  conditions  at  the  face  of 
the  dam.  A  series  of  studies  [15-18]  led  to  a  general  analysis  procedure 
[19]  and  conqputer  program  [20]  for  dynamic  analyses  of  dams,  including 
dam-water  interaction.  This  approach  conveniently  permits  different 
models  to  be  used  for  the  dam  and  water.  The  dam  may  be  idealized  by 
the  finite  element  method,  which  has  the  ability  to  handle  systems  of 
arbitrary  geometry.  At  the  same  time,  the  fluid  domain  may  be  treated 
as  a  continuum,  an  approach  which  is  ideally  suited  to  the  single  geometry 
but  great  upstream  extent  of  the  impounded  water  [19] .  The  fluid  domain 
may  also  be  idealized  by  the  finite  element  method  in  conjunction  with 
infinite  elements  [21].  When  the  substructure  method  is  enployed,  along 
with  transformation  of  the  governing  equations  to  generalized  coordinates 


associated  with  vibration  modes  of  the  dam  alone,  the  analysis  procedure  is 
very  efficient,  and  little  additional  computational  effort  is  required  to  in¬ 
clude  the  hydrodynamic  effects  [19] . 

Utilizing  such  an  analysis  procedure,  it  was  shown  that  the  dam-water 
interaction  and  water  compressibility  have  a  significant  influence  on  the 
dynamic  behavior  of  concrete  gravity  dams  and  their  responses  to  earthquake 
ground  motion  [22].  Because  of  hydrodynamic  effects,  the  vertical  component 
of  ground  motion  is  more  ingwrtant  in  the  response  of  gravity  dams  than  in 
other  classes  of  structures  [18,22].  A  simplified  analysis  procedure  has 
been  developed  which  includes  the  dam-water  interaction  effects  in  the  com¬ 
putation  of  lateral  earthquake  forces  for  dam  design  [23] . 

The  effects  of  dam-foundation  interaction  can  most  simply  be  included 
in  dynamic  analysis  of  dams  by  including,  in  the  finite  element  idealization, 
foundation  rock  or  soil  above  a  rigid  horizontal  boundary.  The  response  of 
such  a  finite  element  system  to  excitation  specified  at  the  bottom,  rigid 
boundary  is  then  analyzed  by  standard  procedures.  Such  an  approach  leads  to 
enormous  computational  requirements  and  the  reliability  of  results  in  some 
cases  is  questionable.  For  sites  where  similar  materials  extend  to  large 
depths  and  there  is  no  obvious  "rigid"  boundary  such  as  a  soil-rock  interface 
the  location  of  the  rigid  boundary  introduced  in  the  analysis  is  often  quite 
arbitrary,  and  it  may  significantly  distort  the  response. 

These  difficulties  can  be  overcome  by  using  the  substructure  method  [24, 
25] ,  wherein  the  dam  and  the  foundation  rock  or  soil  region  are  considered  as 
two  substructures  of  the  complete  system.  The  dam  may  be  idealized  as  a 
finite  element  system  which  has  the  ability  to  represent  arbitrary  geometry 
and  material  properties.  The  foundation  rock  or  soil  region  may  be  idealized 
as  either  a  continuum  (a  viscoelastic  halfspace  for  example) ,  or  as  a  finite 
element  system,  whichever  is  appropriate  for  the  site.  The  half space  ideal¬ 
ization  permits  accurate  modelling  of  sites  where  similar  materials  extend  to 
large  depths.  For  sites  where  layers  of  soil  or  soft  rock  overlie  harder 
rock  at  shallow  depths,  finite  element  idealization  of  the  foundation  region 
would  be  appropriate.  The  governing  equations  for  the  two  substructures  are 
combined  by  imposing  equilibrium  and  compatibility  requirements  at  the  base 
of  the  dam.  These  equations  make  direct  use  of  free  field  ground  motion 
specified  at  the  dam- foundation  interface.  The  resulting  equations  are  trans 
formed  to  generalized  Ritz  coordinates;  the  displacements  are  expressed  as  a 


linear  combination  of  the  first  few  vibration  modes  of  an  associated  dam- 
foundation  system,  thus  leading  to  a  very  efficient  solution. 


The  preceding  review  discussed  the  dam-water  system  and  dam- foundation 
system  separately.  In  each  of  these  two  cases,  the  impounded  water  and  the 
foundation,  respectively,  modify  the  vibration  properties  of  the  dam  and  may 
significantly  affect  its  response.  However,  the  two  problems  actually  are 
coupled  and  the  results  obtained  by  separate  analyses  will,  in  general,  be 
invalid.  There  is  need,  therefore,  for  developing  techniques  for  analysis 
of  complete  dam-water-foundation  systems,  and  assessing  simultaneously,  the 
interaction  due  to  both  water  and  foundation.  Some  work  on  this  problem  has 
been  reported  in  recent  years  [26,27].  The  substructure  methods  developed  for 
separately  considering  hydrodynamic  and  foundation  interaction  effects  in  the 
response  of  dams  appear  to  be  ideally  suited  for  and  extendable  to  conqplete 
analyses  which  simultaneously  include  both  types  of  interaction. 

1.3  Scope  of  this  Report 

Chapter  2  describes  how  the  dam-water-foundation  system  and  earthquake 
ground  motion  are  idealized  in  this  study.  Also  included  in  this  chapter  are 
the  assumptions  underlying  the  procedure,  developed  in  later  chapters,  for 
analysis  of  earthquake  response  of  concrete  gravity  dams. 

For  the  reader's  convenience,  a  summary  is  given  of  the  procedures  avail¬ 
able  for  analyzing  earthquake  response  of  concrete  gravity  dams  under  restric¬ 
tive  conditions.  Chapter  3  summarizes  the  standard  finite  element  method  for 

* 

analysis  of  dams  on  rigid  foundations  with  no  water  stored  in  their  reser¬ 
voirs.  Chapter  4  summarizes  the  substructure  method  for  including  the  effects 
of  dam-water  interaction  in  the  analysis  procedure.  (This  method  treats  the 
impounded  water  and  dam  on  rigid  foundation  as  two  substructures  of  the  total 
system.)  Also  based  on  the  substructure  concept  is  the  procedure  for  includ¬ 
ing  effects  of  interaction  between  the  dam  and  its  flexible  foundation  on  the 
earthquake  response  of  dams  without  water,  summarized  in  Chapter  5. 

Because  the  substructure  concept  has  proven  to  be  effective  in  separately 
including  the  effects  of  dam-water  interaction  and  of  dam-foundation  inter¬ 
action  in  the  analysis,  it  is  extended  in  Chapter  6  to  develop  an  analysis 
procedure  simultaneously  including  various  effects  of  the  water  and  the 
foundation.  These  include  effects  arising  from  interaction  between  the  dam 
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and  foundation,  dam  and  water,  water  and  foundation,  and  interaction  among 
all  three  substructures  —  dam,  water,  and  foundation.  The  computer  program 
developed  to  implement  the  procedure  for  analysis  of  dam-water-foundation 
systems  is  briefly  described  in  this  chapter;  the  user's  guide  and  program 
listing  is  included  in  Appendices  A  and  B. 

Using  this  computer  program,  results  of  several  preliminary  analyses 
are  presented  in  Chapter  7  with  the  aim  of  defining  the  important  system  and 
analysis  parameters  for  the  subsequent  study  of  dynamic  response  behavior  of 
dams.  The  behavioral  study  is  in  two  parts,  separated  into  Chapters  8  and  9. 

The  response  of  idealized  dam  cross-sections  to  harmonic  horizontal  or 
vertical  ground  motion  is  presented  in  Chapter  8  for  a  range  of  the  important 
systems  parameters  characterizing  the  properties  of  the  dam,  foundation  rock, 
and  impounded  water.  With  the  aid  of  these  results,  the  separate  effects  of 
dam-water  interaction  and  dam-foundation  interaction,  and  the  combined  effects 
of  the  two  sources  of  interaction  on  dynamic  response  of  dams  are  investigated. 

Chapter  9  presents  the  responses  of  Pine  Flat  Dam  to  the  S69E  component 
of  the  Taft  ground  motion  only;  and  to  the  S69E  and  vertical  components  act¬ 
ing  simultaneously.  For  each  of  these  excitations,  the  response  of  the  dam 
is  analyzed  four  times  corresponding  to  the  following  four  sets  of  assump¬ 
tions:  (1)  rigid  foundation,  hydrodynamic  effects  excluded;  (2)  rigid 

foundation,  hydrodynamic  effects  included;  (3)  flexible  foundation,  hydro- 
dynamic  effect  excluded;  and  (4)  flexible  foundation,  hydrodynamic  effects 
included.  These  results  provide  insight  into  the  effects  of  dam-water  and 
dam- foundation  interaction,  considered  separately  or  together,  in  the  earth¬ 
quake  response  of  dams. 

Chapter  10  briefly  summarizes  the  significant  conclusions  which  may  be 
drawn  from  this  investigation  of  dam-water  and  dam-foundation  interaction 
effects  in  earthquake  response  of  concrete  gravity  dams. 

1.4  Acknowledgements 

This  study  was  initially  supported  by  the  Office  of  the  Chief  of  Engi¬ 
neers,  Department  of  the  Army,  Washington,  D.C.,  under  Contract  No.  DACW73- 
71-C-0051.  Completion  of  this  study  was  made  possible  by  Grants  ATA74-20554 
and  ENV76-80073  from  the  National  Science  Foundation.  The  authors  are  grateful 
for  the  support  from  both  agencies. 
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2.  SYSTEM  AMD  GROUND  MOTION 


2.1  Dam-Water-Foundation  System 

Vibration  teste  on  Pine  Flat  Dam  [28]  indicate  that  at  small  vibration 
amplitudes  a  concrete  gravity  dam  will  behave  like  a  solid  even  though  there 
is  some  slippage  between  monoliths.  Thus,  at  the  beginning  of  an  earthquake, 
the  behavior  of  a  dam  can  be  best  described  by  a  three-dimensional  model. 
However,  at  large  amplitudes  of  motion,  the  inertia  forces  are  much  larger 
than  the  shear  forces  that  can  be  transmitted  across  joints  between  monoliths. 
The  monoliths  slip  and  tend  to  vibrate  independently,  as  evidenced  by  the 
spalled  concrete  and  increased  water  leakage  at  the  joints  of  Koyna  Dam  during 
the  Koyna  earthquake  of  December  11,  1967.  Consequently,  two-dimensional 
models  of  individual  monoliths  appear  to  be  more  appropriate  than  three- 
dimensional  models  for  predicting  the  response  of  concrete  gravity  dams  to 
the  strong  phase  of  intense  ground  motion.  However,  a  model  more  complicated 
than  either  of  these  two  models  will  be  necessary  to  describe  the  behavior  of 
a  concrete  gravity  dam  through  the  complete  amplitude  range.  . 

Because  the  dimensions  and  dynamic  properties  of  the  various  monoliths 
differ,  the  effects  of  the  dam  on  deformations  and  stresses  in  the  founda¬ 
tion  will  vary  along  the  length  of  the  dam.  Thus  even  with  two-dimensional 
models  for  the  dam,  a  three-dimensional  model  would  seem  necessary  for  the 
foundation.  If  the  dam  were  to  behave  as  a  solid  without  slippage  between 
monoliths,  and  all  of  its  properties  and  the  ground  motion  did  not  vary  along 
the  length,  it  would  be  appropriate  to  assume  the  dam  as  well  as  the  founda¬ 
tion  to  be  in  plane  strain.  A  plane  stress  model  which  applies  to  a  thin 
sheet-like  body  seems  obviously  inappropriate  for  a  continuum  foundation. 

Although  the  basic  concepts  underlying  the  analysis  procedure  presented 
in  this  work  are  applicable  under  more  general  conditions,  the  procedure  is 
specifically  developed  for  two-dimensional  systems;  thus  its  application  is 
restricted  to  systems  in  generalized  plane  stress  or  plane  strain.  Although 
neither  of  the  two  models  are  strictly  applicable,  the  former  is  better  for 
the  dam  and  the  latter  for  the  foundation.  However,  in  order  to  define  the 
dam- foundation  system  on  a  consistent  basis,  the  same  model  should  be 
employed  for  both  substructures.  Results  from  the  two  models  are  compared 
in  Sec.  7.3  to  provide  a  basis  for  choosing  one  for  parameter  studies. 


1 

A  cross-section  of  the  system  considered  is  shown  in  Fig.  2.1.  The 
system  consists  of  a  concrete  gravity  dam  supported  on  the  horizontal  sur¬ 
face  of  a  viscoelastic  half  plane  and  impounding  a  reservoir  of  water.  The 
system  is  analyzed  under  the  assumption  of  linear  behavior  for  the  concrete, 
foundation  soil  or  rock,  and  water.  The  dam  is  idealized  as  a  two-dimen¬ 
sional  finite  element  system,  thus  making  it  possible  to  consider  arbitrary 
geometry  and  variation  of  material  properties.  However,  certain  restric¬ 
tions  on  the  geometry  are  imposed  to  permit  solutions  for  the  foundation  and 
fluid  domains  treated  as  continua.  For  the  purposes  of  determining  hydro- 
dynamic  effects,  and  only  for  this  purpose,  the  upstream  face  of  the  dam  is 
assumed  to  be  vertical.  This  is  reasonable  for  an  actual  concrete  gravity 
dam,  because  the  upstream  face  is  vertical  or  almost  vertical  for  most  of 
its  height,  and  the  hydrodynamic  pressures  on  the  dam  face  are  insensitive 

< 

to  small  departures  of  the  face  slope  from  the  vertical.  For  the  purpose  of 
including  structure-foundation  interaction  effects,  the  foundation  surface 
is  assumed  to  be  horizontal;  thus  the  base  of  the  dam  as  well  as  the  reser¬ 
voir  bottom  is  assumed  to  be  horizontal.  An  actual  system  can  usually  be 
idealized  to  conform  to  this  assumption. 

i 

2.2  Ground  Motion 

The  excitation  for  the  two-dimensional  dam-foundation-water  system  is 
defined  by  the  two  components  of  free-field  ground  acceleration  in  the  plane 
of  a  cross-section  of  the  dam:  the  horizontal  component  transverse  to  the 

dam  axis,  ax(t),  and  the  vertical  component,  a^(t). 
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3.  ANALYSIS  OF  GRAVITY  DAMS 


3.1  Introductory  Note 

The  standard  finite  element  method  for  analysis  of  gravity  dams  on  rigid 
foundations  with  no  water  stored  in  the  reservoir  is  summarized  in  this 
chapter.  In  later  chapters,  the  analysis  will  be  extended  to  include  effects 
of  impounded  water  and  foundation  flexibility. 

3. 2  Governing  Equations 

Consider  a  monolith  of  a  concrete  gravity  dam  on  a  rigid  foundation, 
with  no  water  in  the  reservoir,  subjected  to  earthquake  ground  motion  which 
does  not  vary  across  the  base  of  the  monolith.  The  equations  of  motion  for 
such  a  dam  monolith  idealized  as  a  planar,  two-dimensional  finite  element 
system  ares 

mr  +  cr  +  kr  =  -mlx  aX(t)  -  mly  ay(t)  (3.1) 

—  —  g  —  g 

In  Eq.  3.1,  m,  k  and  c  are  the  mass,  stiffness,  and  viscous  damping  matrices 
for  the  finite  element  system;  r  is  the  vector  of  nodal  point  displacements, 
relative  to  the  free-field  displacement: 


/r*  ry  r  r* 
^  1  1  2  2 


x  y 
r  rJ 
n  n 


N 


where  rX  and  ry  are  the  x-  and  y-  components  of  displacement  of  nodal  point 
n  n 

n  and  the  number  of  nodal  points  above  the  base  is  N;  r  and  r  are,  respec¬ 
tively,  the  nodal  point  velocity  and  acceleration  vectors. 


T 

{lx}  ■  <1  0  1  0 . 1  0 . 1  0> 

T 

{ly}  -  <0  1  0  1 . 0  1 . 0  1> 


ax(t)  and  ay(t)  are  the  x  (horizontal)  and  y  (vertical)  components  of  the 
9  9 

free-field  ground  acceleration. 

The  stiffness  and  mass  matrices  of  the  structure  are  obtained  from  the 
corresponding  matrices  for  the  individual  finite  elements  by  direct  assembly 


procedures.  The  element  stiffness  matrices  are  derived  using  quadratic  in¬ 
terpolation  function  for  displacements,  whereas  the  element  mass  matrices 
are  based  on  a  lumped  mass  approximation. 

Energy  dissipation  in  structures,  even  in  the  linear  range  of  vibration, 
is  due  to  various  complicated  phenomena.  Because  it  is  not  possible  to 
mathematically  describe  each  of  these  sources  of  energy  dissipation,  it  is 
customary  to  define  damping  in  terms  of  damping  ratios  in  the  natural  inodes 
of  vibration.  Damping  ratios  for  similar  structures  obtained  from  both  the 
analysis  of  harmonic  vibration  tests,  and  responses  recorded  during  earth¬ 
quakes,  are  used  as  a  basis  for  assigning  the  modal  damping  ratios.  The 
modal  damping  ratios  provide  a  complete  description  of  damping  properties 
for  purposes  of  linear  analysis  and  the  damping  matrix  need  not  be  defined 
explicitly. 

Traditionally,  analysis  of  dynamic  repsonse  of  structures  has  been 
carried  out  directly  in  the  time  domain.  For  such  analyses,  viscous  damping 
is  the  most  convenient  representation  of  energy  dissipation  in  the  struc¬ 
ture.  When  effects  of  dam-water  interaction  and  dam-foundation  interaction 
are  included,  the  analysis  is  most  effectively  carried  out  by  a  substructure 
method  [3,4].  As  will  be  seen  in  Chapters  4,5  and  6,  the  substructure 
method  is  best  formulated  in  the  frequency  domain.  In  such  a  formulation, 
constant  hysteretic  damping  is  a  preferable  representation  —  for  conceptual 
as  well  as  computational  reasons  —  especially  when  structure-soil  inter¬ 
action  effects  are  included  [5]. 

The  Fourier  transform  of  Eq.  2.1  provides  the  governing  equation  in  the 
frequency  domain 

[-d)2m  +  iuic  +  k]  r  («)  *  -mlX  ax(w)  -  mly  ay(u)) 

-  -  -  -  g  -  g 

A 

in  which  the  Fourier  transform  of  f(t)  is  denoted  by  f (u) . 
hysteretic  damping,  Eq.  2.2  becomes 

[-w2m  +  (1  +  ih)k]  r(w)  -  -mlx  ax(oi)  -  mly  ay(w)  (3.3) 

-  -  g  g 

where  f|  is  the  constant  hysteretic  damping  factor. 

Damping  ratios  as  obtained  from  forced  vibration  tests  on  dams  are 


(3.2) 

For  constant 


essentially  independent  of  the  mode  number  [28].  If  £  is  the  damping  ratio 
appropriate  for  all  the  natural  modes  of  vibration  of  the  dam,  n  ■  2£  would 
be  an  appropriate  value  for  the  constant  hysteretic  damping  factor  [29] . 

With  the  damping  coefficient  so  related,  essentially  the  same  response  will 
be  obtained  for  a  lightly  damped  system  with  either  damping  mechanism. 

3.3  Earthquake  Response  Analysis 

3.3.1  Time  Domain  Analysis 

The  earthquake  response  of  a  dam  is  obtained  by  solving  the  equations 
of  motion  (Eq.  3.1).  These  equations  of  motion  in  nodal  point  coordinates 
may  be  solved  either  directly  or  after  transformation  to  modal  coordinates 
[30].  The  latter  approach,  commonly  known  as  the  mode  superpostion  method, 
is  applicable  if  the  response  is  within  the  linear  range.  This  method  is 
advantageous  for  calculating  the  earthquake  response  of  many  types  of  struc¬ 
tures,  because  the  response  is  essentially  due  to  the  first  few  modes  of 
vibration. 

The  first  step  in  the  mode  superposition  analysis  procedure  is  to  ob¬ 
tain  the  lower  few  natural  frequencies  and  mode  shapes  of  vibration  of  the 
dam  by  solving  the  eigenvalue  problem: 

k  $  »  u2  m  0  (3.4) 


The  equations  of  motion  (Eq.  3.1)  are  uncoupled  by  the  transformation 

2N 

r (t)  -  Z  6  Y  (t)  (3.5) 

-  ,  -n  n 

n*l 

provided  the  damping  matrix  £  satisfies  certain  restrictions  [6] .  In  Eq. 
3.5,  Y  (t)  is  the  generalized  coordinate  and  $  the  mode  shape  for  the  nth 
natural  mode  of  vibration.  The  uncoupled  equation  for  the  n  mode  of  vi¬ 
bration  is 


MY  +CY  +KY  «  -LX  aX(t)  -  Ly  ay(t)  (3.6) 

nn  nn  nn  ng  ng 


I 


where 


£  (2M  bi  ) ,  K 
*n  n  n  n 


2 

u)  M  ,  u  and  £  are  the  natural 
n  n  n  n 
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th  x  T  x 

circular  frequency  and  damping  ratio  for  the  n  mode,  L  *  ml  ,  and 

Ly  -  *T  »ly. 

n  -^n  — 


Equation  3.6  may  be  solved  for  Yn(t)  by  a  step-by-step  integration 
method  [30].  After  computation  has  been  repeated  for  all  modes,  the  nodal 
displacement  vector  r(t)  can  be  obtained  from  Eq.  3.5.  In  practice,  it  is 
generally  sufficient  to  solve  the  equations  of  motion  only  for  the  lower  few 
modes,  because  the  contributions  of  the  higher  modes  to  the  total  response 
are  small. 

The  stresses  a^(t)  *n  element  p  at  any  instant  of  time  are  re¬ 

lated  to  the  nodal  displacement  r^(t)  for  that  element  by 

0  (t)  «T  r  (t)  (3.7) 

-P  -p  “P 


where  the  stress  transformation  matrix  T  is  based  on  the  interpolation 
functions  for  the  el men t  as  well  as  its  elastic  properties.  The  stresses 
throughout  the  dam  at  any  instant  of  time  are  determined  from  the  nodal  point 
displacements  by  application  of  the  above  transformation  to  each  finite 
element . 

The  initial  stresses,  before  the  earthquake,  should  be  added  to  the 
stresses  due  to  earthquake  excitation,  determined  by  the  procedures  pre¬ 
sented  in  the  preceding  sections,  to  obtain  the  total  stresses  in  the  dam. 
Excluding  the  temperature  and  creep  stresses  in  concrete,  only  the  gravity 
loads  need  to  be  considered.  The  equations  of  static  equilibrium  are  for¬ 
mulated: 


K  r  =  R  (3.8) 

where  R  is  the  vector  of  static  loads  due  to  the  weight  of  the  dam  and  hydro¬ 
static  pressures.  Solution  of  these  algebraic  equations  results  in  the  dis¬ 
placement  vector  r.  Static  stresses  are  then  determined  from  the  displace¬ 
ments  by  applying  Eq.  3.7  to  each  finite  element. 

Analytical  predictions,  based  on  the  procedure  summarized  above,  for  the 
perfomance  of  Koyna  Dam  correlated  well  with  the  damage  experienced  by  the 
dam  during  the  Koyna  earthquake  of  December  11,  1976  [2]  . 


3*3.2  Frequency  Domain  Analysis 

An  alternative  approach  to  solving  the  modal  equations  of  motion  (Eq. 
3.6)  exists  in  which,  instead  of  step-by-step  integration  in  the  time 
domain,  the  complex  frequency  responses  are  superposed  in  the  frequency 
domain  utilizing  the  Fourier  integral  [30] .  The  complex  frequency  response 
function  z(w)  for  a  response  quantity  z(t)  has  the  property  that,  when  the 
excitation  is  the  real  part  of  e  ,  then  the  response  is  the  real  part  of 
2(u»)eia,t. 

The  response  in  the  nth  mode  of  vibration  due  to  the  horizontal  ground 
motion  will  be  denoted  by  Y^,  and  that  due  to  the  vertical  ground  motion  by 
Yn*  1116  resP°nse  to  excitation  a^(t)  «  e1U)t  can  then  be  expressed  as 

y£  a  y£  gikfc  where  Z  *  x  or  y. 
n  n 

Substitution  in  Eq.  3.6  leads  to 

<-u>2Mn  +  itoCn  +  Kn)  Y*(u»  =  -  L*;  Z  -  x,y  (3.9) 

from  which 

_LZ 

*„<“)  *  5 - - -  1  Z  -  x,y  (3.10) 

-(*)  M  +  i(t)C  +  K 
n  n  n 

If  the  system  has  constant  hysteretic  damping  instead  of  viscous  danc¬ 
ing,  starting  from  Eq.  3.3,  it  can  be  shown  that  Eqs.  3.9  and  3.10  become 

[-<o2M  +  (1  +  ir,)K  J  (w)  =  -L^j  Z  =  x,y  (3.11) 

n  n  n  n 

and 

-Lz 

y^(oj)  -  — - 2 -  .  l  ,  x , y  (3.12) 

-OJ  H  +  (1  +  in)K 
n  n 

The  response  to  arbitrary  ground  motion  is  the  superpostion  of 
responses  to  individual  harmonic  components  of  the  excitation  through  the 
Fourier  integral: 

9  00 

V*>  "  /  Yn(w)  Vw)  1  “  *'*  (3-13) 
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in  which  A  (u>)  is  the  Fourier  transform  of  a  (t)  : 
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A^(uj)  «  /  a^(t)e  dt;  Z  =  x,y 


(3.14) 


where  d  is  the  duration  of  the  ground  motion. 

The  response  in  the  nfc^  vibration  mode  due  to  simultaneous  action  of  the 
horizontal  and  vertical  components  of  ground  motion  is 

Y  (t)  =  YX(t)  +  Yy(t)  (3.15) 

n  n  n 

Repeating  this  procedure  for  all  the  necessary  values  of  n,  the  displacement 
responses  may  be  obtained  by  superpostion  of  the  modal  responses  (Eq.  3.5) 
and  the  stress  responses  by  calculating  the  stresses  associated  with  these 
displacements  (Eq.  3.7).  The  initial  stresses  are  added  to  determine  the 
total  stresses. 

Until  the  development  of  the  Fast  Fourier  Transform  (FFT)  algorithm 
[31],  numerical  evaluation  of  integrals  such  as  those  in  Eqs.  3.13  and  3.14 
required  prohibitive  amounts  of  computer  time  and  the  errors  in  the  results 
could  not  be  predicted  accurately.  As  a  result,  step-by-step  integration  in 
the  time  domain  has  conventionally  been  used  for  response  analysis.  With 
the  FFT  algorithm,  integrals  of  Eqs.  3.13  and  3.14  can  be  evaluated  accu¬ 
rately  and  efficiently.  As  a  result,  the  frequency  domain  approach  can  now 
be  employed  advantageously  for  analysis  of  dynamic  response  of  structures. 

It  provides  an  alternative  approach  for  analysis  of  dams  without  water.  As 
will  be  seen  in  Chapters  4,5,  and  6,  the  analysis  is  best  formulated  by  the 
frequency  domain  approach  when  dam-water  and  dam-foundation  interaction 
effects  are  included. 
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4.  ANALYSIS  OF  GRAVITY  DAMS  INCLUDING  DAM-WATER  INTERACTION 

4.1  Introductory  Note 

Summarized  in  this  chapter  is  th><  procedure  presented  earlier  [19] ,  for 
including  the  effects  of  dam-water  interaction  in  earthquake  response  anal¬ 
yses  of  concrete  gravity  dams.  The  concept  under lying  this  procedure  is  to 
treat  the  impounded  water  and  the  dam  as  two  substructures  of  the  total 
system.  With  this  approach,  the  dam  may  be  idealized  by  the  finite  element 
method,  which  has  the  ability  to  represent  systems  of  arbitrary  geometry  and 
material  properties  (see  Chapter  3) .  At  the  same  time,  the  impounded  water 
may  be  treated  as  a  continuum,  an  approach  which  is  ideally  suited  to  the 
simple  idealized  geometry  but  great  upstream  extent  of  the  impounded  water. 
The  equations  of  motion  for  the  two  substructures  are  coupled  by  including 
the  forces  of  interaction  between  the  water  and  the  face  of  the  dam. 

4.2  Governing  Equations 
4.2.1  Substructure  1 :  Dam 

A  gravity  dam  supported  on  a  rigid  foundation  and  storing  water  to  a 
given  depth  is  shown  in  Fig.  4.1.  Including  hydrodynamic  effects,  the  equa¬ 
tions  of  motion  are,  for  a  dam  monolith  idealized  as  a  planar  two-dimensional 
finite  element  system  (Chapter  3) ,  subjected  to  horizontal  and  vertical  com¬ 
ponents  of  ground  motion: 

mr  +  cr  +  kr  =  -mlx  ax(t)  -  mly  ay(t)  +  R,  (t)  (4.1) 

- - -  g  9  -n 

This  is  identical  to  Eq.  3.1  except  for  the  inclusion  of  the  hydrodynamic 
loads  R^(t)  on  the  upstream  e  of  the  dam. 

Whereas  Eq.  4.*  is  valid  for  dams  of  arbitrary  geometry,  it  is  assumed 
that  for  the  purpose  of  defining  hydrodynamic  effects ,  and  only  for  this  pur¬ 
pose,  the  upstream  face  of  the  dam  is  vertical.  Considering  that  the  up¬ 
stream  face  of  concrete  gravity  dams  is  vertical  or  almost  vertical  for 
most  of  the  height,  this  is  a  reasonable  assumption.  With  this  assumption 
and  because  the  hydrodynamic  loads  act  only  on  the  upstream  face,  only  those 
loads  in  R^(t)  that  correspond  to  the  x-degrees  of  freedom  of  the  nodal 
points  on  the  upstream  face  are  non-zero.  The  sub-vector  of  the  non-zero 

xf 

loads  is  denoted  by  R£  (t) .  The  superscript  "xf "  indicates  that  only  the 
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x-component  of  the  loads  acting  on  the  face  of  the  dam  are  included  in  the 
vector. 

The  procedure  for  analyzing  the  response  of  dams  without  including  the 
effects  of  water,  presented  in  Chapter  3,  utilized  the  mode- superpostion 
concept.  Although  classical  natural  inodes  of  vibration  do  not  exist  for 
the  dam  when  hydrodynamic  effects  are  included,  expansion  of  displacements, 
including  hydrodynamic  effects,  in  terms  of  natural  modes  of  vibration  with¬ 
out  water  can  still  be  used  advantageously.  Thus, 

J 

r  (t)  *  Z  <p  Y  (t)  (4.2) 

,  -n  n 
n=l 

The  vectors  n  =  1,2...2N  are  linearly  independent  and  span  the  vector 
space  of  dimension  2N.  Thus,  the  expansion  of  Eq.  4.2  is  exact  if  the  con¬ 
tributions  of  all  the  2N  natural  modes  of  vibration  are  included.  However, 
even  when  hydrodynamic  effects  are  included,  the  contributions  of  the  lower 
modes  of  vibration  are  expected  to  be  more  significant  and  relatively  few 
terms,  e.g.  J  «  2N,  would  suffice  in  Eq.  4.2. 


Substituting  Eq.  4.2  into  Eq.  4.1  and  utilizing  the  orthogonality 
property  of  mode  shapes,  the  equation  governing  the  generalized  displace¬ 
ment  Y  (t)  in  the  nth  mode  is 
n 


MY  +  C  Y 
n  n  n  n 


K  Y 
n  n 


lx  xr*.\ 

-L  a  (t) 
n  g 


-  Ly  ay(t) 
n  g 


+  Rhn 


(t) 


(4.3) 


This  is  the  same  as  Eq.  3.6,  except  for  the  additional  term  R^ft)  =  ■ 

(j^R^t)  *  the  vector1  consists  of  the  x-components 

of  the  displacements  of  the  nodal  points  on  the  upstream  face  of  the  dam. 

4.2.2  Substructure  2s  Fluid  Domain 

The  motion  of  the  dam-water  system  is  considered  to  be  two-dimensional, 
i.e.,  it  is  the  same  for  any  vertical  plane  perpendicular  to  the  axis  of  the 
dam.  Assuming  water  to  be  linearly  compressible  and  neglecting  its  internal 
viscosity,  the  small  amplitude,  irrotational  motion  of  the  fluid  is  governed 
by  the  two-dimensional  wave  equation 


3x  3y* 


_1_  3_£ 
c2  3t2 


(4.4) 
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in  which  p(x,y,t)  is  the  hydrodynamic  pressure  (in  excess  of  the  hydrostatic 
pressure)  and  C  is  the  velocity  of  sound  in  water.  The  hydrodynamic  pres¬ 
sure  acting  on  the  upstream  face  of  the  dam  is  pc(y,t)  =  p(o,y,t). 

4.3  Response  Analysis  for  Harmonic  Horizontal  Ground  Motion 

Specializing  the  equations  of  motion  of  the  dam  in  terms  of  generalized 
modal  displacements  (Eq.  4.3)  for  the  case  when  the  ground  motion  acts  only 
in  the  horizontal  (x)  direction  leads  to 


M  YX  +  C  YX  +  K  YX  *  -LX  aX(t)  +  R?  (t) 
nn  nn  nn  ng  nn 


(4.5) 


where  L*  *  ml*  and  (t)  *  {^xf}T  R^f  (t)  .  The  vector  R^f  (t)  consists 
of  forces,  associated  with  hydrodynamic  pressures  acting  in  the  horizontal 
direction,  at  the  nodal  points  on  the  vertical  upstream  face  of  the  dam. 
These  hydrodynamic  pressures  are  governed  by  the  two-dimensional  wave  equa¬ 
tion  (Eq.  4.4)  and  the  following  boundary  conditions: 

•Pressure  at  the  free  surface  is  zero,  implying  that  the  effects  of 
waves  at  the  free  surface  are  ignored; 

•Vertical  motion  at  the  base  of  the  reservoir  is  zero;  and 


•Horizontal  component  of  motion  of  the  fluid  boundary  x  *  0  is  the 
same  as  the  horizontal  motion  of  the  upstream  face  of  the  dam. 

The  responses  to  harmonic  ground  acceleration  ax(t)  =  eia>t  cam  be  ex- 

9 


pressed  as  follows: 

Generalized  Displacements,  Y*  (t)  =  Y x(w)eia,t 

n  n 

Generalized  Accelerations,  ¥x(t)  = 

n  n 


(4.6a) 

(4.6b) 


Horizontal  accelerations  at  the  upstream  face, 

J 


ii(t)  *  {i  +  E  $xf  Yx(u»}e 
j-1  3  3 


iu)t 


(4.6c) 


in  which  1  is  a  vector  with  all  elements  equal  to  unity. 

Hydrodynamic  pressures,  px(y,t)  »  px(y,u))e^wt 

c  c 


(4.6d) 


xf 


Strictly  for  purposes  of  notational  convenience,  the  vector  ^  is  replaced 
by  its  continuous  analog  function  ^  (y) .  The  total  horizontal  acceleration 
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of  the  fluid  boundary  at  x  ■  0,  the  upstream  face  of  the  dam,  may  now  be 
represented  by  the  function 


u(0,y,t)  -  {l  +  Z  <My)  ir.(uj)}e 
j-1  3  3 


iojt 


(4.7) 


With  reference  to  Fig.  4.1,  the  boundary  conditons  for  the  wave  equation 
are 


p(x,H,t)  =  0 
|£  (x,0,t)  -  0 

^  (0  ,y,t)  *  -  £{1  +  Z  <J> .  (y)  Y^wjJe10* 


9x 


j-1 


(4.8) 


Because  the  governing  equation  as  well  as  the  boundary  conditions  are 
linear,  the  principal  of  superpostion  applies.  The  complex  frequency  re- 
sponse  function  for  (y,t)  can  therefore  be  expressed  as 


J  - 


P*(Y.w)  »  Pq  (y,oi)  +  Z  Y*(w)  p.  (y,u>) 

j*l  3  3 


(4.9) 


x  —x  icat 

where  pQ(y,t)  *  PQ(y#w)e  is  the  solution  of  the  wave  equation  at  x  *  0 
(Eq.  4.4)  for  the  following  boundary  conditions: 


p(x,H,t)  =  0 
(x,0,t)  -  0 
3p  ..  w  iut 

3x  ,0'y't)  *  -  5  e 


(4.10) 


and  Pj(y»t)  *  Pj (y,«)e^U)t  is  the  solution  of  the  wave  equation  at  x  *  0  for 
the  following  boundary  condtions: 


p (x,H, t) 

5? 


0 

0 


3£ 

3x 


(o,y,t) 


w  .  ,  ,  iwt 

g  (y)e 


(4.11) 


Solutions  of  Eq.  4.4  for  the  above  two  sets  of  boundary  conditons  are 
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presented  in  Eqs.  4.12  and  4.13 


Pq^.W) 


4w  £  _ (-1;  _ 

*-l  (.21-1) /\Z  -u2/c2 


cos  X^y 


(4.12) 


Pj(y.W) 


"  9"  Jl  /ij  -o2/C2' 


cos  X 


(4.13) 


n 

in  which  X»  =  (2i~l)  TT/2H  and  I..  ■  /  <(>,(y)cos  X*y  dy.  The  complex  frequency 
.  -x  J*-  o  3  -c 

response  function  Pg(y,to)  is  for  the  hydrodynamic  pressures  on  the  dam  when 
the  excitation  to  the  fluid  domain  is  the  horizontal  acceleration  of  the 


ground  and  the  dam  is  rigid.  The  corresponding  function  is  p*(y,<u)  when  the 
excitation  is  the  acceleration  of  the  dam  in  the  j 
there  is  no  motion  of  the  base  of  the  reservoir. 


.th 


mode  of  vibration  and 


The  integrals  I  »  are  computed  as  a  summation  using  the  elements  of  the 
xf 

vector  instead  of  the  continuous  function  (y)  which  was  introduced 

only  for  notational  convenience.  The  nodal  force  vectors  are  RX(t)  * 
iwt  -  1  —  —  ~° 


iwt 


RgfwJe""'''  and  Rj(t)  »  R^  (u))e'*’wl‘/  where  R*(u))  and  R*(u>)  are  the  static  equiv¬ 
alents  of  the  corresponding  pressure  functions  pX(y»w)  and  p_,  (y ,co)  ,  respec¬ 
tively.  They  may  be  computed  directly  by  using  the  principle  of  virtual 
work  with  the  displacements  between  nodal  points  defined  by  the  finite 

element  interpolation  functions.  The  complex  frequency  response  function 

xf 

for  the  total  force  vector  R^  (t)  is.  from  Eq.  4.9, 


=x. 


5f<u> 


j  _ 

RX(U)  +  l  ?*(U))  R  (U» 
j-1  3  "3 


(4.14) 


The  hydrodynamic  forces  on  the  dam  have  thus  been  expressed  in  terms  of  the 
unknown  generalized  coordinate  responses  Y*(iu)  . 

With  the  aid  of  Eq.  4.14  and  substituting  ¥x(w)  ■  -w2  Y*(o>),  Bq.  4.5 

X  ^  {  Mf  ^ 

can  be  written  as  follows  for  the  excitation  a  (t)  *  e  : 

9 


[-A  *  fc*n  .  K„]y>)  .  -1*  .  ^  ?*<«»)<“>} 


(4.15) 
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The  set  of  equations  4.15  for  n  »  1,2  ...  J  may  be  rearranged  and  expressed 
in  matrix  form  as 


"S11(W>  S12(U)  •  ’  *  *  S1J(U)' 

YX(co) 

3 

>Th 

S21  (<o)  s22(co)  ....  S2J(CO) 

YX«o) 

LX(co) 

m 

• 

• 

• 

• 

_SJl(al)  SJ2(U)  •  •  •  •  SJJ(W)^ 

yJ(CO) 

t  < 

LX«0) 

>  . 

in  which 

S  .  (w)  ■  w2(4>!!f}T  R.  (to)  for  n  ^  j 
nj  “j 

S  (CO)  -  (-<02M  +  i<OC  +  K  ]  +  (02{<jiXf}T  R  (to) 

nn  n  n  n  -n  -n 

L  (co)  -  +  id  }  Rrt(co) 

n  n  zn  — 0 


(4.16) 


J  (4.17) 
J 


The  coefficient  matrix  S(co)  in  Eq.  4.17  is  a  frequency  dependent  matrix 
which  relates  the  generalized  displacement  vector  YX(co)  with  the  correspon¬ 
ding  generalized  load  vector  LX(co).  Unlike  classical  modal  analysis,  the 
matrix  S(co)  is  not  diagonal  because  the  vectors  $n  are  not  the  natural  vibra¬ 
tion  modes  of  the  dam-fluid  system;  they  simply  are  the  natural  vibration 
modes  of  the  dam  without  water.  It  can  be  shown  that  S  (co)  is  a  symmetric  matrix. 

With  one  change,  Eqs.  4.16  and  4.17  are  valid  even  if  dancing  for  the 
dam  is  idealized  as  constant  hysteretic  damping  instead  of  viscous  damping. 
The  expression  for  the  diagonal  elements  of  S  (co)  changes  to 

S  (co)  -  [-<o2M  +  (1  +  in)  K  ]  +  co2  UXf}T  R  (co)  (4.18) 

nn  n  n  *n  -n 


Solution  of  Eq.  4.16  for  all  values  of  excitation  frequency  would  give 
the  complete  frequency  response  functions  for  all  the  generalized  coordin¬ 
ates  YX(co),  n  -  1,2  ...  J.  The  frequency  responses  for  generalized  accel- 
n 

erations  may  be  obtained  from  the  generalized  displacements  as 

Vx (to)  -  -to2YX(co)  (4.19) 

n  n 
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The  complex  frequency  response  functions  for  the  vectors  of  nodal  displace¬ 
ments  and  accelerations  relative  to  the  ground  motion  are  obtained  by  super¬ 
position  of  the  contributions  of  the  various  modes  of  vibration  (Eq.  4.2) : 


?(<»>) 


J 

=  Z 

n=l 


Y*(u>) 

n 


(4.20) 


rX{uj)  =  -u)2 


J 

Z 

n*l 


tn  *>> 


(4.21) 


4.4  Response  Analysis  for  Harmonic  Vertical  Ground  Motion 

Specializing  the  equations  of  motion  of  the  dam  in  terms  of  general¬ 
ized  modal  displacements  (Eq.  4.3)  for  the  case  when  the  ground  motion  acts 
only  in  the  vertical  (y)  direction  leads  to: 


•*v 
M  Yy 
n  n 


C  Yy  +  K  Yy 
n  n  n  n 


<  ■;<*>  * 


(4.22) 


where  Ly  *  ^  a  ly  and  R^^t)  *  {$*f}T  R^f (t) .  The  vector  (t)  consiste 
of  forces,  associated  with  hydrodynamic  pressures  acting  in  the  horizontal 
direction,  at  the  nodal  points  on  the  vertical  upstream  face  of  the  dam. 
These  hydrodynamic  pressures  are  governed  by  the  two-dimensional  wave 
equation  (Eq.  4.4)  together  with  the  following  boundary  conditions  (see  Fig. 
4.1)  : 


•Pressure  at  the  free  surface  is  zero,  implying  that  the  effects 
of  waves  at  the  free  surface  are  ignored; 

d 

•Vertical  component  of  motion  of  the  fluid  boundary  y  *  0,  the  base 
of  the  reservoir,  is  prescribed  by  the  vertical  component  of  ground 
acceleration;  and 


•Horizontal  component  of  motion  of  the  fluid  boundary  x  *  0  is  the 
same  as  the  horizontal  motion  of  the  upstream  face  of  the  dam. 

The  responses  to  harmonic  ground  acceleration  ay(t)  *  eiWt  can  be  ex¬ 
pressed  as  follows: 

iu>t 


Generalized  displacements,  Y^(t) 


Yy  (u)e 
n 


(4.23a) 
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(4.23b) 


Generalized  accelerations,  ¥y(t)  *  Yy  (u)  eia) 

n  n 

Horizontal  accelerations  at  the  upstream  face. 


Hydrodynamic  pressures,  py(y,t)  *  py  (y,u>)eiait 

c  c 


(4.23c) 

(4.23d) 


Replacing  the  vector  ^  by  its  continuous  analog  function,  (y) ,  as  before, 
the  horizontal  acceleration  of  the  fluid  boundary  at  x  —  0,  the  upstream 
face  of  the  dam,  is  represented  by  the  function: 


u(0,y,t)  =  {e  <f>.  (y)Yy(w)}siwt 
j-1  3  3 


Then  the  boundary  conditions  for  the  wave  equation  are 


P(x,H,t) 


(x,  0,  t) 
fJ(o,y,t) 


3y 

3 


w  iut 
-  —  e 

g 


Yy«u)}eiuJt 


9  j-l'3 


(4.24) 


(4.25) 


Because  the  governing  equations  as  well  as  boundary  conditions  are  linear, 
the  principle  of  superpositions  applies.  The  complex  frequency  response 
function  for  py(y,t)  can  therefore  be  expressed  as 


p£(y,w) 


J  _ 

Pn<Y*w>  +  E  *?«*>)  P^(y»w) 
j-1  3  3 


(4.26) 


V  “V  i/ilt 

where  p'(y,t)  ■  p*(y,w)e  is  the  solution  of  the  wave  equation  (Eq.  4.4) 
for  the  following  boundary  conditions 


p(x,H,t)  ■  0 


& 


(x,0 , t) 


w  ijjt 

e 


and  Pj  (y,t) 


Pj  (y,w)e 


imt 


is  the  solution  of  the  wave  equation  (Eq.  4.4)  at 
0  for  the  following  boundary  conditions: 

p(x,H,t)  *  0 
3s 


(x,0,t) 


w  ,  ,  .  iut 

-  g 


(4.28) 


Solution  of  Eq.  4.4  for  the  boundary  conditions  of  Eq.  4.27  is 


xy,  ,  w  C 

V*'"5  =  gw 


(1+ct)  sin  -(H(7-y-' 


..  OiH  ^  aiH 

(1+a)  cos  —  +  x(l-a)  sxn  — 


(4.29) 


In  this  equation,  the  reflection  coefficient 

(Crwr/Cw)  -1 


(Cr wr/Cw)  +1 


(4.30) 


where  w^  and  Cr  are  the  unit  weight  and  P-wave  velocity  for  the  rock  material 


at  the  bottom  of  the  reservoir.  For  rigid  rock,  a  =  1  indicating  that  hydro¬ 
dynamic  pressure  waves  are  completely  reflected  at  the  bottom  of  the  reservoir. 
For  deformable  rock, a  <  1  implying  that  hydrodynamic  pressure  waves  imping¬ 
ing  at  the  reservoir  bottom  ar.e  partially  reflected  back  in  the  fluid  domain 
and  partially  refracted  into  rock.  The  complex  frequency  response  function 


Pq  (y,w)  is  for  the  hydrodynamic  pressures  on  the  dam  when  the  excitation  to 


the  fluid  domain  is  the  vertical  acceleration  of  the  ground  and  the  dam  is 
rigid. 


Hie  solution  of  Eq.  4.4  for  the  boundary  conditions  of  Eq.  4.28  has  been 


presented  earlier  in  Eq.  4.13.  As  mentioned  before,  p^ (y,u)  is  the  complex 


frequency  response  function  for  the  hydrodynamic  pressure  when  the  excitation 
is  the  acceleration  of  the  dam  in  the  mode  of  vibration  and  there  is  no 


motion  of  the  base. 


iut 


The  nodal  force  vectors  are  Rq (t)  ■  R^(u)e1U)t  and  (t)  *  Rj  (w)e 
5^  (to)  and  Rj  («)  are  the  static  equivalents  of  the  corresponding  pressure 


where 


functions  pg(y,w)  and  p^ (y,u>)  respectively.  The  complex  frequency  response 
function  for  the  total  force  vector  R?*(t)  is,  from  Eq.  4.25, 


-h 


Rjf(u)  -  5j(w)  +  t  H(u>)  R,  (u>) 


J  - 


(4.31) 


j-1 


N 
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With  the  aid  of  Eq.  4.31  and  substituting  (oj)  »  HiJ25y (u>)  ,  Eq.  4.22  can  be 

written  as  follows  for  the  excitation  ay(t)  -  elU)t: 

9 

j^-u)2Mn  +  itoCn  +  kJy^w)  -  -Ly  +  {$nf}T{S(W)  (4.32) 

The  set  of  Eqs.  4.32  for  n  *  1,2... J  after  rearranging  may  be  expressed  in 
matrix  form  as 


S(u>)Yy(w)  =  Ly(tD)  (4.33) 

in  which  the  coefficient  matrix  s  (u>)  is  identical  to  the  case  when  the  ex¬ 
citation  was  the  horizontal  component  of  ground  motion  (Eqs.  4.16-4.18),  and 

Ly(a>)  -  -£y( u)  +  {4Xf)T  Ry(U»,  n  -  1,2,...J  (4.34) 

n  n  *n  -o 

Solution  of  Eq.  4.33  for  all  values  of  excitation  frequency  would  give 
the  complex  frequency  response  functions  for  generalized  displacements 

•y 

Y^  (u>) .  The  frequency  responses  for  generalized  accelerations 

Yy(w)  «  -  u2Yy(C0)  (4.35) 

n  n 

The  complex  frequency  response  functions  for  the  vectors  of  nodal  displace¬ 
ments  and  accelerations  relative  to  the  ground  motion  are  given  by  expres¬ 
sions  similar  to  Eqs.  4.20  and  4.21: 

J 

ry (w)  =  Z  A  Yy(u»  (4.36) 

-  ,  -n  n 
n*l 

.  J 

P  (U>)  ■  -UZ  Z  4  Yy(u»  (4.37) 

-  ,  *n  n 
n*l 

4.5  Response  to  Arbitrary  Ground  Motion 

Once  the  complex  frequency  response  functions  YX(w)  and  Yy(w),  n  -  1, 

n  n 

2, . . . J,  have  been  determined  by  solving  Eqs.  4.16  and  4.33  for  an  appropri¬ 
ate  range  of  excitation  frequency  u,  the  responses  to  arbitrary  ground 
motion  can  be  obtained  as  the  superposition  of  responses  to  individual  har¬ 
monic  components  of  the  excitation  through  the  Fourier  integral: 
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vf(t)  -  /  ?f(0))  A^faOe^dUi  £  =  X,y  (4.38) 

n  in  —  °» *  n  g 

l  f 

in  which  A  (<u)  is  the  Fourier  Transform  of  a'Mt) : 
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A^(OJ)  «oJ  a^(t)e  iWtdt;  £  =  x,y  (4.39) 

where  d  is  the  duration  of  ground  motion. 

The  combined  response  Y^ (t)  to  horizontal  and  vertical  components  of 
ground  motion  acting  simultaneously  is 

Y  (t)  =  YX(t)  +  Yy(t)  (4.40) 

n  n  n 

Repeating  this  procedure  for  all  the  necessary  values  of  n,  the  displace¬ 
ment  response  is  obtained  by  transforming  back  from  generalized  to  nodal 
point  coordinates 


J 

r (t)  *  E  4  Y  (t)  (4.41) 

—  -  n  n 

n=l 

The  stresses  a^(t)  in  finite  element  p  at  any  instant  of  time  are  re¬ 
lated  to  the  nodal  displacements  r^(t)  for  that  element  by 

a  (t)  =  T  r  (t)  (4.42) 

S'  ¥  XT 

where  T  is  the  stress  transformation  matrix  for  the  finite  element  p.  At 

-p 

any  instant  of  time,  the  stresses  throughout  the  dam  are  determined  from 
the  nodal  point  displacements  by  applying  the  transformation  of  Eq.  4.42  to 
each  finite  element. 

The  total  stresses  are  the  sum  of  the  dynamic  stresses  determined  above 
and  the  intial  static  stresses  due  to  the  weight  of  the  dam  and  hydrostatic 
pressures  (Sec.  3.3). 

A  digital  computer  program  was  developed  [20]  to  implement  the  analysis 
summarized  in  this  chapter.  Studies  using  this  computer  program  have  led  to 
a  better  understanding  of  hydrodynamic  effects  in  earthquake  response  of 
concrete  gravity  dams  [18] . 
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5.  ANALYSIS  OF  DAMS  INCLUDING  DAM-FOUNDATION  INTERACTION 


5.1  Introductory  Note 

This  chapter  summarizes  the  procedure  presented  earlier  [24,25]  for  in¬ 
cluding  the  effects  of  structure-foundation  interaction,  i.e.,  interaction 
between  the  structure  and  underlying  soil  or  rock,  in  earthquake  response 
analyses  of  concrete  gravity  dams  without  water.  This  procedure  to  in¬ 
clude  the  dam- foundation  interaction  effects  is  based  on  the  same  basic 
substructure  concept  utilized  to  include  the  dam-water  interaction.  The 
soil  or  rock  region  and  the  dam  are  considered  as  two  substructures  of  the 
combined  system.  The  dam  may  be  idealized  by  the  finite  element  method, 
which  has  the  ability  to  represent  systems  of  arbitrary  geometry  and 
material  properties  (see  Chapters  3  and  4) .  The  soil  or  rock  region  may  be 
idealized  as  either  a  finite  element  system  or  as  a  continuum,  for  example, 
as  a  viscoelastic  halfspace.  The  halfspace  idealization  permits  accurate 
modelling  of  sites  where  essentially  similar  rocks  extend  to  large  depths. 
For  sites  where  layers  of  soil  or  soft  rock  are  underlain  by  harder  rock  at 
shallow  depth,  finite  element  idealization  of  the  soil  region  would  be  ap¬ 
propriate.  The  governing  equations  for  the  two  substructures  are  combined 
by  imposing  equilibrium  and  compatibility  requirements  at  the  base  of  the 
dam,  which  is  the  structure-foundation  interface. 

5. 2  Frequency  Domain  Equations  in  Nodal  Point  Coordinates 
5.2.1  Substructure  1;  Dam 

The  equations  of  motion  for  the  dam  idealized  as  a  planar,  two-dimen¬ 
sional  finite  element  system  (Chapter  3)  are: 


m  r  +  cr  +kr 
-c-c  -c-c  -c-c 


-m  lx  ax(t)  -  m  ly  ay(t)  +  R  (t) 
-c-c  g  -c-c  9  -c 


(5.1) 


in  which  m  ,  c  ,  and  k  are  the  mass,  damping,  and  stiffness  matrices  for 
— c  — c  -c 

the  finite  element  system;  rc  is  the  vector  of  nodal  point  displacements 
relative  to  the  free-field  earthquake  displacement  at  the  base  (Fig.  5.1) : 


rT  -  <rx  ry  rx  r7 
-c  rl  1  z2  *2 


x  y 
r  r 3 
n  n 


r  r*  > 
N+N  N+N. 

D  b 


in  which  rx  and  ry  are  the  x  and  y  components  of  the  displacement  of  nodal 
n  n 


****** 


point  n;  N  is  the  number  of  nodal  points  above  the  base,  N,  is  the  number  on 

o 

the  base,  and  the  total  number  is  N+N,  ; 

D 


{lX}T  =  <1  0  1  0  ....1  0  .  .  .  .  1  0> 

-c 

T 

{ly}  =  <0  1  0  1  ....0  1  .  .  .  .  0  1> 

-c 


R  (t) ,  the  vector  of  forces  due  to  structure-foundation  interaction,  will 
-c 


X  V 

have  non  zero  elements  only  at  the  base  of  the  dam;  a  (t)  and  a  (t)  are  the 
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x  and  y  components  of  the  specified  free-field  ground  acceleration,  assumed 


to  be  identical  at  all  nodal  points  on  the  dam  base. 


For  harmonic  ground  acceleration  in  the  x  (horizontal)  or  y  (vertical) 


directions  a*"(t)  *  e1U)t,  t  -  x  or  y,  the  displacements  and  interaction  forces 
9 

can  be  expressed  in  terms  of  their  complex  frequency  responses:  r  (t)  = 

-t ,  .  iwt  ,  _  -t .  ,  icot  ... 

r  (o>)  e  and  R  (t)  =  R  (dl)  e  and  the  governing  equation  as 
-c  — c  -c 


_  p 

P-di^m  +  iuc  +  k  lr  (w)  *  -m  1*"  +  RC  (<o) ;  t  *  x  or  y 

*-  Jr  — C  —C  —  r*  * 


-l, 


c-c 


(5.2) 


If  the  structure  has  constant  hysteretic  damping,  Eq.  5.2  becomes 


[-(O^m  +  (1  +  ir|)k  Ir^(a')  =  -m  l't  +  5°  (u)) ;  £  =  x  or  y 

—  c  — cJ  —  c  — c~c  —  c 


t 


(5.3) 


Partitioning  rc  into  r,  the  vector  of  displacements  of  nodal  points 

above  the  base,  and  r.  the  vector  of  interaction  displacements  at  the  base 

“D 

(Fig.  5.1),  Eq.  5.3  can  be  expressed  as 
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5.2.2  Substructure  2:  Foundation  Rock  Region 

The  dynamics  of  the  foundation  soil  or  rock  region  are  considered  sep¬ 
arately.  The  interaction  forces,  R^(t),  at  the  foundation  surface  can  be 
expressed  in  terms  of  the  corresponding  displacements,  r^  (t) ,  relative  to 
the  free-field  earthquake  ground  displacement  [24] : 


PIG. 


SUBSTRUCTURE  2:  FOUNDATION  ROCK  REGION 


5.1  SUBSTRUCTURE  REPRESENTATION  OP  DAM- FOUNDATION  SYSTEM 
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(5.5) 


3f(w)  rf(u))  =  Rf(w) 

where  3-(u»)  is  the  complex  valued,  dynamic  stiffness  matrix  for  the  founda- 
tion  soil  region.  The  ij  element  of  this  matrix,  S^to)^  is  defined  in 
Fig.  5.2  where  displacements  as  shown  have  been  imposed  at  nodal  points  with¬ 
in  the  base  of  the  dam  and  tractions  outside  the  base  are  zero. 

The  foundation  stiffness  matrix  (oj)  is  to  be  determined  by  a  separate 
analysis  of  the  substructure  representing  the  foundation  soil  or  rock  region. 
For  sites  where  formations  of  soft  rock  or  soil  overlie  hard  rock,  it  would 
be  reasonable  to  hypothesize  that  the  boundary  between  the  two  types  of  rock 
is  rigid.  Then  the  foundation  region  may  be  idealized  as  a  finite  element 
system,  in  which  material  non-homogeneity  and  irregular  geometry  can  be  con¬ 
veniently  handled.  Techniques  for  determining  for  a  finite  element 

idealization  of  the  foundation  soil  region  are  discussed  elsewhere  [29]. 

For  sites  where  similar  rock  extends  to  "arge  depths,  it  may  be  more  appro¬ 
priate  to  idealize  the  foundation  as  a  viscoelastic  half-plane.  Assuming 
that  the  dam  is  supported  at  the  surface  of  a  viscoelastic  half-plane  with 
homogeneous  material  properties,  (to)  can  be  determined  from  the  available 
data  and  procedures  [32].  These  results,  in  contrast  to  those  from  many 
earlier  studies,  are  not  restricted  to  a  rigid  plate  at  the  base  of  the 
structure. 

5.2.3  Dam- Foundation  System 

Equilibrium  of  interaction  forces  between  the  two  subs  true  tur  e  s  —  dam 
and  foundation  —  at  the  base  of  the  dam  requires  that 


-Rf  (u) 


(5.6) 


Similarly,  compatibility  of  interaction  displacements  in  the  two  substruc¬ 
tures  at  the  dam  base  requires  that 


r£(u>)  *  rf  (w) 


(5.7) 


Consequently,  Eq.  5.5  becomes 


-£f(w)  r£(u>) 


(5.8) 
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FIG.  5.2 


DEFINITION  OF  THE  ij -ELEMENT  OF  THE  DYNAMIC 

STIFFNESS  MATRIX  £,(u>)  FOR  THE  FOUNDATION  REGION 


which  is  substituted  in  Eg,  5.4  to  obtain 


[-“2Sc 


+  (1  + 


in)k 

— C 


2f  (oi) 


]2 


-l  t 
;(w)  =  -m  lc 
c  -c-c 


(5.9) 


where 


|  f(u» 


0 

if  (w> 


(5.10) 


For  a  particular  excitation  frequency  w,  Eq.  5.9  represents  2 (N+N,  )  alge- 

b 

braic  equations  with  complex-valued  coefficients,  in  the  unknown  nodal  point 
-l 

displacements  r^. 

5.2.4  Computational  Requirements 

Analysis  of  structural  response  to  earthquake  motion  via  the  frequency 
domain  requires  that  the  complex  frequency  response  functions  be  determined 
for  a  range  of  frequencies  over  which  the  ground  motion  and  structural  re¬ 
sponse  have  significant  components.  For  each  excitation  frequency,  this 
entails  solution  of  Eq.  5.9,  a  set  of  as  many  algebraic  equations  as  the 
number  of  degrees  of  freedom  for  the  dam.  Finite  element  systems  with  a  few 
hundred  degrees  of  freedom  are  typically  employed  to  represent  concrete 
gravity  dams.  Enormous  computational  effort  would  be  required  for  repeated 

solutions  of  these  equations  for  many  values  of  the  excitation  frequecy. 

I 

Furthermore,  the  entire  displacement  vector,  r  ,  at  each  u>  has  to  be  stored 

— c 

for  future  Fourier  synthesis  ,  requiring  very  large  amounts  of  computer 
storage.  Fourier  synthesis  calculations  would  have  to  be  performed  for  each 
of  the  few  hundred  displacements,  requiring  large  computational  effort  even 
when  the  very  efficient  FFT  procedure  is  employed.  The  need  for  judiciously 
reducing  the  number  of  degrees  of  freedom  in  the  analysis  is  therefore 
apparent . 

5.3  Reduction  of  Degrees  of  Freedom 

At  any  instant  of  time,  the  structural  displacements  in  the  DOF  above 
the  base  may  be  separated  into  static  displacements  due  to  the  interaction 
displacements  at  the  base  and  the  remaining  dynamic  displacements.  The 


f 
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latter  can  be  expressed  as  a  linear  combination  of  the  first  few  natural 
modes  of  vibration  of  the  structure  on  a  fixed  base.  This  has  been  demon¬ 
strated  to  be  effective  in  reducing  the  number  of  DOF  for  one-dimensional 
structural  systems  on  rigid  footing  supported  on  the  surface  of  a  halfspace 
and  also  for  complex  structures  on  a  deformable  base.  If  J  vibration  modes 
are  adequate  to  represent  the  structural  response,  the  number  of  unknowns 
in  Eq.  5.9  would  be  reduced  to  J+2N^#  the  latter  being  the  number  of  DOF  at 
the  base.  While  this  approach  is  conceptually  appealing,  the  number  of  un¬ 
knowns  can  not  be  reduced  below  2N^,  the  number  of  DOF  at  the  base,  which  in 
finite  element  models  of  dams  might  be  of  the  order  of  20. 


For  such  systems  the  number  of  unknowns  in  Eq.  5.9  can  be  most  effec¬ 
tively  reduced  by  use  of  the  Ritz  concept.  The  displacements  rc  are  ex¬ 
pressed  as  linear  combinations  of  the  Ritz  vectors,  chosen  as  the  normal 
modes  of  an  assoicated  undamped  structure-foundation  system.  The  associated 
system  considered  is  one  in  which  ^(u)  is  replaced  by  a  frequency- 
independent  value,  say  the  static  value  3^(0).  The  frequencies  X^  and  mode 
shapes  of  the  associated  system  are  solutions  of  the  eigenvalue  problem 

|k  +  } -(0)1  $  -  X2  m  to 
I -c  -f  rn  n  -c  rn 


(5.11) 


.  to,  are  considered  as  the  Ritz  vectors, 

”  J 


(5.12) 


where  £f(0)  was  defined” in  Eq.  5.10. 

If  the  first  J  modes  .. 

the  structural  displacements  are  expressed  as 

rf  (t)  =  f  zf(t)  to  l  -  x,y 
j=l  3  "3 

where  Z/s  are  the  generalized  coordinates.  In  terms  of  the  complex 
frequency  response  functions,  Eq.  5.12  is 

*  l  fc. 

C  j-1  3  3 


(5.13) 


Introducing  the  transformations  of  Eq.  5.13  into  Eq.  5.9,  premultiply- 
T 

ing  by  and  utilizing  the  orthogonality  property  of  modes  of  the  associ¬ 
ated  structure-foundation  system  with  respect  to  the  "stiffness"  and  mass 
matrices  of  Eq.  5.11,  results  in 

S(w)  !*(«)  -  t*(w)  (5.14) 
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In  Eq.  5.14,  the  diagonal  elements  of  S (w)  are 


S„  (w) 
nn 


-CO2  +  (1  +  in)X2UT  m  f  +  vTf|.((0)  -(l+in)iL(0)]  ib  (5.15a) 
!_  nj -n  -c  ■‘•n  nl-f  -t  -n 


whereas  the  off-diagonal  elements  are 


S  .  (co) 


IX 

*n 


(ui)  -  (l  +  in)  3f (0) 


•]*j 


(5.15b) 


where  n  and  j  =  1,2,3, Z(oo)  is  a  column  vector  of  the  complex  frequency 
response  functions  for  the  generalized  coordinates  Z^;  I  is  a  column  vector 
with  the  nfck  element 


n 


(5.15c) 


For  a  particular  excitation  frequency  to,  Eq.  5.14  represents  J  simultaneous 
algebraic  equations  in  the  generalized  coordinates  Z^Oo),  n  -  1,2,...,J. 

If  all  the  2(N+Nfa)  vibration  modes  of  the  associated  structure- 
foundation  system  were  included  in  the  Ritz  method,  the  results  obtained  by 
solving  the  simultaneous  equations  5.14  would  be  identical  to  the  "exact” 
results  from  solution  of  Eq.  5.9.  There  would,  however,  be  no  advantage  in 
choosing  to  solve  Eq.  5.14  instead  of  Eq.  5.9  because  the  number  of  equa¬ 
tions  is  not  reduced  nor  are  they  uncoupled.  On  the  other  hand,  if  suffi¬ 
ciently  accurate  results  can  be  obtained  with  J«2(N+Nb),  there  would  be 
profound  computational  advantages  in  working  with  Eq.  5.14. 

That  this  is  indeed  the  case  has  been  demonstrated  for  one-dimensional 
structural  systems  on  rigid  footings  supported  on  the  surface  of  a  halfspace 
[33]  and  also  for  concrete  gravity  dams  [24].  The  number  of  generalized 
coordinates  that  should  be  included  in  analysis  depends  on  the  properties  of 
the  structure-foundation  system,  the  response  quantities,  and  the  frequency 
range  of  interest.  The  number  necessary  in  analysis  of  concrete  gravity 
dams  studied  in  this  report  is  determined  in  Section  7.6. 

5.4  Response  to  Arbitrary  Ground  Motion 

-4C  — y 

The  complex  frequency  response  functions  Z  (u>)  and  ZJ  (w) ,  n  »  1, 

n  n 

2,...,  J,  are  determined  by  solving  Eq.  5.14  for  an  appropriate  range  of 
excitation  frequency  w.  The  response  to  arbitrary  ground  motion  can  be 
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computed  using  the  equations  of  Sec.  4.5  with  the  following  changes  in 

l  / 

notation:  replace  Y  by  Z  ,  r  by  r  , and  d>  by  W. 

n  n  -  -c'  in  *-n 
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6.  ANALYSIS  OF  DAMS  INCLUDING  HYDRODYNAMIC  AND  FOUNDATION 
INTERACTION  EFFECTS 


6.1  Introductory  Note 

The  substructure  concept  was  employed  to  separately  include  the  effects 
of  dam-water  interaction  (Chapter  4)  and  dam-foundation  interaction  (Chapter 
5)  in  the  analysis.  Using  the  same  concept  the  analysis  procedure  is  ex¬ 
tended  to  analyze  earthquake  response  of  dams,  simultaneously  including 
various  effects  of  the  water  and  the  foundation.  These  include  effects 
arising  from  interaction  between  the  dam  and  foundation,  dam  and  water, 
water  and  foundation,  and  from  interaction  among  all  three  substructures  — 
dam,  water  and  foundation. 


6. 2  Frequency  Domain  Equations 
6.2.1  Substructure  1 :  Dam 


A  cross-section  of  the  dam-water-foundation  system  is  shown  in  Fig.  2.1. 
The  equations  of  motion  for  the  dam  idealized  as  a  planar,  two-dimensional 
finite  element  system  ares 


mr  +  cr  +kr 
-c-c  -c-c  -c-c 


-m  lxax (t)  -  m  lyay (t) 

-C-C  CJ  — c— c  ^ 


+  R  (t) 

“C 


(6.1) 


in  which  m  ,  c  ,  and  k  are  the  mass,  danping  and  stiffness  matrices  for 
— c  -c  — c 

the  finite  element  system;  rc  is  the  vector  of  nodal  point  displacements 
relative  to  the  free-field  earthquake  displacement  at  the  base  (Fig.  6.1) : 


T 


r 

-c 


x  y 
r  r 
r2  2 


x  y 
r  r 
n  n 


in  which  rX  and  ry  are  the  x  and  y  components  of  the  displacement  of  nodal 
point  n;  N  is  the  number  of  nodal  points  above  the  base,  N^  is  the  number  of 
nodal  points  or  the  base,  and  the  total  number  is  N+Nfa; 


{lX}T  -  <1010 . 10 . 10> 

“C 

{ly}T  -  <0101 . 01 . 01> 

*"C 


The  force  vector  Rc (t)  includes  hydrodynamic  forces  R^ (t)  at  the  upstream 
face  of  the  dam  and  forces  R^ (t)  at  the  base  of  the  dam  due  to  structure- 


x  y 

foundation  interaction;  a  (t)  and  a  (t)  are  the  x  and  y  components  of  the 
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specified  free-field  ground  acceleration,  assumed  to  be  identical  at  all 


nodal  points  on  the  dam  base. 


For  harmonic  ground  acceleration  in  the  x  (horizontal)  or  y  (vertical) 

be 
iojt 


directions  a^(t)  *  e*Ut,  t  *  x  or  y,  the  displacements  and  forces  can  be 

g  -Li 

expressed  in  terms  of  their  complex  frequency  responses:  r^(t)  *  £c  (u>)  e 

Rc(t)  =  R^(t)  =  R^  (to)  and  Rj^(u))  =  •  The  governing 


equation,  Eq.  6.1,  then  becomes 

1-e 


2  1-i 

-oj  m  +  itoc  +  k  r 
-c  -c  -cj  -c 


(to) 


/  -l 

-m  l*'  +  R*-(u) 
-c-c  “C 


x  or  y 


(6.2) 


If  the  structure  has  constant  hysteretic  damping,  Eq.  6.2  becomes 


-c 


-U)*m^  +  (l+in)kcJr^  (a) )  =  -m_l^  +  R^(oj)  ;  L  =  x  or  y 


■U 


l  .  at. 


-c-c  -c 


(6.3) 


Partitioning  r  into  r,  the  vector  of  displacements  of  nodal  points 
~c  — 


above  the  base,  and  r,  ,  the  vector  of  interaction  displacements  at  the  base 

“D 

(Fig.  6.1),  Eq.  6.3  can  be  expressed  as 


m  0 

2  2b 


+  d+in) 


*b 


*bb 


-l  1 

£  (W) 

if 

5(“l 

-l 

'  s  —  ^ 

,1 

•  +  • 

1  -i,  „ 

•  (6.4) 

rb(u)) 

sA. 

[vu)j 

-l, 


The  forces  at  the  base  of  the  dam,  (w)»  are  next  expressed  in  terms  of 
interaction  displacements  by  appropriate  analyses  of  substructure  2,  the 
foundation  (Sec.  6.2.2).  Later  (Sec.  6.2.5)  the  hydrodynamic  forces  R^w) 
are  expressed  in  terms  of  displacements  or  accelerations  at  the  upstream 
face  of  the  dam  by  appropriate  analyses  of  substructure  3,  the  fluid 
region. 


6.2.2  Substructure  2:  Foundation  Rock  Region 

The  forces  acting  at  the  surface  of  the  foundation  include  the  forces 
Rf  at  the  base  of  the  dam  due  to  structure-foundation  interaction  and  the 
hydrodynamic  forces  at  the  bottom  of  the  reservoir.  Beyond  a  certain 
distance  upstream  of  the  dam,  the  hydrodynamic  pressures  will  become  small 
enough  to  be  negligible.  Only  the  significant  forces  are  included  in  J>^. 
For  unit  harmonic  ground  acceleration  these  forces  can  be  expressed  in 
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terms  of  their  complex  frequency  response  functions 
Q^(t)  ■  Q^(u)ela)t.  The  corresponding  displacements 
field  ground  dispalcement  cure  rf(t)  =  r^(U))eia)t  and 
forces  and  displacements  can  be  related  through  the 
foundation  stiffness  matrix  $ (u) : 


Rf(t)  -  Rf  (oj)e1Wt  and 
relative  to  the  free- 
g(t)  =  g(o;)ela,t.  The 
complex  valued,  dynamic 


"3  „(w) 

3_  (a))' 

•mm 

r. 

r— 

r. 

-rr 

-rq 

4 

-f 

■  = . 

-f 

£ ,«“> 

3  (u) 

-qq  J 

A 

(6.5) 


The  ij^  element  of  this  matrix,  ^  (w) ,  is  defined  in  Pig.  6.2  where  dis¬ 
placements  as  shown  have  been  imposed  at  nodal  points  contained  in  r^  and  g 
and  tractions  outside  these  nodal  points  are  zero. 


The  second  matrix  equation  from  Eq.  6.5  can  be  expressed  as 


3 


3-1  (a>) 

-qq 


(w)  rf 


(6.6) 


Substituting  this  expression  into  the  first  matrix  equation  from  Eq.  6.5  leads 
to 

I*  (CO)  -  8  (<o)  3*  <u>f|  -  R.  -  $  (u)  3_1(u»  Q. 

|_-rr  -rq  -qq  -rq  J  -f  -f  -rq  -qq  -h 

or 

'  3_(u)r-  =  R,  -  3  (w)  3-1  (w)  Q.  (6.7a) 

-£  -t  -f  -rq  -qq  -h 

where 

3-<w)  =  3  (w)  -  3  (w)  3“i (u>)  3*  (w)  (6.7b) 

— r  -rr  -rq  -qq  -rq 

The  matrix  3f  (w)  of  Eq.  6.7b  is  the  same  as  the  matrix  3f(a))  in  ^3*  5.5.  In 
particular,  without  hydrodynamic  forces,  gh  ®  0  and  Eq.  6.6  reduces  toEq. 5.5. 

The  dynamic  stiffness  matrices  for  the  foundation  appearing  in  Eqs.  6.5, 
6.7a,  and  6.7b  are  to  be  determined  by  a  separate  analysis  of  the  substruc¬ 
ture  representing  the  foundation  soil  or  rock  region.  At  many  sites  of  con¬ 
crete  gravity  dams ,  similar  rocks  extend  to  large  depths ,  and  it  is  therefore 
appropriate  to  idealize  the  foundation  as  a  viscoelastic  half  plane.  The 
surface  of  the  half  plane  is  assumed  to  be  horizontal  and  its  material  prop¬ 
erties  to  be  homogeneous.  Available  data  and  methods  [32)  are  utilized  in 
this  investigation  to  determine  dynamic  stiffness  for  such  idealized  foundations. 


FIG.  6.2  DEFINITION  OF  Sf(u>)  ,  THE  ij-ELEMENT  OF  THE  DYNAMIC 
STIFFNESS  MATRIX  2f(u>)  FOR  THE  FOUNDATION  REGION 
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6.2.3  Dam- Foundation  System 

Equilibrium  of  interaction  forces  between  the  two  substructures  —  dam 
and  foundation  —  at  the  base  of  the  dam  requires  that 

(6.8) 


|£(u» 


-Rf  (W) 


Similarly,  compatibility  of  interaction  displacements  at  the  dam  base  in  the 
two  substructures  requries  that 


r£(w)  =  rf (u) 


Consequently,  Eq.  6.6  because 

-l 


(w)  -  -  S,(ui)  r^(oj)  -  i  Q 

-f  — b  -rq  -qq  -n 


(6.9) 


(6.10) 


—  n 

Thus,  the  forces  R£(<d)  at  the  base  of  the  dam  due  to  structure-foundation 

interaction  have  been  expressed  in  terms  of  the  interaction  displacements 

rj-(ui)  and  hydrodynamic  forces  Q.  through  dynamic  stiffness  matrices  for  the 
-b  -h 

foundation  region.  Substituting  Eq.  6.10  into  Eq.  6.4  leads  to 


m 


0  2b 


+  (l+in) 


*b‘ 

0  0 

\ 

r£(u)' 

«T» 

+ 

■ 

-l 

S  *bbj 

o  |f(u>) 

/ 

or 


where 


m  1 

.1 


+  S 


s£‘“> 

’-rq-qq-hj 


f-w2  m  +  (l+in)k  +  3-(w)]r^(ca)  =  -m  +  R  (u) 
— c  -c  -f  J  -c  — c— c  — e 


|f(id) 


0  0 


0  2f(w) 


(6.11) 


(6.12) 


(6.13) 


The  hydrodynamic  forces  R^  and  will  be  expressed  in  terms  of  acceler¬ 
ations  of  the  upstream  face  of  the  dam  and  the  bottom  of  the  reservoir  by 
analyses  of  the  fluid  region,  to  be  described  later. 


6.2.4  Reduction  of  Degrees  of  Freedom 

For  each  excitation  frequency  Eq.  6.11  is  a  set  of  2(N+Nb)  algebraic 
equations,  as  many  as  the  number  of  DOF  of  the  dam,  which  typically  would 
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be  a  few  hundred.  Enormous  computational  effort  would  be  required  for  re¬ 
peated  solutions  of  these  equations  for  many  values  of  the  excitation  fre¬ 
quency  and  subsequent  Fourier  synthesis  of  the  harmonic  responses  in  each 
DOF.  Thus,  it  is  important  to  judiciously  reduce  the  number  of  unknowns  in 
the  analysis.  As  mentioned  in  Chapter  5,  two  approaches  have  proven  to  be 
effective  in  reducing  the  number  of  DOF. 

I 

In  the  first  approach  [25,  33],  the  structural  displacements  r  (t)  at 

some  time  t  in  the  DOF  above  the  base  are  separated  into  two  parts:  static 

l 

displacements  due  to  base  displacements  r,  (t)  at  the  same  time  plus  the  re- 

—  D 

maining  "dynamic"  displacements;  and  the  latter  are  expressed  as  a  linear 
combination  of  the  first  few  natural  modes  of  vibration  of  the  structure  on 
a  fixed  base.  If  J  vibration  modes  are  adequate  to  represent  the  structural 
response,  the  number  of  unknowns  in  Eq.  6.11  would  reduce  to  J+2N,  ,  where 

D 

2Nfa  =  the  number  of  DOF  at  the  base.  While  this  approach  is  conceptually 

appealing,  the  number  of  uknowns  cannot  be  reduced  below  2N  ,  which  in  finite 

b 

element  models  of  dams  may  be,  say,  20. 

The  second  approach  [24,33],  based  on  the  Ritz  concept,  is  more  effective 
for  such  systems.  The  displacements  rc  are  expressed  as  linear  combinations 
of  the  Ritz  vectors,  chosen  as  the  normal  modes  of  an  associated  undamped 
structure- foundation  system.  The  associated  system  chosen  here  and  in  earlier 
studies  is  one  in  which  (w)  is  replaced  by  a  frequency  independent  value, 
say  the  static  value  ^(0).  The  vibration  frequencies  X^  and  mode  shapes 
of  the  associated  system  are  solutions  of  the  eigenvalue  problem 


0= +  if(o>]!„ 

where  (0)  was  defined  in  Eq.  6.13. 


X^m  to 
n-cxn 


(6.14) 


If  the  first  J  modes  .  .  .  $ ^ 

the  structural  displacements  are  expressed  as 

J 


are  considered  as  the  Ritz  vectors. 


r*(t) 

“C 


j 


(6.15) 


where  Zj's  are  the  generalized  coordinates.  In  terns  of  the  complex  frequency 
response  functions,  Eq.  6.15  is 


r*>)  -  l  zf  (0J)^. 
c  j-1  3  3 


(6.16) 
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Equation  6.16  contains  the  following  equation  for  the  displacements  at  the 
base  of  the  structure: 


rb(to)  =  l  Z  <U3) ijpb  (6.17) 

j“l  3  3 

where  4  .  is  the  subvector  of  ip.  corresponding  to  the  base  DOF.  Because  r  =r  , 
~°3  “3  -f  -b 

and  4>.  were  determined  from  Eq.  6.14  involving  (0),  the  static  foundation 

stiffness  matrix,  g  may  be  expressed  in  terms  of  generalized  coordinates  by 

setting  to  =  0  in  Eq.  6.6  and  ignoring  the  hydrodynamic  loads  Q;  thus 

q  =  -3~*(0)  S*(0)r.  (6.18) 

-qq  -rq  -b 

Combining  Eqs.  6.17  and  6.18  leads  to 

g=  l  Z  (to))L.  (6.19) 

j=l  3  “ 

where 

5,‘o>*bj  <6-2°> 

Introducing  the  transformation  of  Eq.  6.16  into  Eq. 6. 12, premultiplying 
T 

by  j|j  and  utilizing  the  orthogonality  property  of  eigenvectors  of  the  associ¬ 
ated  structure-foundation  system  with  respect  to  the  stiffness  and  mass 
matrices  of  Eq.  6.14,  results  in 

-l  l 

S(00)  Z  (to)  =  L  (to)  (6.21) 

In  Eq.  6.21,  the  diagonal  elements  of  §  (to)  are 


nn 


(CO)  =  -oo2  +  (1+iri)*2  £Tm  \j)n  +  ii)T L2 _ (co)  -  (l+in)iL(0)  (6.22a) 

nj  n-c1  n  ■Ln|j-t  -t  J  -n 


whereas  the  off-diagonal  elements  are 

,f(co)  -  (l+in)|f  (0) 


S  .  (oo) 
nj 


- 


* 


(6.22b) 


rl. 


where  n  and  j  *  1,2,3 - J;Z  (oo)  is  a  column  vector  of  the  complex  frequency 

l  i 

response  functions  for  the  generalized  coordinates  Z  ;  L  is  a  column  vector 

th  ^  — 

with  the  n  element  »  » 

L  *  -LC  +  4#  R  (oo)  (6.23) 

n  n  -n-c 


where 


.T  .1 
L  ■  ii/  m  1 
n  zn-c-c 


(6.24) 
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J*. 


...  u* 


For  a  particular  excitation  frequency  oj,  Eq.  6.21  represents  J  simultaneous 

/ 

algebraic  equations  in  the  generalized  coordinates  Z  (w),n*l,2---J. 

n 

These  equations  need  to  be  solved  for  several  hundred  values  of  the  excita¬ 
tion  frequency  to  determine  the  complex  frequency  responses . 

Equation  6.21  would  contain  2(N+N.  )  equations  if  all  the  vibration 

D 

modes  of  the  associated  structure-foundation  system  where  included  as  Ritz 
vectors.  Solution  of  these  equations  would  then  be  identical  to  the  solu¬ 
tion  of  Eq.  6.11  in  nodal  point  coordinates.  There  would,  however,  be  no 
advantage  in  choosing  to  solve  Eq.  6.21  instead  of  Eq.  6.11  because  the 
number  of  equations  is  not  reduced  nor  are  the  equations  uncoupled  by  the 
transformation.  However,  there  would  be  profound  computational  advantages 
in  working  with  Eq.  6.21  if  sufficiently  accurate  results  can  be  obtained 
by  including  only  a  few  Ritz  vectors.  The  number  necessary  depends  on  the 
properties  of  the  structure-foundation  system,  the  response  quantities,  and 
the  frequency  range  of  interest;  for  concrete  dams  this  number  is  determined 
in  Section  7.6. 


6.2.5  Substructure  3;  Fluid  Domain 

The  motion  of  the  dam-water  system  is  assumed  to  be  two-dimensional, 
and  the  same  for  any  vertical  plane  perpendicular  to  the  axis  of  the  dam 
(Chapter  2 ) .  Assuming  water  to  be  linearly  compressible  and  neglecting  its 
internal  viscosity,  the  small  amplitude  irrotational  motion  of  water  is 
governed  by  the  two-dimensional  wave  equation 


_i_ 

C2  3t2 


(6.25) 


in  which  p(x,y,t)  is  the  hydrodynamic  pressure  (in  excess  of  the  hydrostatic 
pressure)  and  the  velocity  of  sound  in  water  C  =  ►'gK/w  where  K  and  w 
are  bulk  modulus  and  unit  weight  of  water,  respectively,  and  g  =  accelera¬ 
tion  of  gravity.  For  harmonic  ground  motion  a^(t)  =  eiwt,  £  *  x  or  y, 


p(x,y,t) 


g(x,y,u))e 


io)t 


(6.26) 


where  p  is  the  complex  frequency  response  function  for  the  pressure  and  Eq. 
6.25  becomes  the  Helmholtz  equations 


3x2  3y2 


+ 


0 


(6.27) 
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Whereas  Eqs.  6.11  and  6.17  are  valid  for  dams  of  arbitrary  geometry,  it  is 
assumed  that  for  purposes  of  defining  hydrodynamic  loads,  and  only  for  this 
purpose,  the  upstream  face  of  the  dam  is  vertical.  This  is  a  reasonable 
assumption,  given  that  the  upstream  face  of  concrete  gravity  dams  is  verti¬ 
cal  or  almost  vertical  for  most  of  the  height. 

Dynamic  water  pressures  are  generated  by  horizontal  motions  of  the 
vertical  upstream  face  of  the  dam  and  by  vertical  motions  of  the  horizontal 
bottom  of  the  reservoir.  The  vector  of  horizontal  accelerations  at  the 
nodal  points  on  the  upstream  face  of  the  dam  is 


u(t>  -  lx  +  i  ZmL* 

l  j=l  3  3  > 


(6.28a) 


if  the  excitation  is  horizontal  ground  motion  a*(t)  =  elwt;  and 


u(t)  =  |  \  Z^(oj)j-€ 


(6.28b) 


if  the  excitation  is  vertical  ground  motion  ay(t)*e  .  In  Eq.  6.28  the 

xf  ® 

vector  is  a  subvector  of  ^  containing  only  the  elements  corresponding 

to  the  x-DOF  of  the  nodal  points  on  the  upstream  face  of  the  dam.  From 

Eq.  6.19,  the  vector  of  vertical  accelerations  at  the  nodal  points  at  the 

bottom  of  the  reservoir  is 


=  {  I  Xj  Z^(0))Je1( 


(6.29a) 


if  the  excitation  is  horizontal  ground  motion;  and 


v(t)  -  (l  +  l  £  ^(ai)}e1( 
(  j=l  3  3  > 


(6.29b) 


if  the  excitation  is  vertical  ground  motion,  where  is  a  subvector  of  ^ 
containing  only  the  elements  corresponding  to  the  y-DOF  of  the  nodal  points 


at  the  bottom  of  the  reservoir. 


XX 

Strictly  for  purposes  of  notational  convenience,  the  vectors  ^  and 
are  replaced  by  their  continuous  analogue  functions  ^ (y)  and  Xj (x) t  re¬ 
spectively  in  Eq.  6.28  to  obtain 


* 


u(0 ,y,t)  *  fl  +  l  i|>.  (y)  ZX (u)l< 
L  j-i  3  3  J 


iut 


(6.30a) 


u(0,y,t)  *  j”  I  .  (y)  Z^(u)le 

Lj=l  3  3 


iut 


and  in  Eq.  6.29  to  obtain 


v(x,0,t)  =  Xj  (x)  Z*(u)J< 


iut 


(6.:J0b) 


(6.31a) 


;(x,o  »t)  -  [l  +  1  Xj(x)  2^(0))], 


iut 


(6.31b) 


x  iut 

When  the  excitation  is  horizontal  ground  motion  a  (t)  =  e  ,  neglect¬ 
ing  effects  of  the  waves  at  the  free  surface  of  water  [11] ,  the  boundary 
conditions  for  Eq.  6.27  governing  p  are  (with  reference  to  Fig.  6.1) 


(6.32) 


For  the  case  of  vertical  ground  motion  a 
are 


y 

g 


(t) 


p(x,H,U)  -  0 


iut 

e 


the  boundary  conditions 


|i(x,0,u)  -  -  |  [l  +  J  Xj  (x)  Z ^  (u)J 

§§<o,y,<*»  -  -  |[  I  Vyl  *j(u)] 


(6.33) 


Because  the  governing  equ^ion  as  well  as  the  boundary  conditions  are 
linear,  the  principle  of  superposition  applies.  The  complex  frequency  re¬ 
sponse  functions  px(x,y,u>)  and  Py(x,y,u),  associated  respectively  with  x  and 
y  components  of  ground  motion,  can  therefore  be  expressed  as 
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(6.34) 


■t  t  ^  ~z.  r  f  _k  "] 

5  (x,y,w)  ■  "p^(x,y,u)  +  £  Zj(w)  [p^  (x,y,w)  +  P?(x,y,o))J 


In  Eq.  6.34,  p*(x,y,u))  is  the  solution  of  Eq.  6.27  for  the  following  bound¬ 
ary  conditions : 


p (x,H,ii>)  -  0 


ii 

3y 


(x,0,0i)  -  0 


(6.35) 


-  -  j 


— y 

pi(x,y,oi)  is  the  solution  of  Eq.  6.27  for  the  following  boundary  conditions: 


|5 (x,H,u>)  *  0 


3y 


(x,0,W) 


w 

9 


(6.36) 


|^(0, y,u)  =  0  ;  J 


pf (x,y,w)  is  the  solution  of  Eq.  6.27  for  the  following  boundary  conditions: 

p(x,H,w)  =  0 


3j3 

3y 


(x,0,u>)  -  0 


•^(°,y,«)  “  "  g  f  j 


(6.37) 


-b, 


and  Pj (x,y,«)  is  the  solution  of  Eq.  6.27  for  the  following  boundary  condi¬ 
tions  : 


The  complex  frequency  response  functions  p*  and  py  are  for  the  hydrodynamic 
pressures  due  to  accelerations  of  the  ground  in  the  horizontal  and  vertical 
directions,  respectively,  and  the  dam  is  rigid;  p^ (x,y,cu)  is  the  correspond¬ 
ing  function  when  the  excitation  is  the  horizontal  acceleration  (jj .  (y)  of  the 

.  th  ^ 

dam  in  the  j  vibration  mode  of  the  associated  dam- foundation  system,  with 

no  motion  of  the  reservoir  bottom;  p.(x,y,td)  is  due  to  vertical  acceleration 

Xj (x)  of  the  reservoir  bottom,  associated  with  the  j  vibration  mode  of  the 

associated  dam- foundation  system,  and  with  no  motion  of  the  dam. 


Solution  of  Eq.  6.27  in  fluid  domains  of  infinite  extent  in  the  up¬ 
stream  direction  (Fig.  6.1)  and  the  boundary  conditions  of  Eqs.  6.35,  6.37 
and  6.36  have  been  reported  earlier  and  are  presented  in  Eqs.  6.39,  6.40  and 
6.41,  respectively: 

X 


-x,  ,  4w  r 

p  (x,y,<D)  -  —  l 
0  it  l 


(-1) 


(2l-l)Jx2t  -  ^ 


exp^x/A^  -  ^  J  cos  A^y 


(6.39) 


-f, 


2w 


pt  (x,y,ui)  -  -  l 
3  9H  £=1 


Jxl-  7 


rH 

in  which  A^  *  (2-£-l)  V2H  and  1^  “  J  ^  (y)  cos  A^y  dy. 


(6.40) 


Pq  (x,y,cu) 


wc  (1+g)  sin 

gw  , ,  U)H 

(1+a)  cos  — 


where  the  reflection  coefficient 

(CrVCW)  -  1 

°  "(Crwr/Cw)  +  1 


(6.41) 


(6.42) 


where  w^  and  Cr  are  the  unit  weight  and  P-wave  velocity  for  the  foundation 
rock  at  the  bottom  of  the  reservoir;  a  =  1  for  rigid  rock  indicating  that 
hydrodynamic  pressure  waves  are  completely  reflected  at  the  reservoir  bottom; 
and  a<l  for  deformable  rock  implying  that  hydrodynamic  pressure  waves  im¬ 
pinging  at  the  reservoir  bottom  are  partially  reflected  back  in  the  fluid 
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domain  and  partially  refracted  into  rock. 


Solution  of  Eq.  6.27  in  fluid  domains  of  infinite  extent  in  the  up¬ 
stream  direction  (Fig.  6.1)  and  the  boundary  conditions  of  Eg.  6.38  may  be 
obtained  by  procedures  similar  to  those  used  in  obtaining  the  above  solu¬ 
tions  of  Eqs.  6.39  and  6.41  [34]  and  by  employing  the  Fourier  transform  with 
respect  to  the  spatial  coordinate  X.  Such  a  general  solution  is  not  neces¬ 
sary,  however,  because  the  resulting  pressures  have  little  influence  on  the 
response  of  the  dam  (see  discussion  preceeding  Eqs.  6.52). 

The  complex  frequency  response  function  for  R^(t),  the  vector  of  total 
hydrodynamic  forces  at  the  upstream  face  of  the  dam,  is,  from  Eq.  6.34, 


R«  (w)  +  V  zf  (oj) 
-°  j-1  3 


K 


~5(w)  +  (co)J 


(6.43) 


in  which  R  (u>)  ,  R  (u>) ,  and  R?(w)  are  the  nodal  forces  statically  equivalent 
-o  -  j  -j 

to  the  corresponding  pressure  functions  at  the  dam  face:  p  (0,y,(i)) , 

f  b  ® 

Pj (0,y,ui) ,  and  p?(0,  y,u) ,  respectively.  They  are  computed  by  using  the 

principle  of  virtual  work  with  the  displacements  between  nodal  points 
defined  by  the  interpolation  functions  used  in  the  finite  element  ideali¬ 
zation  of  the  dam.  Because  the  pressures  act  in  the  horizontal  direction 
on  a  vertical  upstream  face,  elements  of  R^  that  correspond  to  the  y-DOF 
will  be  zero.  Similarly,  from  Eq.  6.34,  the  complex  frequency  response 
function  for  (t) ,  the  vector  of  hydrodynamic  forces  at  the  reservoir 
bottom,  is 


2h(w)  *  2q(u>)  +  l  Z^(u») 
j-1  ■> 


[q*(u»  +  Q*(u)) 


(6.44) 


in  which  2q^'  2j  ^  2j  are  t*ie  nodal  forces  statically  equivalent 

to  the  corresponding  pressure  functions  at  the  reservoir  bottom:  p^(x,0,u)), 

Pj(x,0,oj)  and  p^(x,0,o)),  respectively.  Because  the  pressures  act  in  the 

vertical  direction  on  a  horizontal  bottom  of  the  reservoir,  elements  of  Q 

-n 

that  correspond  to  the  x-DOF  will  be  zero. 


Thus,  by  analyses  of  the  fluid  domain  for  appropriate  conditions,  the 
hydrodynamic  forces  in  Eq.  6.21  have  been  expressed  in  terms  of  the  unknown 

-l 

generalized  accelerations  Z^ (u)  and  by  substituting 


_r>uy. 


n/  2 

Z7"  (CO)  *  -CO  Z  ^  (CO) 


in  terms  of  the  unknown  generalized  displacements. 
6.2.6  Dam-Water-Foundation  System 

The  ntJl  equation  contained  in  Eq.  6.21  is 


(6.45) 


S  «o)Z  <C0)  +  l  S  .  (CO)  Z.<C0) 
nn  n  .  ,  ni  i 

n=l  J 


(6.46) 


'  -l'„  *  (ln,T  ?rqOh(“> 

t  f 

in  which  S  .  (co)  were  defined  in  Eq.  6.22,  L  in  Eq.  6.24,  and  ){)  is  the  sub- 
n]  n  n 

vector  of  containing  only  the  elements  corresponding  to  the  face  (f)  DOF. 

n  -f  -  . 

Expressing  hydrodynamic  forces  R^(co)  and  g^  (to)  in  terms  of  generalized  dis¬ 
placements  from  Eqs.  6.43  -  6.45  leads  to 


S  (co)zf;  «o)  +  l  S  (co)Z^(tt)) 
nn  n  nj  j 

3*n 


=  - 1 *  +  {j£}T(j£  (CO) -co2  l  Z^«0){rJ(CO)  +  r*(w)}J  (6*47) 

-  {4bn>T  lrq<W)?qq<U))  Z^(W){§J(W)+§j(U)}J 

The  set  of  equations  6.47  for  n  ■  1,2  ...  J  may  be  rearranged  and  expressed 
in  matrix  form  as 


SU  (co)  S12(co)  . SlJ^u) 

§21  (co)  §22  (co)  . §2J(co) 


SJ1(W)  SJ2(U) 
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\JL  ' 

Zx  (co) 


z2(w) 


L£  ' 
Lx  (co) 

J, 

VW) 


s„(u)  z*<«»  £**> 


(6.48) 


in  which 


S  .  (u>)  =  S  .  (o>)  +  0) 
n]  n: 


{4T{sJ<">  +  »>>} 

-“2{«bn}T|rq<“)?^t“,{^(“>  +2j'“>} 

valid  for  diagonal  as  well  as  off-diagonal  elements  of  S  (u>) 
and 

Ll(u)  =  -Ll  +  (/]■*  R f(0»  -  U.  )T  Sm(u»S~*(u)Q*(a)) 

-n  n  (,-nJ  -0  [-bnj  -rg  -qq  -0 

Equation  6.49,  after  substitution  of  Eq.  6.22,  becomes 

§  (in)  -  r-<o2+(i+in)X2|  djTm  ^  +  4/t  GL(<o)-u+ n)If(o)  L 

nn  nj  rn-c-n  2n  [_-f  r  J-n 

-“2{tb„}T  ?«S?>  '“><<“>} 

.to)  -  t^[|f +  iu2|*^|t|rJ(w)+5j(u)J 

*2?' <“>} 


and 


(6.49) 


(6.50) 


(6.51a) 


(6.51b) 


Equations  6.48  and  6.49  contain  effects  of  the  foundation  and  water  in  vari¬ 
ous  forms: 

•  Dam- foundation  interaction  effects  appear  in  eigenvalues  and 
eigenvectors  of  the  associated  dam- foundation  system  and  through  the 
foundation  stiffness  matrix  ^(u)  (Eqs.  6.7b  and  6.13). 

/ 

•  Additional  hydrodynamic  loads  R^  on  the  dam  due  to  the  free-field 
ground  motions  computed  on  the  assumption  of  a  rigid  dam. 

•  Dam-water  interaction  effects  appear  through  the  hydrodynamic  forces 
R^  on  the  face  of  the  dam  and  g*  on  the  reservoir  bottom  due  to 
deformational  motions  of  the  dam. 

•  Water- foundation  interaction  effects  appear  through  the  hydrodynamic 
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forces  R.  on  the  face  of  the  dam  and  Q.  on  the  reservoir  bottom  due 
-3  s3 

to  deformational  motions  of  the  reservoir  bottom. 


Dam-water- foundation  interaction  effects  appear  through  the  hydro- 
i  £  b 

dynamic  forces  g^,  g^,  and  at  the  bottom  of  the  reservoir  which 
influences  the  forces  at  the  dam- foundation  interface. 


It  can  be  argued  and  shown  through  numerical  examples  that  several  terms 
in  Eqs.  6.49  -  6.51  are  relatively  small  and  can  be  dropped  without  intro¬ 
ducing  any  significant  errors.  One  group  of  such  terms  arises  from  the  hydro 


hydrodynamic  forces  at  the  upstream  face  of  the  dam  due  to  deformational 
motions  of  the  reservoir  bottom. 


Dropping  these  terms  from  Eqs.  6.49  -  6.51  leads  to  their  simplified 
version : 


in  which 

_  _/> 

SttirtZMGj)  ■  L  (ai) ,  b*x,  y 

(6.52) 

=  +  fnf|f(“>)“(l+in)!f(0)]in 

«2{4t 

(6.53a) 

snj  <00 

-  j*||f(«i»-(i+in>If  (0)]^  +  u)2{^jT  Ij(u) 

(6.53b) 

E*(o0 

-  -64  *  {4}T  5o'“> 

(6.53c) 

n 

These  are  the  final  equations  governing  the  response  of  the  dam  to  horizon¬ 
tal  or  vertical  ground  motions,  applied  separately.  Included  in  these  equa¬ 
tions  are  the  hydrodynamic  effects  and  dam- foundation  interaction  effects 
that  are  significant  in  the  response  of  the  dam.  Equation  6.52  represents 
J  algebraic  equations  in  the  unkown  complex  frequency  response  functions  in 

the  generalized  coordinates  (co) ,  j  -  1,2, - J,  corresponding  to  the 

J  Ritz  vectors  included  in  the  analysis.  The  coefficient  matrix  S(u>)  is  to 
be  determined  for  each  excitation  frequency  and  simultaneous  solution  of  the 
equations  leads  to  the  values  of  5*  at  that  excitation  frequency.  Repeated 
solution  for  several  hundred  excitation  frequencies,  covering  the  frequency 

range  over  which  the  ground  motion  and  structural  response  have  significant 

—l 

components,  leads  to  the  complete  functions  Z,  (<*>) . 
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The  number  of  Ritz  vectors  that  need  be  included  in  the  analysis 
depends  on  the  system  properties,  the  response  quantities,  and  frequency 
range  of  interest.  The  number  necessary  in  analysis  of  concrete  gravity 
dams  is  determined  in  Sec.  7.6.  As  will  be  seen  later,  a  few  generalized 
coordinates  (ten  to  twenty)  are  sufficient  in  the  analysis  of  dams  ideal¬ 
ized  as  finite  element  systems  with  say  200  DOF,  thus  leading  to  very 
efficient  procedures. 

6.2.7  Singularities  of  Response 

Elements  of  matrix  S  (a))  and  load  vector  are  to  be  computed  for  a 

frequency  u>  by  using  Eq.  6.53  together  with  the  pressure  functions  p0  and 

p.  of  Eqs.  6.39  and  6.40.  The  pressure  functions,  and  hence  the  corres- 

3  -l  r  r 

ponding  nodal  forces  Rn  and  R„  are  unbounded  at  oj  =  oj  where  ui  =  (2m-l)TC/2H, 

“3  mm 

is  the  m  resonant  frequency  of  the  fluid  domain.  Consequently,  the 
elements  of  S(w)  and  Lx(oj)  are  unbounded  at  these  frequencies. 

When  J  =  1,  i.e.  when  only  one  vibration  mode  of  the  dam  is  considered 

in  the  analysis,  the  J  equations  of  Eq.  6.53  reduce  to  one  equation  and  the 
r 

response  at  u  *  can  be  obtained  through  a  limiting  process.  However, 

when  J  >  1,  i.e.  when  more  than  one  mode  is  consider eed  in  the  analysis, 

the  limiting  process  yields  a  system  of  equations  such  that  J  (co)  is  singular 

at  0)  =  u)  .  In  particular,  all  the  M  equations  become  identical  to  each 
m 

other  and  no  solution  can  be  obtained. 

This  degeneracy  of  the  equations  may  be  considered  as  a  limitation  of 

the  substructure  method  of  analysis.  However,  this  limitation  is  of  no 

practical  consequence  in  obtaining  earthquake  responses.  Numerical  values 

for  the  frequency  responses  may  be  obtained  for  values  of  u)  arbitrarily 
T  T 

close  to  w  .  The  singularities  at  u)  constitute  a  discrete  set  and  will 
m  m 

therefore  not  affect  the  values  of  the  Fourier  integrals  which  will  lead  to 
earthquake  responses  (Sec.  6.3). 

6.3  Response  to  Arbitrary  Ground  Motion 

-X  — y 

Once  the  complex  frequency  response  functions  Z  (u>)  and  Z  (oj)  ,  n  *  1,2, 

n  n 

...,  J,  have  been  determined  by  solving  Eq.  6.52  for  an  appropriate  range  of 
excitation  frequency  w,  the  responses  to  arbitrary  ground  motion  can  be  ob¬ 
tained  as  the  superposition  of  responses  to  individual  harmonic  components 
of  the  excitation  through  the  Fourier  integral: 
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ZZlt)  =  ~  |  ZZ0ti)  AZ  (w)  eil0tda>;  Z  *  x,y 

n  z»  j  n  g 


Z  Z 

in  which  A  (<»))  is  the  Fourier  transform  of  a  (t)  : 

g  g 


A*  to) 


fd 

j  a^(t)e  i0Jtdt;  Z  *  x,y 


(6.54) 


(6.55) 


where  d  is  the  duration  of  ground  motion. 

The  combined  response  Z  (t)  to  horizontal  and  vertical  components  of 

n 

ground  motions  acting  simultaneously  is 


Z  (t)  =  ZX(t)  +  Z*(t) 
n  n  n 


(6.56) 


Repeating  this  procedure  for  all  the  necessary  values  of  n,  the  displacement 
response  is  obtained  by  transforming  back  from  generalized  to  nodal  point 


coordinates 


rjt)  »  I 
n*l 


(6.57) 


The  stresses  0p(t)  in  finite  element  p  at  any  instant  of  time  are  re¬ 
lated  to  the  nodal  point  displacements  r^(t)  for  that  element  by 


0  1 t)  -  T  r  (t) 

-p  -p  -p 


(6.58) 


where  T  is  the  stress  transformation  matrix  for  the  finite  element  p.  At 
— P 

any  instant  of  time,  the  stresses  throughout  the  dam  are  determined  from  the 
nodal  point  displacements  by  applying  the  transformations  of  Eq.  6.58  to 
each  finite  element. 

6.4  Static  Stress  Analysis 

The  standard  procedure  for  analysis  of  stresses  due  to  static  loads 
(Chapter  3)  applies  to  structures  on  rigid  base.  Including  foundation  flex¬ 
ibility,  the  equilibrium  equations  for  a  dam  subjected  only  to  hydrostatic 
pressures  and  forces  associated  with  its  own  weight  can  be  obtained  as  a 
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special  case  of  the  equations  of  dynamics  (Eq.  6.11) 


- 1 

S' 

_ 1 

r 

’  R 

* 

*bb  +  4f<0>_ 

Sb 

0 

(6.59) 


The  vector  of  static  loads  is  denoted  by  R  and  the  effects  of  foundation 
deformations  are  represented  by  the  foundation  stiffness  matrix  at  zero 
frequency. 


In  principle,  Eq.  6.59  may  be  directly  solved  for  the  desired  nodal 
point  displacements  resulting  from  the  gravity  loads.  In  practice,  this 
direct  approach  may  be  unsatisfactory  because  the  displacements  r  and  rfe 
include  large  rigid  body  components  which  have  no  effect  on  the  stresses 
in  the  dam.  To  avoid  this  difficulty,  the  total  displacements  are  expressed 
as  the  sum  of  rigid  body  displacments  plus  an  increment  due  to  deformations, 

1  •  6  m  / 


r 

r 

-0 

r 

.  = . 

•  +  • 

-rbj 

Sabj 

-b; 

(6.60) 


in  which  the  subscript  refers  to  rigid  body  displacements,  and  the 
"tilde"  identifies  the  relative  displacement  components.  The  rigid  body 
and  relative  displacements  at  the  base  of  the  dam  are  shown  in  Fig.  6.3. 

Introducing  Eq.  6.60  into  Eq.  6.59  yields: 

'5  Sb  ' 

%,T  Sbb  +  Sf ,0) 

But  since  the  rigid  body  displacements  produce  no  forces  in  the  structure, 
this  may  be  written  as 


/  A 

*0 


l -Ob) 


S  Sb 

VT  ^bb  +  ^f(0) 


r  ~ 

r 


v£b. 


(6.61) 


0  $f(0) 


'  r  ' 


-Ob 


>  + 


*  Sb 

%/  ^b+?£<0' 
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R  ' 
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*  s  • 

o 

(6.62) 
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Designating  the  foundation  forces  associated  with  rigid  body  displacements 
by  R Qb ,  where 


-Ob  =  -f(0)  -Ob 


(6.63) 


Eq.  6.62  may  be  written  as  follows: 


5b  5bb  +  -f(0)  4 


(6.64) 


Equation  6.64  thus  provides  a  means  to  solve  directly  for  the  relative 
displacements;  the  rigid  body  displacements  enter  only  in  the  evaluation  of 
the  rigid  body  foundation  forces  RQb .  To  evaluate  these  forces,  it  is  con¬ 
venient  to  express  the  rigid  body  displacements  r0b  in  terms  of 


-g  =  cg 


the  three  components  of  displacement  of  the  center  of  gravity  (c.g)  of  the 
base  (Fig.  6.4): 


r  *  D  r 
-Ob  -  -eg 


(6.65) 


where  the  transformation  matrix  D  involves  only  the  geometry  of  the 
location  of  the  base  nodal  points  and  the  c.g  of  the  base. 

Similarly,  the  resultant  base  forces  R  corresponding  with  the  dis¬ 
placements  rcg  can  be  expressed  in  terms  of  the  nodal  forces  R^  as  follows: 


R  =  D  R 
-eg  -  -Ob 


(6.66) 


where 


FIG.  6.3  RIGID  BODY  AND  RELATIVE  DISPLACEMENTS  AT  BASE  OF  THE  DAM 


6 


in  which 


U  “  resultant  horizontal  force 
V  =  resultant  vertical  force 

M  =  resultant  moment  about  the  center  of  gravity 
Introducing  Eqs.  6.63  and  6.65  into  Eq.  6.66  yields 

hg  ■  ST  -Ob  "  ?T  S-C l0)  Eob  •  2T  If  l0>5  Eog 

from  which 


R 

-eg 


r 

-eg 


(6.67) 


if  the  symbol 

|f0  =  DT  Sf (0)  D  (6.68) 

is  used  to  represent  the  rigid  body  foundation  stiffness  matrix. 

Solving  Eq.  6.67  for  the  rigid  body  displacements: 

r  =  R 

-eg  -fO  -eg 

and  making  use  of  Eqs.  6.65  and  6.63,  the  nodal  forces  R^  associated  with 
rigid  body  displacements  are  expressed  in  terms  of  R  ,  the  resultant  forces 
at  the  base  due  to  all  applied  loads  R: 


R-, 

-Ob 


=  Sf  (0) 


5  C 


-eg 


(6.69) 


With  this  result  the  right  hand  side  of  Eq.  6.64  is  known.  This  equation 
can  be  solved  for  the  relative  nodal  displacements,  from  which  all  desired 
stresses  in  the  dam  cam  be  found,  without  any  specific  consideration  of  the 
rigid  body  displacements. 
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6.5  Special  Cases 


6.5.1  Dam- Foundation  System,  No  Water 

The  frequency-domain  equations  for  the  dam,  without  any  impounded 

water  and  including  dam- foundation  interaction  effects,  can  be  obtained 

from  the  general  Eqs.  6.52  and  6.33  by  simply  dropping  the  hydrodynamic 
L  t 

terms  R^  (u>)  and  R .  (to) .  The  resulting  equations  are 


Stto^to)  =  (to);  l  =  x,y 


(6.70) 


in  which 


(to)  =  |-co2+(l+in)A2 ip  +  (to)- (l+in)£-(0)li|> 

i  nj  -n-c-n  -“-n  |_-t  -x  j  -n 


(6.71a) 


S  .(to) 
nj 


jL(io)  -  (l+in)i,(0) 
•*-n  — r  -r 


(6.71b) 


Tl  .  T 

L  =  -lit  m  1 
n  £n-c-c 


(6.71c) 


Equations  6.70  and  6.71  are  the  same  as  those  presented  in  Chapter  5  for 
dynamic  response  of  dams  including  structure-foundation  interaction. 

6.5.2  Dam  on  Rigid  Foundation,  with  Water 

In  this  case,  the  base  DOF  can  be  eliminated  from  the  analysis;  the 
associated  structure- foundation  system  defined  in  Sec.  6.2.4  is  simply  the 
dam  on  rigid  foundation.  Denoting  the  mass  and  stiffness  matrices  of  the 
dam  on  rigid  base  by  m  and  k,  the  natural  vibration  frequencies  u)^  and 
corresponding  mode  shapes  <£n  of  this  system  are  solutions  of  the  eigen¬ 
value  problem 


kd  =  a)  md 
-xn  n--n 


(6.72) 


The  structural  displacements  are  expressed  as  a  linear  combination  of  the 
first  J  mode  shapes: 


r*(t)  =  l  Y*(t>$.  and  r^(co)  =  £  Y*(<i))$. 

■j=l  3  3  i=l  3  3 


(6.73) 


Thus,  the  frequency  domain  equations  for  the  dam  on  a  rigid  foundation, 
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including  hydrodynamic  effects/  can  be  obtained  by  appropriately  speciali¬ 
zing  Eqs.  6.52  and  6.53.  The  resulting  equations  are 


S((D) 


Y^Cid) 


I^(u>) 


Z  =  x,y 


(6.74) 


in  which 

S  (ui)  *  |-a)2+(l+in)u  +a)2Uf}TRf (w) 

nn  (_  nj  *n-*n  £n  -n 

S  .(w)  =  u>2(4>f >T5f  (u) 
n]  -n  -3 

LZ(u>)  =  -<j>Tml£  +  {<pf}T^;  (w) 
n  -n-  —  -n  -u 


(6.75a) 

(6.75b) 

(6.75c) 


These  equations  are  equivalent  to  those  in  Chapter  4  (Eqs.  4.16-4.18 
and  4.34-4.35)  developed  directly  for  dams  with  water  on  rigid  foundation. 


6.5.3  Dam  on  Rigid  Foundation,  No  Water 

Equations  governing  conplex  frequency  responses  of  a  dam  without  water 

are  obtained  by  further  specializing  Eqs.  6.74  and  6.75,  simply  by  setting  to 

-L 

zero  the  hydrodynamic  forces  R_(«J)  and  R .  (w) .  The  coefficients  of  Eq.  6.74 
are  then  given  bys 

S  (w)  =  -m2+(l+in)W2  UTm<f>  (6.76a) 

nn  I  nj  ■Ln--n 


S  .«*»  =  0  (6.76b) 

nj 

Z  t  Z 

u  (W)  =  -  *  m  l  (6.76c) 

n  £n — 


The  matrix  S (oj)  is  now  diagonal  since  the  chosen  Ritz  vectors  are  the  natural 
inodes  of  vibration  of  the  system  considered. 

Comparing  Eqs.  6.76  and  6.75,  it  is  apparent  that  §  (o>)  and  £>(oj)  are 
modified  due  to  the  presence  of  water.  The  diagonal  terms  of  S(ui)  are  modi¬ 
fied  by  an  additional  frequency-dependent  mass.  Off-diagonal  mass  terms 
appear  because  the  Ritz  vectors,  the  natural  modes  of  vibration  of  the  dam 
(without  water) ,  are  not  the  vibration  modes  of  the  combined  dam-water  system, 
and  they  become  coupled  due  to  hydrodynamic  interaction  effects.  An  additional 
l 

term  appears  in  L  (co)  due  to  hydrodynamic  loads. 
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These  equations  are  equivalent  to  those  in  Chapter  3  (Eq.  3.11)  derived 
directly  for  the  dam  without  water  on  rigid  foundation. 

6.6  Computer  Program 

Based  on  the  analytical  procedures  developed  in  this  chapter,  a  com¬ 
puter  program  has  been  written  in  FORTRAN  IV  to  numerically  evaluate  the 
responses  of  concrete  gravity  dams,  including  various  effects  of  both  the 
water  and  the  foundation.  These  include  effects  arising  from  interaction 
between  the  dam  and  foundation,  dam  and  water,  water  and  foundation,  and 
from  interaction  among  all  three  substructures — dam,  water,  and  foundation. 
The  dynamic  stiffness  matrix  for  the  foundation  is  not  computed  within  the 
computer  program.  Thus  it  should  be  computed  externally  and  supplied  to 
the  computer  program.  Available  data  and  methods  132]  may  be  utilized  to 
determine  the  dynamic  stiffness  matrix  for  a  foundation  region  idealized  as 
a  viscoelastic  half  plane. 

The  computer  program  is  capable  of  analyzing  the  response  of  a  dam  for 
four  conditions:  Dam  on  rigid  foundation  excluding  hydrodynamic  effects, 
dam  on  rigid  foundation  including  hydrodynamic  effects,  dam  on  flexible 
foundation  excluding  hydrodynamic- effects,  and  dam  on  flexible  foundation 
including  hydrodynamic  effects. 

The  response  of  dam-water- foundation  systems ,  idealized  as  described 
in  Chapter  2,  to  horizontal  and  vertical  components  of  ground  motion  can  be 
analyzed.  The  dam  is  treated  as  am  assemblage  of  two-dimensional  finite 
elements;  the  impounded  water  and  foundation  are  treated  as  continua. 

The  output  from  the  computer  program  includes  the  complex  frequency 
response  functions  describing  the  response  to  harmonic  ground  motions  and 
the  complete  time-history  of  displacements  and  stresses  throughout  the  dam. 
Because  the  program  is  capable  oi  including  any  number  of  Ritz  shape  vectors, 
results  can  be  obtained  to  any  desired  degree  of  accuracy. 

The  users  guide  and  listing  of  the  program  are  included  in  Appendices 


A  and  B. 


7.  PRELIMINARY  RESULTS  FOR  PARAMETER  SELECTION 


7.1  Scope  of  the  Chapter 

Several  preliminary  analyses  were  carried  out  to  assist  in  the  selection 
of  parameters  defining  the  system  and  also  those  necessary  to  carry  out  the 
analyses  for  the  dynamic  response  study  presented  in  Chapter  8.  Results 
from  two  analyses  of  a  dam  foundation  system,  assuming  the  system  to  be  in 
plane  stress  in  one  and  in  plane  strain  in  the  other,  are  presented  to  pro¬ 
vide  a  basis  for  selecting  one  of  the  two  assumptions,  similarly,,  response 
of  the  dam  is  presented  and  compared  for  two  damping  models,  viscous  damping 
and  constant  hysteretic  damping.  With  the  aid  of  response  results,  the 
minimum  number  of  parameters  —  those  that  have  the  most  important  influence 
on  the  dam  response  —  defining  the  elastic  modulii  of  the  dam  and  founda¬ 
tion  are  identified.  The  number  of  generalized  coordinates  necessary  to 
obtain  accurate  responses  over  the  desired  range  of  excitation  frequency  is 
determined  with  the  aid  of  response  results  obtained  with  varying  numbers  of 
these  coordinates. 

7.2  Systems,  Ground  Motions,  and  Response  Quantities 
7.2.1  Systems 

The  dam  is  idealized  as  having  a  triangular  cross  section  with  a  ver¬ 
tical  upstream  face  and  a  downstream  face  with  a  slope  of  0. 8  :1 .  The  dam 
is  assumed  to  be  homogeneous  with  linearly  elastic  and  isotropic  proper¬ 
ties  for  mass  concretes  Young's  modulus  Es  *  2  or  4  million  psi;  Pois¬ 
son's  ratio  *  0.2;  and  unit  weight  ■  155  lbs/ft3.  Representative  values 
for  the  latter  two  parameters  have  been  chosen.  The  possible  small  varia¬ 
tion  around  the  chosen  values  will  have  little  influence  on  the  dam  res¬ 
ponse.  Two  damping  assumptions  are  considered:  Viscous  danping  with 
danping  ratio  £g  *  0.05  in  each  natural  mode  of  vibration  of  the  dam 
alone  (without  water)  on  a  fixed  base,  and  constant  hysteretic  damping 

with  the  energy  loss  coefficient  n  ■  0.1. 

s 

The  finite  element  idealization  for  a  monolith  of  the  dam  (Fig.  7.1) 
consists  of  20  quadrilateral  elements  and  26  nodal  points  which  provide 
42  degrees  of  freedom  on  a  fixed  base  and  52  on  a  flexible  base. 

The  dam  monolith  is  supported  on  the  surface  of  a  homogeneous,  isotro- 
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1  FINITE  ELEMENT  IDEALIZATION  OF  A  DAM  MONOLITH  WITH  TRIANGULAR 
CROSS-SECTION 
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pic i  linearly  viscoelastic  half  space.  The  following  properties  are  assumed 

for  the  foundation  materials:  E_  *  2E  ;  Poisson's  ratio  »  1/3?  and  unit 

£  s 

weight,  wf  =  165  lbs/  ft3.  The  latter  two  parameter  values  are  not  varied 

because  the  response  is  insensitive  to  the  small  variations  in  those 

parameter  values  that  cam  occur  in  practical  problems.  The  foundation 

soil  or  rock  is  idealized  as  a  constant  hysteretic  solid  with  several 

values  of  energy  loss  coefficient  nf  ■  0.1,  0.25,  and  0.5,  or  as  a  Voigt 

solid  with  energy  loss  coefficient  *  0.1.  Impounded  water  in  the 

reservoir  has  a  constant  depth  H  and  is  idealized  as  extending  to  infinity 

in  the  upstream  direction.  The  unit  weight  of  water  w  *  62.5  lbs/ft3  and 

the  velocity  of  sound  in  water  ■  4720  ft/ sec.  Two  values  of  the  depth  of 

water  H,  relative  to  the  height  of  the  dam  H  ,  are  considered:  H/H  =  0 

s  s 

(no  water)  and  1  (full  water) . 

The  dam  and  foundation  both  are  assumed  to  be  in  a  state  of  generalized 
plane  stress  or  plane  strain. 

7.2.2  Ground  Motions 

The  excitation  for  the  dam-water- foundation  system  is  defined  by  the 

two  components  of  free-field  ground  motion  in  the  plane  of  a  monolith  of 

the  dam:  the  horizontal  component  transverse  to  the  dam  axis?  and  the 

vertical  component .  Each  component  of  ground  acceleration  is  assumed  to 

be  harmonic,  aX(t)  *  ay(t)  =  e1U)t,  with  the  excitation  frequency  to  be 
9  9 

varied  over  a  relevant  range. 

7.2.3  Response  Quantities 

The  response  of  the  dam  to  harmonic  free  field  ground  motion  was  de¬ 
termined  by  the  analytical  procedures  presented  in  Chapter  6.  Five  mode 
shapes  on  fixed  base  or  generalized  coordinates,  as  appropriate,  were 
included  in  the  analysis.  Although  the  analysis  led  to  response  at  all 
nodal  points,  only  the  results  for  a  few  selected  degrees-of-freedom  are 
presented.  For  dams  on  rigid  foundations,  the  complex  frequency  response 
functions  for  horizontal  acceleration  at  the  crest  of  the  dam  are  presented. 
In  addition,  the  conplex  frequency  response  functions  for  the  horizontal 
acceleration  at  the  center  of  the  base  of  the  dam  and  for  vertical  acceler¬ 
ations  at  the  upstream  and  downstream  edges  of  the  base  are  presented  when 
the  interaction  between  the  structure  and  flexible  foundation  material  is 

included. 
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The  absolute  value  (or  modulus)  of  the  complex -valued  acceleration 
response  is  plotted  against  normalized  excitation  frequency  parameter 
co/uj^,  where  =  fundamental  fixed  base  natural  frequency  of  the  alone. 

7.3  Comparison  of  Plane  Strain  and  Plane  Stress  Assumption 

Because,  at  the  present  time,  the  analysis  procedure  presented  in 
Chapter  6  has  been  implemented  only  for  two-dimensional  systems,  its  appli¬ 
cation  is  restricted  to  systems  in  plane  stress  or  plane  strain.  Although, 
as  mentioned  in  Chapter  2,  neither  of  the  two  models  are  strictly  valid, 
the  former  is  better  for  the  dam  and  the  latter  for  the  foundation.  However, 
in  order  to  define  the  dam- foundation  system  on  a  consistent  basis,  the 
same  model  should  be  enployed  for  both  substructures. 

The  results  of  two  analyses  of  the  dam- foundation  system  described  in 
Section  7.2.1,  assuming  the  system  to  be  in  plane  stress  in  one  and  in  plane 
strain  in  the  other,  are  presented  in  Fig.  7.2.  Constant  hysteretic  damping 
with  energy  loss  coefficient  r)  *  0.1  was  assumed  for  the  dam  and  foundation 
materials.  The  results  presented  are  for  systems  with  Ef/Es  *  2,  indepen¬ 
dent  of  the  Ef  and  Eg  value  separately.  It  is  obvious  from  Fig.  7.2  that 
the  response  of  the  dam,  including  effects  of  dam- foundation  interaction, 
is  essentially  the  same  under  two  assumptions  —  plane  stress  or  plane 
strain. 

The  dam-foundation  system  is  assumed  to  be  in  plane  stress  for  all 
the  results  in  this  and  later  chapters.  As  discussed  in  Chapter  2,  based 
on  physical  behavior,  this  assumption  is  appropriate  for  describing  the 
behavior  of  a  concrete  gravity  dam  vibrating  at  large  anplitudes  of  motion, 
but  not  for  the  foundation.  However,  based  on  the  above  results,  the  same 
assumption  may  be  made  for  the  foundation  behavior  without  significant  loss 
in  accuracy. 

7.4  Elastic  Moduli  Parameters 

The  complex  frequency  acceleration  response  functions  describe  the 

steady-state  acceleration  response  of  dams  to  harmonic  excitations.  These 

functions  for  a  dam  on  a  rigid  foundation  and  no  water,  plotted  against 

normalized  excitation  frequency  parameter  w/u^,  apply  to  dams  of  all  heights 

provided  they  all  have  the  same  cross-sectional  shape.  Furthermore,  they 

are  independent  of  the  unit  weight  w  and  Young's  modulus  E  for  the  dam 

s  s 
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RESPONSE  OF  A  DAM  WITHOUT  WATER,  INCLUDING  DAM-FOUNDATION 
INTERACTION, TO  HARMONIC  HORIZONTAL  GROUND  ACCELERATION 
OBTAINED  UNDER  TWO  ASSUMPTIONS  FOR  THE  SAME  SYSTEM: 

(I)  PLANE  STRESS:  AND  (2)  PLANE  STRAIN 


concrete.  It  can  be  concluded  from  the  equations  of  motion  for  the  dam, 

including  the  dam- foundation  interaction  but  excluding  hydrodynamic  effects, 

that  the  complex  frequency  response  functions  plotted  against  normalized 

excitation  frequency  do  not  depend  separately  on  Eg  and  but  only  the 

ratio  E,/E  .  These  functions  for  dams  on  rigid  foundations  vary  signi- 
I  s 

ficantly  with  Eg  when  hydrodynamic  effects  are  considered  [22].  They  vary 

even  for  dams  on  flexible  foundations:  for  a  fixed  value  of  E_/E  ,  the 

r  s 

value  of  Eg  has  influence  on  the  response  of  the  dam  including  dam-water 
interaction  effects.  Figure  7.3  demonstrates  that  this  influence  is  minor 
in  the  response  to  horizontal  ground  motion  but  significant  in  the  response 
to  vertical  ground  motion. 

7. 5  Damping  Models  and  Parameters 

7.5.1  Dam 

The  analytical  procedures  presented  in  Chapters  3-6  included  two 
models  —  viscous  damping  and  constant  hysteretic  dancing  —  for  energy 
dissipation  in  the  dam.  Viscous  damping  had  been  selected  earlier  [22], 
in  studying  hydrodynamic  effects  in  earthquake  response  of  dams.  Constant 
hysteretic  dancing  is  preferable  for  conceptual  as  well  as  computational 
reasons,  especially  when  structure-soil  interaction  is  included  in  the 
analysis  [29].  Therefore,  all  numerical  results  presented  in  Chapters  8 
and  9  are  for  dams  with  constant  hysteretic  damping.  In  order  to  provide 
a  basis  for  comparing  these  results  with  those  presented  earlier  [22] , 
response  results  using  both  damping  models  are  presented  in  this  section 
and  compared. 

Using  the  analysis  procedure  in  Chapter  6,  the  steady-state  response 
of  the  idealized  dam  monolith  (Fig.  7.1)  to  harmonic  ground  acceleration 
in  the  horizontal  direction  is  determined  for  several  cases.  Assuming 
the  concrete  properties  listed  in  Section  7.2.1,  with  Young's  modulus  * 

4  x  106  psi,  and  considering  the  foundation  to  be  rigid,  the  complex 
frequency  response  functions  for  the  dam  were  determined  for  the  following 
four  cases.  (The  first  five  vibration  modes  were  included  in  the  analysis.) 

1.  Viscous  damping,  Cs»  0.05  for  all  vibration  modes;  no  water 

2.  Constant  hysteretic  damping,  ns=  0.10?  no  water 

3.  Viscous  damping,  £g*  0.05  for  all  vibration  modes;  full  water 

4.  Constant  hysteretic  damping,  r)s=  0.10?  full  water 
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FIG.  7.3  INFLUENCE  OF  YOUNG'S  MODULUS  Es  OF  DAM  CONCRETE  ON  RESPONSE 
OF  A  DAM,  INCLUDING  EFFECTS  OF  WATER  AND  FOUNDATION,  TO 
HARMONIC  GROUND  MOTION.  Ef/Es  IS  THE  SAME  IN  BOTH  CASES. 
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The  selected  values  of  damping  factor  n  in  constant  hysteretic  damping 

s  # 

and  viscous  damping  ratio  £sfor  all  modes  of  vibration  are  related  by  ns= 

2CS.  As  a  result,  the  energy  dissipated  per  vibration  cycle  in  each  natural 
mode  of  vibration  at  the  resonant  excitation  frequency  would  be  the  same  for  • 

the  two  damping  models.  The  energy  dissipated  per  cycle  in  constant  hystere¬ 
tic  damping  is  independent  of  excitation  frequency,  but  varies  linearly  with 
frequency  for  viscous  damping. 

Because  the  energy  dissipated  at  the  fundamental  resonant  frequency  of 
the  dam  alone  is  the  same  for  the  two  damping  models  and  the  selected  para¬ 
meters,  the  dam  response  is  independent  of  the  type  of  damping  (Fig.  7.4). 

At  the  fundamental  resonant  frequency  of  the  dam,  including  hydrodynamic  ^ 

effects,  the  energy  dissipated  in  constant  hysteretic  damping  will  be  the 
same  but  reduced  in  viscous  damping,  compared  to  the  energy  dissipated  at 
the  natural  frequency  of  the  dam  alone.  As  a  result,  the  resonant  response  • 

is  significantly  smaller  for  the  dam  with  constant  hysteretic  damping,  com¬ 
pared  to  the  results  for  viscous  damping  (Fig.  7.4) .  The  differences  in 
resonant  responses  for  the  two  damping  models  are  directly  related  to  the 
reduction  in  resonant  frequencies  due  to  dam-water  interaction  effects. 

Because  such  reduction  is  relatively  small  for  higher  resonant  frequencies,  a 

the  corresponding  resonant  responses  are  affected  less  by  the  damping  model 
(Fig.  7.4) .  At  excitation  frequencies  not  close  to  resonant  frequencies, 
damping  has  little  effect  on  the  response  and  it  is  essentially  the  same  with  * 

the  two  damping  models. 

Constant  hysteretic  damping  with  n  =  0.1  is  assumed  for  the  dam  in  ob- 

s 

taining  all  subsequent  results  in  Chapters  7,  8,  and  9. 

7.5.2  Foundation 

Numerical  results  and  analytical  procedures  have  been  presented  to 
determine  dynamic  stiffness  matrices  for  a  viscoelastic  half-plane  of  Voigt 
or  constant  hysteretic  materials  [32].  In  this  section,  computed  responses  of 
the  dam  including  structure-soil  interaction  effects  are  presented  to  aid  in 
selecting  the  damping  model  and  parameters  for  the  foundation  material. 

The  response  of  the  idealized  dam  monolith  (Fig.  7.1)  without  water  but 

including  dam- foundation  interaction  effects  is  computed  by  the  procedures  of 

Chapter  6.  The  following  properties  are  selected  for  the  dam  in  addition  to 

those  listed  in  Section  7.2.1:  Young's  modulus  Eg  =  4  x  106  psi,  constant 

hysteretic  damping  factor  n  =  0.1.  The  foundation  modulus  E  is  twice  that 

s  f 

of  the  dam  and,  in  addition  to  the  properties  of  the  foundation  matarial  listed 
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CURVE  DAMPING  MODEL  WATER 

1  VISCOUS  NO 

2  CONSTANT  HYSTERITIC  NO 

3  VISCOUS  FULL 

4  CONSTANT  HYSTERITIC  FULL 
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in  Section  7.2.1,  the  following  are  selected:  E^/E^  =  2,  energy  loss  coef¬ 
ficients  =  0.1  for  the  Voigt  solid;  and  nf  =  0.1  for  the  constant  hystere- 
tic  solid.  The  two  dam- foundation  system^  analyzed  have  identical  properties 
with  one  exception:  the  foundation  is  of  Voigt  material  in  one  case,  and  con¬ 
stant  hysteretic  material  in  the  other  case.  Including  five  generalized  co¬ 
ordinates  in  the  analysis,  the  complex  frequency  response  functions  for  the 
dam  were  obtained  for  the  two  systems.  The  results  presented  in  Fig.  7.5a 
demonstrate  that  the  damping  model  for  the  foundation  has  little  effect  on 
the  resonant  frequencies  of  the  dam  or  on  the  response  except  at  the  fundamental 
resonant  frequency.  The  fundamental  resonant  response  of  the  dam  is  somewhat 
larger  with  Voigt  foundation  material. 

For  a  Voigt  solid,  the  energy  loss  per  cycle  of  harmonic  vibration  is 
proportional  to  the  excitation  frequency.  However,  over  a  considerable  range 
of  frequencies,  rocks  and  soils  exhibit  energy  loss  essentially  independent 
of  the  frequency  of  vibration.  Such  materials  are  better  idealized  as  con¬ 
stant  hysteretic  solids.  This  is  the  model  adopted  for  the  foundation  mater¬ 
ial  for  the  dam  responses  presented  in  Chapters  8  and  9.  Based  on  the  above 
discussion,  however,  those  results  are  also  indicative  of  the  response  of 
dams  on  Voigt  foundation  materials. 

Using  the  properties  mentioned  above  for  the  dam  and  foundation,  the 
analysis  was  repeated  for  two  other  values  of  damping  coefficient  for  a  con¬ 
stant  hysteretic  foundation:  nf  =  0.25  and  0.50.  The  complex  frequency 
response  functions  for  the  dam  with  constant  hysteretic  damping,  ng  =0.1, 
on  a  constant  hysteretic  foundation  are  presented  in  Fig.  7.5b  for  three 
values  of  the  damping  factor  =  0.1,  0.25,  and  0.50.  The  responses  are 
insensitive  to  the  foundation  damping  except  in  a  neighborhood  of  the  fun¬ 
damental  resonant  frequency.  Even  this  resonant  response  is  not  influenced 
greatly  by  damping;  increasing  =  0.1  to  0.25  reduces  the  resonant  response 
by  only  13%.  Damping  in  the  foundation  is  modeled  as  a  constant  hysteretic 
solid  with  nf  =  0.1  in  obtaining  all  subsequent  results  in  Chapters  7,  8, 
and  9.  The  assumed  value  of  damping  factor  is  approp.-i  te  for  sound  rocks 
at  sites  of  major  concrete  dams.  Because  the  dam  respunse  is  not  very  sen¬ 
sitive  to  foundation  damping,  however,  those  results  are  approximately  indi¬ 
cative  of  tne  dam  response  with  other  levels  of  foundation  damping. 
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7.6  Number  of  Generalized  Coordinates 

The  governing  equations  in  the  frequency  domain  were  transformed  from 
the  physical  nodal  point  displacement  coordinates  to  generalized  coordinates 
(Chapters  3-6) .  As  many  of  these  generalized  coordinates  should  be  included 
in  the  analysis  as  necessary  to  obtain  accurate  responses  over  a  desired 
range  of  excitation  frequencies.  The  required  number  of  generalized  coor¬ 
dinates  will  generally  be  a  small  fraction  of  the  number  of  degrees-of- 
freedom  in  the  finite  element  idealization  of  the  dam. 

The  natural  modes  of  vibration  of  the  dam  alone  were  selected  as  the 
generalized  coordinates  for  analysis  of  the  dam  on  rigid  foundation,  with 
or  without  water  (Chapters  3  and  4) .  In  order  to  accurately  obtain  the 
response  of  the  dam  without  water  on  a  rigid  base,  all  the  modes  of  vibra¬ 
tion  should  be  included  which  have  natural  frequencies  within  and  close  to 
the  range  of  excitation  frequenices  of  interest,  and  contribute  signifi¬ 
cantly  to  the  response  in  this  frequency  range.  Accelerograph  records 
accurately  reproduce  ground  motion  components  with  frequencies  up  to  approx¬ 
imately  25  cps.  The  first  few,  say  five,  natural  vibration  modes,  would 
generally  be  sufficient  to  obtain  accurate  responses  of  concrete  gravity 
dams  —  they  have  relatively  high  natural  frequencies  —  in  this  frequency 
range.  Dam-water  interaction  significantly  reduces  the  fundamental  reso¬ 
nant  frequency  but  has  little  influence  on  the  higher  resonant  frequencies. 
Thus,  usually,  the  same  number  of  natural  vibration  modes  would  be  sufficient 
for  analysis  of  the  response  of  the  dam,  with  or  without  water. 

The  eigenvectors  (vibration  mode  shapes)  of  an  associated  dam-fOundation 
system,  the  dam  with  the  foundation  characterized  by  the  stiffness  matrix 
3f(0)  equal  to  the  dynamic  stiffness  matrix  l^(ui)  at  u  ■  0,  were  selected 
as  the  generalized  coordinates  for  response  analysis  of  the  dam  including 
dam-foundation  interaction  effects  (Chapter  5) .  Unlike  the  natural  modes 
of  vibration  of  the  dam  on  a  rigid  base,  these  generalized  coordinates  are 
not  physically  meaningful.  Thus  the  number  of  these  generalized  coordinates 
that  need  to  be  included  in  the  analysis  is  not  as  obvious.  The  complex 
frequency  response  functions  for  the  dam  with  horizontal  ground  accelera¬ 
tion  as  the  excitation  was  obtained  from  four  analyses  using  the  procedure 
of  Chapters  5  and  6,  including  5,  7,  19,  and  15  generalized  coordinates, 
respectively.  The  results  presented  in  Fig.  7.6  are  for  the  idealized  dam 
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monolith,  without  water,  but  including  dam-foundation  interaction  effects. 
The  first  8  eigenvalues  (vibration  frequenices)  of  the  assoicated  dam- 
foundation  system  are  identified  along  the  frequency  axis.  It  is  apparent 
that,  in  order  to  obtain  satisfactory  results  for  responses  over  some  fre¬ 
quency  range,  all  the  generalized  coordinates  associated  with  eigenvalues 
within  the  frequency  range  plus  a  few,  say  three,  more  should  be  included 
in  the  analysis.  Based  on  this  criterion,  the  frequency-range  represented 
in  earthquake  ground  motion,  and  the  fact  that  coup  lex  frequency  response 
functions  decrease  at  larger  excitation  frequencies,  ten  generalized  coordi¬ 
nates  were  considered  sufficient  for  the  results  presented  in  Chapters  8 
and  9. 
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8.1  Scope  of  Chapter 

The  analysis  procedure  developed  in  Chapter  6  can  be  used  to  evaluate 
the  response  of  concrete  gravity  dams  subjected  to  earthquake  ground  motion. 
In  this  procedure,  the  complex  frequency  response  functions  describing  the 
response  to  harmonic  ground  motion  are  determined,  followed  by  the  Fourier 
transform  procedures  to  compute  responses  to  arbitrary  ground  motion.  In 
this  chapter,  the  complex  frequency  response  functions  for  an  idealized 
cross  section  representative  of  concrete  gravity  dams  are  determined  for  a 
range  of  the  important  system  parameters  characterizing  the  properties  of 
the  dam,  foundation  soil  or  rock,  and  impounded  water.  The  effects  of  dam- 
water  interaction  and  of  dam- foundation  interaction  on  the  dynamic  response 
of  the  dam  are  studied. 

8.2  Systems,  Ground  Motions,  Cases  Analyzed  and  Response  Quantities 
8.2.1  Systems 

The  idealized  cross  section  considered  as  representative  of  concrete 

gravity  dams  is  a  triangle  with  a  vertical  upstream  face  and  a  downstream 

face  with  a  slope  of  0.8;1.  The  dam  is  assumed  to  be  homogeneous  with 

linearly  elastic  and  isotropic  properties  for  mass  concrete:  Young's 

modulus  E  *  3,  4  or  5  million  psi;  Poisson’s  ratio  =•  0.2;  and  unit  weight 
s 

=  155  lbs/cu.ft.  Energy  dissipation  in  the  structure  is  represented  by 
constant  hysteretic  damping,  with  the  energy  loss  coefficient  ng  selected  as 
0.1.  This  is  equivalent  to  a  damping  ratio  of  0.05  in  all  the  natural  modes 
of  vibration  of  the  dam  alone  (without  water)  on  a  fixed  base. 

The  finite  element  idealization  for  a  monolith  of  the  dam  is  shown  in 
Fig.  8.1,  consisting  of  20  quadrilateral  elements  and  26  nodal  points,  which 
provide  42  degrees  of  freedom  on  a  fixed  base  and  52  on  a  flexible  base. 

The  dam  is  supported  on  the  surface  of  a  homogeneous,  isotropic, 
linearly  viscoelastic  half  space.  For  the  foundation  material,  several 
values  for  the  Young's  modulus  Ef  are  considered,  Poisson's  ratio  *  1/3,  and 
unit  weight  *  w^  «  165  lbs/cu.ft.  The  foundation  soil  or  rock  is  idealized 
as  a  constant  hysteretic  solid  with  energy  loss  coefficient  nf  ■  0.1. 
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Impounded  water  in  the  reservoir  has  constant  depth  H  and  is  idealized 
as  extending  to  infinity  in  the  upstream  direction.  Several  different  values 
for  H,  the  depth  of  water,  relative  to  Hg/  the  height  of  dam,  are  considered: 
H/Hg  =>  0  (empty  reservoir),  0.8  (partially  filled  reservoir),  and  1  (full 
reservoir).  The  unit  weight  of  water  w  =  62.5  lbs/cu.ft.  and  the  velocity 
of  sound  in  water  =  4720  ft/sec. 

The  dam  and  foundation  are  assumed  to  be  in  a  state  of  generalized  plane 
stress.  This  assumption,  although  not  most  appropriate  for  the  foundation 
material,  is  dictated  by  the  expected  behavior  of  the  dam  (  see  Chapter 
2). 

8.2.2  Ground  Motions 

The  excitation  for  the  dam-water-foundation  system  is  defined  by  the 
two  components  of  free-field  ground  motion  in  the  plane  of  a  monolith  of 
the  dam:  the  horizontal  component  transverse  to  the  dam  axis,  and  the 
vertical  component.  Each  component  of  ground  acceleration  is  assumed  to  be 
harmonic,  a*(t)  =  ^ (t)  =  e*^,  with  the  excitation  frequency  to  be  varied 
over  a  relevant  range. 

8.2.3  Cases  Analyzed 

As  concluded  in  Chapter  7,  the  most  important  parameters  influencing 

the  complex  frequency  response  functions  for  dams  on  rigid  foundations  are 

Eg,  the  Young's  modulus  of  elasticity  for  mass  concrete  in  the  structure, 

and  H/H  ,  the  ratio  of  water  depth  to  dam  height.  When  dam-foundation  inter- 
s 

action  is  considered,  an  additional  parameter  is  important:  E^/Eg,  the  ratio 

of  Young's  moduli  for  the  foundation  and  dam  materials.  For  a  fixed  value 

of  eTe  ,  the  value  of  E  influences  the  response  of  the  structure  to  a 
r  s  s 

minor  extent  for  horizontal  ground  motion  and  significantly  for  vertical 
ground  motion.  However,  to  keep  the  number  of  cases  analyzed  to  a  minimum, 
the  Eg  value  is  not  varied.  Several  dam-water- foundation  systems  were 
defined  by  values  for  parameters  shown  in  Table  8.1.  The  complex  frequency 
response  functions  for  each  case  were  determined.  The  following  figures  are 
presented  in  a  manner  to  facilitate  study  of:  the  effects  of  structure-water 
interaction  and  structure- foundation  interaction,  considered  separately;  and  the 
combined  effects  of  these  two  sources  of  interaction,  on  the  response  of  the 
dam. 
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RESPONSE!  TO  HORIZONTAL  AND  VERTICAL  GROUND  MOTION  SEPARATELY:  CASES  ANALYZED 


if  elasticity  of  concrete  in  the  dam 
>f  elasticity  of  foundation  material 


The  hydrodynamic  pressures  due  to  vertical  ground  motion  depend  on, 

among  other  parameters,  the  reflection  coefficient  at  the  reservoir  bottom 

a  =*  [ (C^w^/Cw)  -  l]/[ (C^w^/Cw)  +  1]  where  and  are  the  unit  weight  and 

P-wave  velocity  for  the  rock,  w  and  C  are  the  unit  weight  and  sound  velocity 

in  water.  It  can  be  easily  seen  that  a  <  1  with  the  maximum  value  correspond- 

£ 

ing  to  rigid  rock.  For  E  *  4  x  10  psi  and  E  /E  *  00 ,  4,  2,  1,  and  the  values 

s  r  s 

of  w  and  w^  mentioned  earlier,  the  values  of  a  are  1,  0.85,  0.80  and  0.72, 
respectively. 

8.2.4  Response  Quantities 

The  response  of  the  dam  to  harmonic  free-field  ground  motion  was  deter¬ 
mined  by  the  analytical  procedures  presented  in  Chapter  6.  Ten  generalized 
coordinates  were  included  in  the  analysis.  Based  on  Sec.  7.6,  the  resulting 
conplex  frequency  response  functions  should  be  accurate  for  excitation  fre¬ 
quencies  up  to  approximately  25  cps.  Although  the  analysis  led  to  response 
at  all  nodal  points,  only  the  results  for  a  few  selected  degrees-of-freedom 
are  presented.  For  dams  on  rigid  foundations,  the  coup  lex  frequency  re¬ 
sponse  functions  for  horizontal  and  vertical  accelerations  at  the  crest  of 
the  dam  are  presented.  In  addition,  the  complex  frequency  response  functions 
for  horizontal  acceleration  at  the  center  of  the  base  of  the  dam  and  for  the 
vertical  accelerations  at  the  upstream  edge  of  the  base  are  presented  when 
interaction  between  the  dam  and  flexible  foundation  material  is  included. 

All  the  conplex  frequency  response  functions  presented  are  for  accel¬ 
erations  relative  to  the  prescribed  free-field  ground  acceleration.  They 
are  not  direct  measures  of  deformations  within  the  dam.  The  absolute  value 
(or  modulus)  of  the  complex  valued  acceleration  response  is  plotted  against 
the  normalized  excitation  frequency  parameter  where  *  fundamental 

fixed-base  natural  frequency  of  the  dam  alone.  When  presented  in  this  form, 
the  results  apply  to  dams  of  all  heights  with  the  idealized  triangular  cross 
section  and  chosen  Es  value.  Furthermore,  if  the  reservoir  is  empty,  the 

plotted  results  are  also  independent  of  the  Eg  value. 

8.3  Complex  Frequency  Responses 

8.3.1  Dam-Water  Interaction  Effects 

In  order  to  identify  the  effects  of  water  on  the  dynamic  response  of 
the  dam,  the  results  of  analysis  cases  1,  2,  and  3  (Table  8.1)  to  horizontal 
and  vertical  ground  motion,  separately,  are  presented  in  Figs.  8.2  and  8.3, 


) 
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8.2  HYDRODYNAMIC  EFFECTS  IN  THE  RESPONSE  OF  DAMS  TO  HARMONIC,  HORIZONTAL  GROUND 
MOTION.  RESULTS  ARE  PRESENTED  FOR  VARYING  WATER  DEPTH:  CASES  1,  2,  AND  3 
OF  TABLE  8.1 
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respectively.  The  response  curve  for  the  dam  without  water  (H/Hg  =  0)  is 
representative  of  a  multi-degree-of-freedom  system  with  constant  mass,  stiff¬ 
ness,  and  damping  parameters.  However,  dam-water  interaction  introduces 
frequency  dependent  terms  in  the  equations  of  motion  of  the  dam,  resulting 
in  complicated  shapes  of  the  response  curves.  The  response  behavior  is 
especially  complicated  at  excitation  frequencies  in  the  neighborhood  of  the 
natural  frequencies  of  water  in  the  reservoir.  In  particular,  the  response 
curve  has  a  double  resonant  peak  near  and  which  is  especially  pro¬ 
nounced  for  the  case  of  vertical  ground  motion.  At  these  frequencies  the 
hydrodynamic  terms  in  the  equations  of  motion  become  unbounded  but  the  re¬ 
sponse  has  finite  limits  (see  Section  6.2.7).  Based  on  the  results  of  Figs. 
8.2  and  8.3,  the  following  observations  can  be  made. 

Presence  of  water  results  in  a  significant  decrease  in  the  fundamental 
resonant  frequency  of  the  dam  but  relatively  little  decrease  in  the  higher 
resonant  frequencies.  At  an  excitation  frequency  w  smaller  than  oj£,  the 
fundamental  resonant  frequency  of  water  in  the  reservoir,  the  effect  of  dam- 
water  interaction  is  equivalent  to  an  added  mass  and  real-valued  load;  their 
magnitude  depends  on  the  excitation  frequency.  This  added  mass,  in  addition 
to  reducing  the  fundamental  resonant  frequency  of  the  dam,  has  the  indirect 
effect  of  reducing  the  apparent  damping  ratio  for  the  fundamental  mode,  re¬ 
sulting  in  narrower  band  width  at  resonance  and  larger  resonant  response. 

Some  of  the  reduction  in  the  resonant  band  width  is,  however,  a  consequence 
of  the  double  resonant  peak.  At  (0  >  wj,  the  effect  of  dam-water  interaction 
is  equivalent  to  frequency  dependent  additional  mass,  damping  and  load;  the 
load  is  complex  valued  for  horizontal  and  vertical  ground  motion  except  that 
it  is  real  valued  for  vertical  ground  motion  with  a  =  1.  The  added  mass  is 
relatively  small,  which  has  relatively  little  effect  on  the  higher  resonant 
frequencies,  but  the  added  damping  is  significant,  resulting  in  decrease  in 
the  response  at  the  higher  resonant  frequencies.  The  fundamental  resonant 
frequency  of  the  dam,  including  hydrodynamic  effects,  is  less  than  both  the 
fundamental  natural  frequency  of  the  dam  with  no  water  and  of  water  in 
the  reservoir,  U)J. 

Comparing  the  results  for  cases  1  (H/Hs  =  0) ,  2  (H/Hs  =  0.8),  and 
3  (H/Hg  *  1),  it  is  seen  that  the  fundamental  resonant  frequency  is  affected 
most  by  water  in  the  upper  parts  of  the  dam  height.  At  higher  excitation 
frequencies,  there  is  little  difference  in  the  responses  for  H/Hs  *  0.8  and  1, 
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except  that  the  response  to  vertical  ground  motion  in  the  neighborhood  of  the 
second  resonant  frequency  is  strongly  affected  by  increasing  H/Hg  from  0.8 
to  1. 

The  vertical  and  horizontal  motions  at  the  crest  of  the  dam  are  obvi¬ 
ously  related  by  the  horizontal  and  vertical  components  of  the  vibration 
mode  shapes.  Vertical  motion  at  the  crest  of  the  dam,  while  smaller  than 
the  horizontal  motion,  as  expected,  is  not  negligible  for  either  horizontal 
or  vertical  ground  motions. 

Comparing  the  responses  of  the  dam  to  horizontal  and  vertical  ground 
motions  (Figs.  8.2  and  8.3),  it  is  apparent,  consistent  with  common  view, 
that  without  water  the  response  to  vertical  ground  motion  is  relatively 
small.  However,  with  a  full  reservoir  the  response  to  vertical  ground  motion 
is  almost  as  large  as  that  due  to  horizontal  ground  motion  when  a  <  1  and 
it  is  even  larger  when  a  -  1.  The  response  to  vertical  ground  motion  is 
large  because  the  hydrodynamic  forces  act  in  the  horizontal  direction  on  the 
vertical  upstream  face  of  the  dam,  causing  significant  lateral  response  of 
the  dam.  Hydrodynamic  effects  are  seen  to  have  especially  large  influence 
on  the  response  to  vertical  ground  motion. 

The  response  of  the  dam  alone  without  water  (H/Hs  »  0) ,  when  presented 
in  the  form  of  Figs.  8.2  and  8.3,  is  independent  of  the  modulus  of  elastic¬ 
ity  Es  of  the  concrete.  Similarly,  the  response  curves  including  hydro- 
dynamic  effects,  do  not  vary  with  Es  if  water  compressibility  is  neglected  [22]. 
However,  results  from  analysis  of  cases  3,  4,  and  5  (Table  8.1)  presented 
in  Figs.  8.4  and  8.5  demonstrate  that  the  Eg  value  affects  the  response 
functions  when  compressibility  of  water  is  included.  This  effect  is 
primarily  on  the  fundamental  resonant  frequency  and  on  the  response  in 
the  neighborhood  of  this  frequency.  The  fundamental  resonant  frequency  of 
the  dam  decreases  due  to  hydrodynamic  effects  to  a  greater  degree  for  the 
larger  values  of  Es.  Increase  in  Es  causes  larger  amplification  of  response 
but  over  a  narrower  frequency  band.  The  resonant  response  to  horizontal 
ground  motion  is  affected  little  by  variations  in  Es  but  is  influenced  sub¬ 
stantially  for  vertical  ground  motion.  At  higher  excitation  frequencies, 
the  response  functions  for  both  components  of  ground  motion  are  insensitive 
to  E3.  With  decreasing  Eg  the  effects  of  water  compressibility  on  response 
become  smaller  and  the  response  approaches  the  incompressible  case. 
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8.3.2  Dam-Foundation  Interaction  Effects 

Complex  frequency  response  functions  for  the  dam,  including  effects  of 
dam- foundation  interaction,  Jure  presented  in  Figs.  8.6  and  8.7  for  varying 
foundation  stiffness  (cases  1,  6,  8  and  11  in  Table  8.1).  When  presented  in 
this  form,  these  functions  do  not  depend  separately  on  Es  and  Ef  but  only  on 
the  ratio  Ef/Eg.  Results  are  presented  for  four  values;  Ef/Es  =  00 ,  4,  2  and 
1.  The  first  represents  rigid  foundation  material,  thus  reducing  the  system 
to  dam  on  fixed  base,  and  the  smallest  value  represents  the  same  elastic 
moduli  for  the  foundation  material  and  dam  concrete. 

Unlike  water  in  the  reservoir,  the  half  plane  does  not  have  any  resonant 
frequencies;  the  foundation  impedances  are  smooth,  slowly  varying  functions 
of  the  frequency  (32].  As  a  result,  structure-foundation  interaction  affects 
the  response  of  the  dam  in  a  simpler  manner  than  does  structure-water  inter¬ 
action.  As  the  Ef/Es  ratio  decreases,  which  for  a  fixed  Es  implies  decrease 
of  foundation  modulus,  the  fundamental  resonant  frequency  of  the  dam 
decreases;  the  response  at  the  crest  of  the  dam  at  this  frequency  decreases 
and  the  frequency  band  width  at  resonance  increases,  implying  an  increase  in 
the  apparent  damping  of  the  structure  (Figs.  8.6  and  8.7).  Accompanying 
this  decrease  in  the  resonant  response  at  the  crest  of  the  dam  is  an  in¬ 
creasing  response  at  the  base  of  the -dam — which  is,  however  a  small  fraction 
of  the  response  at  the  crest--  with  increasingly  flexible  soils  (Figs.  8.6  and 
8.7).  Similar  influences  of  decreasing  the  Ef/Es  ratio  is  observed  at  higher 
resonant  frequencies,  resulting  in  reduced  response  at  the  crest  of  the  dam 
and  increased  response  at  the  base  over  a  wide  range  of  excitation  frequen¬ 
cies;  however,  the  higher  resonant  frequencies  are  decreased  to  a  lesser 
degree  by  structure- foundation  interaction.  The  above  mentioned  effects  of 
structure- foundation  interaction  are  similar  in  response  to  horizontal  and 
to  vertical  ground  motions. 

8.3.3  Dam-Water  and  Dam-Foundation  Interaction  Effects 

In  order  to  understand  the  effects  of  structure-water  interaction  and 
of  structure- foundation  interaction  on  the  complex  frequency  response 
functions  for  the  dam,  the  results  of  several  cases  in  Table  8.1  are  pre¬ 
sented  in  different  ways. 

The  results  of  analysis  cases  8,  9,  and  10  for  horizontal  and  vertical 
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!  OF  THE  RATIO  Ef/Es  ON  RESPONSE  OF  DAMS,  WITHOUT  WATER  BUT  INCLUDING 
lATION  INTERACTION,  TO  HARMONIC,  HORIZONTAL  GROUND  ACCELERATION. 

6,  8,  AND  11  OF  TABLE  8.1 


ground  motion  are  presented  in  Figs.  8.8  and  8.9.  These  results,  for  a 
fixed  Ef/Es  ■  2,  indicate  the  influence  of  varying  depths  of  water  in  the 
reservoir.  Whereas  the  response  function  for  the  dam  without  water 
(H/Hs  *  0)  is  influenced  by  structure-foundation  interaction  effects,  for 
the  other  two  cases  (H/H3  *  0.8,  1)  dam-water  interaction  is  also  involved. 
Qualitatively,  the  effects  of  water  on  the  response  at  the  crest  of  the  dam 
are  generally  similar  whether  the  foundation  material  is  rigid  or  flexible. 
However,  some  of  the  effects  are  relatively  small.  Whereas  the  fundamental 
resonant  frequency  is  decreased  by  roughly  the  same  degree  as  before,  the 
apparent  damping  for  the  fundamental  vibration  mode  is  now  dominated  by 
effects  of  structure- foundation  interaction  and  varies  little  with  the  depth 
of  water.  The  effects  of  dam-water  interaction  on  the  response  of  a  dam  on 
flexible  foundation  are  significant  primarily  at  excitation  frequencies  in 
a  neighborhood  of  the  fundamental  resonant  frequency.  At  higher  frequencies, 
dam-water  interaction  has  little  influence  on  the  response  of  the  dam  on 
flexible  foundation  medium,  except  locally  near  the  resonant  frequencies  of 
the  reservoir.  The  principal  effect  of  dam-water  interaction  on  the 
response  at  the  base  of  the  dam  permitted  by  flexible  foundation  medium  is 
to  decrease  the  fundamental  resonant  frequency  and  to  increase  the  amplitude 
of  response  in  the  neighborhood  of  this  resonant  frequency.  At  higher  fre¬ 
quencies  the  response  at  the  base  of  the  dam  is  not  affected  strongly  by  the 
water,  except  locally  at  the  resonant  frequencies  of  the  water  contained  in 
the  reservoir. 

The  conclusion  presented  in  Sec.  8.3.1  is  also  confirmed  by  the  results 
of  Figs.  8.8  and  8.9.  Without  water  the  response  of  the  dam  to  vertical 
ground  motion  is  a  small  fraction  of  the  response  to  horizontal  ground 
motion,  but  with  water  the  two  responses  at  the  fundamental  resonant  fre¬ 
quency  are  similar;  at  some  higher  frequencies  the  response  to  vertical 
ground  motion  exceeds  that  due  to  horizontal  ground  motion.  Because  of  the 
hydrodynamic  effects,  the  vertical  ground  motion  is  especially  significant 
in  the  response  of  concrete  gravity  dams  to  earthquakes . 

Complex  frequency  response  functions  for  the  dam  with  the  reservoir 
full  of  water  and  varying  modulus  for  foundation  rock  (cases  3,  7,  10,  and 
12  in  Table  8.1)  are  presented  in  Figs.  8.10  and  8.11  for  horizontal  and 
vertical  ground  motion,  respectively.  Results  are  presented  for  a  fixed 
value  of  Es  -  4  x  106  psi  and  four  values  of  Ef/Es  *  00 ,  4,  2,  and  1.  By 
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INFLUENCE  OF  WATER  DEPTH  ON  THE  RESPONSE  OF  DAMS,  INCLUDING  DAM-FOUNDATION 
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INFLUENCE  OF  HATER  DEPTH  ON  THE  RESPONSE  OF  DAMS,  INCLUDING  DAM-FOUNDATION 
INTERACTION,  TO  HARMONIC,  VERTICAL  GROUND  MOTION.  CASES  8,  9,  AND  10  OF 
TABLE  8.1 


DAM  CREST 


INFLUENCE  OF  Ef/Eg  RATIO  ON  RESPONSE  OF  DAMS  WITH  FULL  RESERVOIR  TO  HARMONIC 
HORIZONTAL  GROUND  MOTION.  CASES  3,  7,  10,  AND  12  OF  TABLE  8.1 


HORIZONTAL  ACCLN. 


H/H1«IO 


comparing  the  results  with  those  presented  in  Figs.  8.6  and  8.7,  it  is 
apparent  that  the  effects  of  dam- foundation  interaction  are  generally  sim¬ 
ilar,  independent  of  the  effects  of  dam-water  interaction. 

As  the  Ef/Es  ratio  decreases,  which  for  a  fixed  Es  implies  am  increas¬ 
ingly  flexible  foundation,  the  fundamental  resonant  frequency  decreases. 

The  response  at  the  crest  of  the  dam  at  this  frequency  decreases  and  the 
frequency  bandwidth  at  resonance  increases,  implying  an  increase  in  the 
apparent  damping  of  the  structure.  The  response  behavior  near  the  funda¬ 
mental  resonant  frequency  is  complicated,  especially  for  vertical  ground 
motion,  because  of  the  double  resonant  peak  created  by  hydrodynamic  effects. 
The  response  at  the  base  of  the  dam,  although  only  a  small  fraction  of  the 
response  at  the  crest  of  the  dam,  is  also  affected  by  Ef/Es .  In  the  case  of 
horizontal  ground  motion,  the  base  accelerations  at  and  near  the  fundamental 
resonant  frequency  increase  with  decreasing  Ef/Es,  but  there  is  no  clear 
trend  in  the  case  of  vertical  ground  motion.  Similar  influence  of  increas¬ 
ingly  flexible  soils  is  observed  at  higher  resonant  frequencies,  resulting 
in  reduced  response  at  the  crest  of  the  dam  and  increased  response  at  the 
base  over  a  wide  range  of  excitation  frequencies.  Structure-foundation 
interaction  has  little  influence  on  the  higher  resonant  frequencies  and 
tends  to  suppress  the  response  at  the  dam  crest  at  these  frequencies. 

At  excitation  frequencies  u>  =  ui^,  the  resonant  frequencies  of  the  water 
in  the  reservoir,  the  response  at  the  crest  of  the  dam  is  controlled  by  dam- 
water  interaction  and  is  essentially  independent  of  the  properties  of  the 
foundation  rock  (Figs.  8.10  and  8.11). 

In  order  to  compare  the  effects  of  structure-water  interaction  and 
structure-foundation  interaction  on  the  complex  frequency  response  functions 
for  the  dam,  results  from  analyses  of  four  systems  are  presented  together: 
dam  on  rigid  foundation  with  no  water,  dam  on  rigid  foundation  with  full 
reservoir,  dam  on  flexible  foundation  with  no  water,  dam  on  flexible  founda¬ 
tion  with  full  reservoir.  The  results  for  the  first  two  systems  (cases  1 
and  3  in  Table  8.1)  along  with  those  for  the  latter  two  systems  for  Ef/Es*  4 
(cases  6  and  7) ,  Ef/Es  =  2  (cases  8  and  10) ,  and  Ef/Eg  *  1  (cases  11  and  12) 
are  presented  in  Figs.  8.12  to  8.14  for  horizontal  ground  motion  and  Figs. 
8.15  to  8.17  for  vertical  ground  motion. 
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PIG.  8.12  RESPONSES  OP  DAMS  TO  HARMONIC,  HORIZONTAL  GROUND  MOTION  FOR  FOUR  CONDITIONS: 

DAM  WITH  NO  WATER  ON  RIGID  FOUNDATION  (CASE  1) >  DAM  WITH  FULL  RESERVOIR  ON 
RIGID  FOUNDATION  (CASE  3) ;  DAM  WITH  NO  WATER  ON  FLEXIBLE  FOUNDATION  (Ef/Eg  = 
(CASE  6) }  AND  DAM  WITH  FULL  RESERVOIR  ON  FLEXIBLE  FOUNDATION  (CASE  7) 


FIG.  8.13  RESPONSE  OF  DAMS  TO  HARMONIC,  HORIZONTAL  GROUND  MOTION  FOR  FOUR  CONDITIONS 
DAM  WITH  NO  WATER  ON  RIGID  FOUNDATION  (CASE  1);  DAM  WITH  FULL  RESERVOIR  ON 
RIGID  FOUNDATION  (CASE  3) ;  DAM  WITH  NO  WATER  ON  FLEXIBLE  FOUNDATION  (Ef/Eg 
(CASE  8);  AND  DAM  WITH  FULL  RESERVOIR  ON  FLEXIBLE  FOUNDATION  (CASE  10) 


RESPONSE  OF  DAMS  TO  HARtypNIC,  HORIZONTAL  GROUND  MOTION  FOR  FOUR  CONDITION! 
DAM  WITH  NO  WATER  ON  RIGID  FOUNDATION  (CASE  l)f  DAM  WITH  FULL  RESERVOIR  ON 
RIGID  FOUNDATION  (CASE  3) j  DAM  WITH  NO  WATER  ON  FLEXIBLE  FOUNDATION  (Ef/E( 
(CASE  11) t  AND  DAM  WITH  FULL  RESERVOIR  ON  FLEXIBLE  FOUNDATION  (CASE  12) 


DAM  CREST 


FIG.  8.17  RESPONSE  OF  DAMS  TO  HARMONIC,  VERTICAL  GROUND  MOTION  FOR  FOUR  CONDITIONS; 

DAM  WITH  NO  WATER  ON  RIGID  FOUNDATION  (CASE  1);  DAM  WITH  FULL  RESERVOIR  ON 
RIGID  FOUNDATION  (CASE  3) ;  DAM  WITH  NO  WATER  ON  FLEXIBLE  FOUNDATION  (Ef/Eg 
(CASE  11) ;  AND  DAM  WITH  FULL  RESERVOIR  ON  FLEXIBLE  FOUNDATION  (CASE  12) 


The  fundamental  resonant  frequency  of  the  dam  decreases  because  of  dam- 

water  interaction  and  dam- foundation  interaction,  with  the  influence  of  the 

water  being  larger.  For  the  most  flexible  foundations  considered  (E.  *  E  ) , 

r  s 

the  decrease  in  the  fundamental  frequency  due  to  dam- foundation  interaction 
is  about  the  same  as  that  due  to  dam-water  interaction.  This  trend  may  or 
may  not  exist  for  higher  resonant  frequencies. 

As  mentioned  earlier,  the  response  at  the  crest  of  the  dam  is  increased 

at  the  fundamental  resonant  frequency,  but,  in  many  cases,  decreased  at  higher 

resonant  frequencies  because  of  dam-water  interaction.  These  effects  of 

dam-water  interaction  sure  similar  whether  the  foundation  is  rigid  or  flexible 

for  all  values  of  E./E  considered. 

f  s 

As  mentioned  earlier,  the  response  at  the  crest  of  the  dam  is  reduced  and 
motions  at  the  base  of  the  dam,  relative  to  the  free-field  motion,  are  pos¬ 
sible  due  to  flexibility  of  the  soil.  These  effects  of  dam- foundation  inter¬ 
action  are  similar  whether  the  reservoir  is  empty  or  full.  The  above  mentioned 

effects  of  dam-foundation  interaction  tend  to  increase  as  E./E  decreases,  i.e., 

r  s 

for  a  fixed  Eg,  as  the  foundation  rock  becomes  increasingly  flexible. 

Comparing  the  effects  of  dam-water  interaction  and  dam-foundation  interac¬ 
tion  on  the  resonant  responses  at  the  crest  of  the  dam  to  horizontal  ground 
motion  (Figs.  8.12  to  8.14),  it  is  apparent  that  the  effects  of  the  stiff 
rock  (Ef/Eg  ■  4)  are  similar  to  those  of  the  water,  but  for  softer  rocks 
(E^/Eg  *  2  and  1)  the  effects  of  dam-foundation  interaction  are  relatively 
more  significant.  Whereas  the  amplitude  of  the  response  at  the  crest  of  the 
dam  at  the  fundamental  resonant  frequency  is  increased  by  dam-water  inter¬ 
action  but  reduced  due  to  dam- foundation  interaction,  both  sources  of  inter¬ 
action  contribute  towards  a  decrease  in  the  response  at  higher  resonant  fre¬ 
quencies. 

Because,  as  mentioned  earlier,  dam- water  interaction  influences  the  dam 
response  to  vertical  ground  motion  to  a  much  greater  degree  than  it  influences 
the  response  to  horizontal  ground  motion,  the  relative  effect  of  dam-water 
interaction  and  dam- foundation  interaction  on  dam  response  to  vertical  ground 
motion  (Figs.  8.15  to  8.17)  is  different  than  the  above  observations  from  res¬ 
ponses  to  horizontal  ground  motion  (Figs.  8.12  to  8.14).  At  the  fundamental 
resonant  frequency,  the  response  of  the  dam  is  affected  by  dam-water  interaction 
to  a  greater  degree  than  it  is  affected  by  dam-foundation  interaction.  However, 
at  higher  resonant  frequencies,  the  influence  of  dam-foundation  interaction 
is  greater.  Whereas  the  amplitude  of  the  response  at  the  crest  of  the  dam 
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at  the  fundamental  resonant  frequency  is  increased  by  dam-water  interaction 
but  reduced  due  to  dam-foundation  interaction,  both  sources  of  interaction 
contribute  towards  a  decrease  in  the  response  at  higher  resonant  frequencies, 
but  with  one  exception:  Vertical  acceleration  response  at  the  crest  of  the 
dam  is  increased  by  dam-water  interaction  even  at  the  second  resonant  fre¬ 
quency. 

The  effects  of  water  on  the  response  at  the  base  of  the  dam,  permitted 
by  a  flexible  foundation,  are  closely  tied  to  the  corresponding  effects  on 
the  responses  at  the  crest  of  the  dam.  Dam-water  interaction  has  the  effect 
of  decreasing  the  resonant  frequencies  for  the  base  response  and  modifying  the 
response  amplitudes.  These  effects  are  relatively  large  in  the  responses  to 
vertical  ground  motion. 
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9.  EARTHQUAKE  RESPONSE  OF  PINE  FLAT  DAM 

9.1  Scope  of  the  Chapter 

Chapter  8  presented  the  effects  of  dam  water  Interaction  and  of  dam- 
foundation  interaction,  separately  or  together,  on  the  complex  frequency 
response  functions  of  idealized  gravity  dams  with  triangular  cross-section. 
This  chapter  presents  the  responses  of  Pine  Flat  Dam  to  a  selected  earth¬ 
quake  ground  motion,  analyzed  using  the  procedures  of  Chapter  6  under 
veurious  assumptions.  The  objective  of  this  chapter  is  to  identify  the 
effects  of  dam-water  interaction  and  dam- foundation  interaction  on  the 
response  of  the  dam. 

9. 2  Pine  Flat  Dam,  Ground  Motion,  Cases  Analyzed  and  Response  Results 
9.2.1  Pine  Flat  Dam 

Located  on  the  King's  River  near  Fresno,  California,  Pine  Flat  Dam 
consists  of  thirty-six  50  ft  and  one  40  ft  wide  monoliths.  The  crest 
length  is  1840  ft  and  the  height  of  the  tallest  monolith  is  400  ft.  The 
dam  is  shown  in  Fig.  9.1  and  its  downstream  elevation  is  shown  in  Fig.  9.2. 

The  tallest,  non-overflow  monolith  of  the  dam  with  water  at  El.  951.00 
is  selected  for  purposes  of  this  study.  The  two-dimensional  finite  element 
idealization  for  this  monolith  in  plane  stress  is  shown  in  Fig.  9.3,  con¬ 
sisting  of  136  quadrilateral  elements  with  162  nodal  points.  The  nine  nodal 
points  at  the  base  of  the  dam  are  equally  spaced,  as  required  for  developing 
the  dynamic  stiffness  matrix  for  the  foundation  [321.  This  finite  element 
system  has  306  degrees  of  freedom  analysis  of  the  dam  on  rigid  base.  In 
analyses  considering  the  foundation  flexibility  and  the  resulting  motion  at 
the  dam- foundation  interface,  the  system  has  324  degrees  of  freedom.  The 
mass  concrete  in  the  dam  is  assumed  to  be  a  homogeneous,  isotropic,  elastic 
solid  with  the  following  properties:  Young's  modulus  of  elasticity  -  3.25 
million  psi,  unit  weight  “  155  pcf,  and  Poisson's  ratio  ■  0.2.  The  selected 
elastic  modulus,  determined  by  forced  vibration  tests  on  the  dam  [28],  is 
different  from  that  used  in  earlier  analyses  [22].  The  deunping  factors 
eMsociated  with  the  modes  of  vibration  of  the  dam,  determined  from  the 
above  mentioned  forced  vibration  tests,  were  in  the  range  of  2  to  3.5  per¬ 
cent  of  critical  damping.  A  constant  hysteretic  damping  coefficient  of  0.1, 
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which  corresponds  to  a  5  percent  da&ping  factor  in  all  modes  of  vibration 
of  the  dam,  has  been  selected.  This  is  considered  appropriate  for  the  much 
larger  motions  and  stress  levels  expected  during  strong  earthquake  ground 
shaking. 

The  foundation  rock  is  idealized  as  a  homogenous,  isotropic,  visco¬ 
elastic  half  plane  (in  plane  stress)  of  constant  hysteretic  solid  with  the 
following  properties:  Young's  modulus  of  elasticity  =  10  million  psi,  a  value 
which  is  appropriate  for  the  granites  and  basalts  at  the  site;  Poisson's  ratio 
=  1/3,  unit  weight  *  165  pcf ,  and  damping  coefficient  for  constant  hysteretic 
solid  =»  0.05.  The  P-wave  velocity  for  this  rock  =  10,266  ft/sec. 

The  following  properties  are  assumed  for  water:  unit  weight  *  62.5  pcf, 
wave  velocity  ■  4720  ft/sec.  For  these  properties  of  water  and  rock,  the 
wave  reflection  coefficient  at  the  bottom  of  the  reservoir  is  a  *  0.817. 

9.2.2  Ground  Motion 

The  ground  motion  recorded  at  the  Taft  Lincoln  School  Tunnel  during 
the  Kern  County,  California,  earthquake  of  July  21,  1952,  is  selected  as 
the  excitation  for  analyses  of  Pine  Plat  Dam.  The  ground  motion  acting 
transverse  to  the  axis  of  the  dam  and  in  the  vertical  direction  is  defined 
by  the  S69E  and  vertical  components  of  the  recorded  motion,  respectively. 
These  two  components  of  the  recorded  ground  motion  and  the  maximum  values 
of  acceleration  are  shown  in  Fig.  9.4. 

9.2.3  Cases  Analyzed 

Using  the  confute r  program  described  in  Appendix  A  based  on  the  anal¬ 
ysis  procedure  presented  in  Chapter  6,  responses  of  the  tallest  monolith  of 
Pine  Flat  Dam  to  the  selected  ground  motion  are  analyzed.  The  dam  monolith 
and  the  foundation  are  assumed  to  be  in  plane  stress. 

Responses  of  the  dam  to  the  S69E  component  of  the  Taft  ground  motion 
only  and  to  the  S69E  and  vertical  components  acting  simultaneously  are 
analyzed.  For  each  of  these  excitations  the  response  of  the  dam  is  sep¬ 
arately  analyzed  four  times,  corresponding  to  the  following  four  sets  of 
assumptions  for  the  foundation  and  hydrodynamic  effects. 
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GROUND  MOTION  AT  TAFT  LINCOLN  SCHOOL  ' 
CALIFORNIA,  EARTHQUAKE,  JULY  21,  1952 


CURVE  WATER  FOUNDATION 
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FIG.  9.5  RESPONSES  OF  PINE  FLAT  DAM  TO  HARMONIC,  HORIZONTAL  GROUND  MOTION 

COMPUTED  FOR  FOUR  CONDITIONS:  DAM  WITH  NO  WATER  ON  RIGID  FOUNDATION, 
DAM  WITH  FULL  RESERVOIR  ON  RIGID  FOUNDATION,  DAM  WITH  NO  WATER  ON 
FLEXIBLE  FOUNDATION,  AND  DAM  WITH  FULL  RESERVOIR  ON  FLEXIBLE  FOUNDATION 
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Rigid 

Included 
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Flexible 

Excluded 

4 

Flexible 

Included 

In  order  to  consider  hydrodynamic  effects  realistically,  compressibil¬ 
ity  of  water  is  included.  The  displacements  and  stresses  due  to  the  weight 
of  the  dam  and  to  hydrostatic  pressures  acre  included  in  all  analyses.  As 
mentioned  in  Chapter  6,  the  rigid  body  displacements  of  the  dam  due  to  a 
deformable  foundation  have  been  excluded. 

All  the  vibration  modes  or  generalized  coordinates,  as  appropriate, 
necessary  to  obtain  accurate  results  for  complex  frequency  response  up  to 
excitation  frequencies  of  approximately  20  cps  were  included  in  the  analysis. 
The  first  five  natural  modes  of  vibration  of  the  dam  were  included  in  the 
analyses  assuming  a  rigid  base.  The  first  ten  generalized  coordinates, 
were  included  in  analyses  considering  dam- foundation  interaction  effects 
(see  Sec.  6.2.6) . 

9.2.4  Response  Results 

Hie  complex  frequency  response  functions  for  the  horizontal  accelera¬ 
tion  at  the  crest  of  Pine  Flat  Dam,  with  horizontal  ground  motion  as  the 
excitation,  for  the  four  sets  of  assumptions  for  foundation  and  hydrodynamic 
effects  (Sec.  9.2.3)  are  presented  in  Fig.  9.5.  From  these  results,  the 
fundamental  resonant  period  of  vibration  and  the  effective  damping,  deter¬ 
mined  by  the  half-power  band  width  method,  are  as  follows: 


Hydrodynamic 
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Fundamental 

Resonance 

Case 

Foui  ;ion 
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0.317 
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Excluded 

0.341 

6.7 
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Flexible 

Included 

0.429 

5.2 
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These  vibration  periods  are  identified  on  the  response  spectrum  of  the 
S69E  component  of  the  Taft  ground  motion  (Fig.  9.6),  The  spectrum  ordinates 
at  these  vibration  periods  corresponding  to  damping  ratio  of  5  percent  can 
be  observed  in  the  same  figure. 

The  response  of  Pine  Flat  Dam  was  determined  for  each  of  the  two  exci¬ 
tations  listed  in  Sec.  9.2.2  and  the  four  sets  of  assumptions  in  Sec.  9.2.3. 
In  each  cue,  the  conputer  results  of  dynamic  analysis  represent  the  total 
response,  including  the  effects  of  the  weight  of  the  dam  and  hydrostatic 
pressures.  These  results  consisted  of  the  complete  time-history  of  the 
horizontal  and  vertical  components  of  displacement  of  all  the  nodal  points 
and  of  three  cotqponents  of  stress  in  all  the  finite  elements.  For  each  of 
the  cases  analyzed  only  a  small  part  of  the  total  result  is  presented.  The 
maximum  crest  displacement  and  maximum  tensile  stresses  in  three  critical 
parts  of  the  monolith  are  summarized  in  Table  9.1.  Presented  in  Figs.  9.7 
to  9.17  are  the  time-history  of  displacements  at  nodal  points  1,  73,  and 
118,  located  at  different  levels  on  the  upstream  face,  and  at  nodal  points 
154,  158,  and  162  at  the  base  when  foundation  flexibility  is  considered; 
and  the  distribution  of  envelope  values  of  maximum  principal  stress  (maxi¬ 
mum  tensile  stress  or  minimum  compressive  stress)  during  the  earthquake . 

9.3  Dam-Water  Interaction  Effects 

Without  hydrodynamic  effects,  the  response  of  the  dam  is  typical  of  a 
multi-degree-of-freedom  system  with  mass,  stiffness,  and  damping  properties 
independent  of  excitation  frequency.  Dam-water  interaction  introduces 
frequency  dependent  terms  in  the  equations  of  motion,  resulting  in  compli¬ 
cated  response  curves,  especially  in  the  neighborhood  of  the  resonant 
frequencies  for  the  impounded  water  (Fig.  9.5).  A  decrease  in  the  funda¬ 
mental  resonant  frequency  of  the  dam  is  apparent  from  a  comparison  of  the 
responses  with  and  without  hydrodynamic  effects. 

When  the  excitation  is  the  S69E  component  of  Taft  ground  motion,  the 
inclusion  of  hydrodynamic  effects  increases  the  maximum  displacement  at  the 
crest  of  the  dam  from  1.38" to  1.83" (Fig.  9.7).  The  maximum  tensile  stresses 
in  the  dam  are  increased  at  the  upstream  face  from  153  to  223  psi,  from  208 
to  254  psi  at  the  downstream  face,  and  from  257  to  366  psi  at  the  heel 
(Fig.  9.8).  The  area  enclosed  by  a  particular  stress  contour  increases 
due  to  hydrodynamic  effects,  indicating  that  stresses  exceed  the  value 
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FIG.  9.6  RESPONSE  SPECTRUM  FOR  THE  S69E  COMPONENT  OF  TAFT  GROUND 
MOTION.  THE  FUNDAMENTAL  VIBRATION  PERIOD  OF  PINE  FLAT 
DAM  COMPUTED  FOR  THE  FOUR  SETS  OF  ASSUMPTIONS  (SEC.  9.2.4) 
FOR  FOUNDATION  AND  HYDRODYNAMIC  EFFECTS  ARE  AS  NOTED 
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DISPLACEMENT  RESPONSE  OF  PINE  FLAT  DAM  ON  RIGID  FOUNDATION  TO  S69E 
PONENT  OF  TAFT  GROUND  MOTION:  (A)  HYDRODYNAMIC  EFFECTS  EXCLUDED,  i 
(B)  HYDRODYNAMIC  EFFECTS  INCLUDED 


9.8  ENVELOPE  VALUES  OF  MAXIMUM  PRINCIPAL  (MAXIMUM  TENSILE  OR  MINIMUM 
COMPRESSIVE)  STRESSES  IN  PINE  FLAT  DAM  ON  RIGID  FOUNDATION  DUE  T< 
S69E  COMPONENT  OF  TAFT  GROUND  -MOTION:  (A)  HYDRODYNAMIC  EFFECTS 
EXCLUDED,  AND  (B)  HYDRODYNAMIC  EFFECTS  INCLUDED. 


corresponding  to  that  stress  contour  over  a  larger  portion  of  the  monolith. 

Consistent  with  the  observations  in  Sec.  8.3.1,  when  hydrodynamic 
effects  sure  excluded,  the  response  of  the  dam  is  only  slightly  increased  by 
the  contribution  of  the  vertical  component  of  ground  motion  (compare  Figs. 
9. 7- 9. 8  with  9.9-9.10).  But,  contrary  to  the  observations  in  Sec.  8.3.1 
srnd  earlier  results  for  Pine  Flat  Dam  (18),  the  displacements  and  stresses 
in  the  dam  with  hydrodynamic  effects  included  are  slightly  decreased  by 
the  contributions  of  the  vertical  component  of  ground  motion  (compare  Figs. 
9. 7- 9. 8  with  9.9-9.10). 

Earlier  results  (18)  for  the  response  of  Pine  Flat  Dam  subjected  to 
the  same  Taft  ground  motion  but  assuming  a  different  value  of  Young's 
modulus,  demonstrated  that  the  vertical  component  of  ground  motion  causes 
considerable  increase  in  the  response  of  the  dam.  Ifce  vertical  ground 
motion  causes  lateral  hydrodynamic  forces  resulting  in  significant  lateral 
displacements  and  associated  stresses.  Because  the  principal  change  in  the 
system  properties  and  ground  motion  from  the  earlier  results  to  those  pre¬ 
sented  here  is  in  the  value  of  assumed  Young's  modulus  for  the  dam  (the 
other  change  is  the  damping  model  discussed  in  Sec.  7.5.1);  it  is  surpris¬ 
ing  that  the  results  presented  here  are  not  consistent  with  the  earlier 
conclusions.  In  order  to  resolve  what  appears  to  be  an  anomaly,  the  earth¬ 
quake  stresses  (excluding  initial  static  stresses)  due  to  horizontal  and 
vertical  ground  motions,  separately,  and  including  hydrodynamic  effects,  for 
selected  finite  elements  (Nos.  40,  41,  129  in  Fig.  9.3),  are  presented  in 
Fig.  9.11.  Although  considerable  stresses  are  caused  by  vertical  ground 
motion,  they  partially  cancel  the  stresses  due  to  horizontal  ground  motion, 
resulting  in  reduced  response  when  both  ground  motion  components  are  con¬ 
sidered  simultaneously.  The  contribution  of  the  vertical  conponent  of 
ground  motion  to  the  total  response  of  a  dam,  including  hydrodynamic 
effects,  therefore  depends  on  the  relative  phasing  of  the  responses  to 
horizontal  and  vertical  ground  motion,  which  in  turn,  depends  on  the 
phasing  of  the  ground  motion  components  and  the  vibration  properties  of  the 
dam. 

9.4  Dam-Foundation  Interaction  Effects 

The  foundation  impedances  for  a  half  plane  are  smooth,  slowly  varying 
functions  of  the  excitation  frequency.  However,  the  hydrodynamic  terms  are 
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FIG.  9.9  DISPLACEMENT  RESPONSE  OF  PINE  FLAT  DAM  ON  RIGID  FOUNDATION  TO  S69E  AND 
VERTICAL  COMPONENTS,  SIMULTANEOUSLY,  OF  TAFT  GROUND  MOTION:  (A)  HYDRO- 
DYNAMIC  EFFECTS  EXCLUDED ,  AND  (B)  HYDRODYNAMIC  EFFECTS  INCLUDED 


MAXIMUM  PRINCIPAL  STRESS 


GROUND  MOTION 
S69E  COMPONENT 


ELEMENT-41 


ELEMENT-129 


0  I  2  3  4  5  6  7  8  9  10  II  12  13 

TIME-SECONDS 


FIG.  9.11  STRESSES  AT  SELECTED  LOCATIONS  OF  PINE  FLAT  DAM  ON  RIGID 

FOUNDATIONS  DUE  TO  THE  SEPARATE  ACTION  OF  S69E  AND  VERTICAL 
COMPONENT  OF  TAFT  GROUND  MOTION.  HYDRODYNAMIC  EFFECTS 
INCLUDED.  STATIC  STRESSES  EXCLUDED. 
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unbounded  at  the  natural  frequencies  of  water  in  the  reservoir.  As  a 
result,  structure- foundation  interaction  influences  the  response  of  the 
in  a  siitpler  manner  than  the  hydrodynamic  effects  do  (Fig.  9.5).  The  funda¬ 
mental  resonant  frequency  is  decreased  and  the  corresponding  damping  ratio 
is  increased  due  to  structure- foundation  interaction  (Fig.  9.5).  This  de¬ 
crease  in  frequency  is,  however,  smaller  than  the  frequency  reduction  due 
to  hydrodynamic  effects.  The  ordinate  of  the  pseudo-acceleration  response 
spectrum  for  the  S69E  component  of  Taft  ground  motion  is  essentially  unaf¬ 
fected  by  the  increase  in  vibration  period  but  would  be  reduced  due  to  the 
increased  damping  (Fig.  9.6).  This  leads  to  increase  in  displacements 
(compare  Figs.  5.7a  and  5.12)  and  reduction  in  stresses  (compare  Figs.  9.8a  4 

and  9.13a).  The  displacements  increase  due  to  lengthening  of  the  vibration 
period  (Figs.  9.7a  and  9.12).  Flexibility  of  the  soil  permits  motions  at 
the  base  of  the  dam  (Fig.  9.12)  but  these  are  much  smaller  than  the  motions  * 

at  the  crest  of  the  dam.  Comparison  of  Figs.  9.8a  and  9.13a  indicates  that 
the  stresses  near  the  base  of  the  dam  are  relaxed  because  of  soil  flexibil¬ 
ity. 

The  effects  of  structure- foundation  interaction  on  the  response  of  the  ^ 

dam  to  horizontal  and  vertical  components  of  ground  motion  acting  simulta¬ 
neously  can  be  observed  by  coiq?aring  Figs.  9.9a  and  9.14  with  Figs.  9.10a 
and  9.15a,  respectively.  These  effects  are  generally  similar  to  those  t 

observed  above  in  the  responses  to  horizontal  ground  motion  alone.  The 
contributions  of  the  vertical  component  of  ground  motion  to  the  total 
response  of  the  dam  sure  rather  small  (compare  Fig.  9.13a  and  9.15a). 

9.5  Dam-Water  and  Dam- Foundation  Interaction  Effects 

Structure- foundation  interaction  affects  the  response  of  the  dam  in  a 
similar  manner,  reducing  the  fundamental  resonant  frequency  and  increasing 
the  effective  damping,  whether  hydrodynamic  effects  are  included  or  not 
(Fig.  9.5).  The  fundamental  resonant  frequency  is  reduced  by  dam- foundation 
interaction  and  by  dam-water  interaction.  The  two  reductions  are  additive, 
resulting  in  further  reduction  in  frequency  when  both  interaction  effects 
are  considered  (Fig.  9.5),  and  in  considerably  smaller  resonant  response  as 
compared  with  either  the  response  of  the  dam  alone  or  the  response  includ¬ 
ing  hydrodynamic  effects. 
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FIG.  9.12  DISPLACEMENT  RESPONSE  OF  PINE  FLAT  DAM,  INCLUDING  DAM- FOUND ATI i 
INTERACTION,  BUT  EXCLUDING  HYDRODYNAMIC  EFFECTS,  TO  S69E  COMPO 
OF  TAFT  GROUND  MOTION 


ENVELOPE  VALUES  OF  MAXIMUM  PRINCIPAL  (MAXIMUM  TENSILE  OR  MINIMUM 
COMPRESSIVE)  STRESSES  IN  PINE  FLAT  DAM  ON  FLEXIBLE  FOUNDATION 
DUE  TO  S69E  COMPONENT  OF  TAFT  GROUND  MOTION:  (A)  HYDRODYNAMIC 
EFFECTS  EXCLUDED,  AND  (B)  HYDRODYNAMIC  EFFECTS  INCLUDED. 
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FIG.  9.14  DISPLACEMENT  RESPONSE  OF  PINE  FLAT  DAM,  INCLUDING  DAM-FOUNDATION 

INTERACTION  BUT  EXCLUDING  HYDRODYNAMIC  EFFECTS,  TO  S69E  AND  VERTICAL 
COMPONENTS,  SIMULTANEOUSLY ,  OF  TAFT  GROUND  MOTION 


However,  the  displacements  and  stresses  due  to  either  excitation — 
horizontal  ground  motion  only  or  horizontal  and  vertical  ground  motions 
simultaneously-- are  considerably  increased  due  to  structure -foundation 
interaction,  (compare  Figs.  9.7b  and  9.16,  9.8b  and  9.13b,  9.9b  and  9.17, 

9.10b  and  9.19b).  Compared  to  the  response  of  the  dam  including  only 
hydrodynamic  effects,  the  stresses  in  upper  parts  of  the  dam  are  increased 
due  to  structure- foundation  interaction  when  the  excitation  is  only  the 
horizontal  component  of  ground  motion  (Figs.  9.8b  and  9.13b);  also,  when 
the  excitation  included  the  vertical  component  of  ground  motion  (Figs.  9.10b 
and  9.15b).  Stresses  at  the  heel  of  the  dam  are  increased  to  a  much  lesser 
extent  because  of  the  stress-relaxation  due  to  foundation  flexibility.  The 
area  enclosed  by  a  particular  stress  contour  increases  due  to  dam- foundation 
interaction,  indicating  that  the  dam  is  stressed  beyond  that  contour  value  over 
a  larger  portion.  This  increase  in  earthquake  response  occurs  in  spite  of 
the  decreased  response  indicated  by  complex  frequency  responses,  because 
structure- foundation  interaction  shifts  the  fundamental  resonant  period  to 
correspond  with  a  peak  of  the  response  spectrum  (Fig.  9.6).  Thus,  the 
effects  of  dam- foundation  interaction  and  dam-water  interaction  on  earth¬ 
quake  response  of  the  dam  depend  partly  on  the  relative  ordinates  of  the 
response  spectrum  at  resonant  periods  of  the  dam  with  and  without  these 
interaction  effects. 

The  effects  of  structure-foundation  interaction  on  the  response  of  the 
dam  are  generally  similar  with  or  without  the  vertical  component  of  ground 
motion.  For  reasons  mentioned  in  Sec.  9.3,  the  contributions  of  the  vert¬ 
ical  conponent  of  ground  motion  to  the  response  of  the  dam,  including  both 
sources  of  interaction,  are  rather  small. 
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FIG.  9.16  DISPLACEMENT  RESPONSE  OF  PINE  FLAT  DAM,  INCLUDING  DAM-FOUNDATION 
INTERACTION  AND  HYDRODYNAMIC  EFFECTS,  TO  S69E  COMPONENT  OF  TAFT 
GROUND  MOTION 
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PIG.  9.17  DISPLACEMENT  RESPONSE  OF  PINE  FLAT  DAM,  INCLUDING  DAM-FOUNDATION 
INTERACTION  AND  HYDRODYNAMIC  EFFECTS,  TO  S69E  AND  VERTICAL  COMPOI 
SIMULTANEOUSLY,  OF  TAFT  GROUND  MOTION 
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10.  CONCLUSIONS 


The  method  presented  for  two-dimensional  linear  analysis  of  earthquake 
response  of  gravity  dam  monoliths  is  effective  for  practical  problems. 

Included  in  this  analysis  are  the  effects  of  dynamic  interaction  among  dam, 
impounded  water,  and  foundation  rock  on  the  dam  response.  The  dam,  impounded 
water,  and  foundation  rock  are  considered  as  three  substructures  of  the  com¬ 
plete  system.  Displacements  of  the  dam  jure  expressed  as  linear  combinations 
of  generalized  coordinates,  which  are  selected  as  the  normal  modes  of  an 
associated  dam-foundation  system.  Results  can  be  obtained  to  any  desired 
degree  of  accuracy  by  including  the  necessary  number  of  generalized  coordi¬ 
nates.  The  substructure  approach,  combined  with  transformation  of  displace¬ 
ments  to  generalized  coordinates,  leads  to  an  efficient  analysis  procedure. 

The  analysis  procedure  is  developed  specifically  for  two-dimensional 
analyses  of  gravity  dam  monoliths  supported  on  the  surface  of  a  viscoelastic 
halfplane  and  inpounding  water  in  the  reservoir  with  horizontal  bottom.  The 
dam  is  discretized  as  a  two-dimensional  finite  element  system,  but  the  founda¬ 
tion  and  fluid  domains  are  treated  as  continua.  The  general  approach  and 
concepts  embodied  in  the  substructure  method  are  also  applicable  to  more  com¬ 
plex  systems  where,  because  of  irregular  geometry  and/or  nonhomogeneous  mate¬ 
rial  properties,  the  foundation  and/or  fluid  domains  must  also  be  discretized. 

In  general,  dam-water  interaction,  including  water  compressibility  and 
dam-foundation  interaction,  have  significant  influence  on  the  response  of  the 
dam  and  should  be  considered  in  the  analysis.  Because  the  foundation  impedances 
are  smooth,  slowly  varying  functions  of  excitation  frequency,  but  the  hydro- 
dynamic  terms  are  unbounded  at  the  resonant  frequencies  of  the  fluid  domain, 
structure-foundation  interaction  affects  the  response  of  the  dam  in  a  simpler 
manner  than  does  structure-water  interaction. 

Dam-foundation  interaction  effects  in  the  response  of  dams  depend  on  several 

factors,  the  most  iaportant  of  which  is  the  ratio  Ef/Eg  of  the  elastic  modulii 

of  the  foundation  and  dam  materials.  As  the  E./E  ratio  decreases  (which  for 

r  s 

a  fixed  Eg  implies  a  decrease  of  the  foundation  modulus) ,  the  fundamental  res¬ 
onant  frequency  of  the  dam  decreases;  the  response  at  the  crest  of  the  dam  at 
this  frequency  decreases;  and  the  frequency  bandwidth  at  resonance  increases, 
inlying  an  increase  in  the  apparent  damping  of  the  structure.  The  influence 
of  decreasing  the  E./E  ratio  is  similar  at  higher  resonant  frequencies;  however, 
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the  higher  resonant  frequencies  are  decreased  to  a  lesser  degree  by  structure- 
foundation  interaction.  The  effects  of  structure- foundation  interaction  are 
essentially  independent  of  the  effects  of  dam-water  interaction,  except  at  < 

the  resonant  frequencies  of  the  fluid  domain.  The  dam  response  at  these 
frequencies  is  controlled  by  dam-water  interaction  and  is  essentially  inde¬ 
pendent  of  the  properties  of  the  foundation  rock. 

The  influence  of  structure- foundation  interaction  on  the  earthquake 
response  of  a  dam  depends  in  part  on  the  change  in  the  earthquake  response 
spectrum  ordinate  due  to  the  decrease  in  frequency  and  increase  in  damping. 

In  addition,  the  stresses  near  the  base  of  the  dam  are  relaxed  because  of 
foundation  flexibility. 

The  frequency  response  curves  are  complicated  in  the  neighborhood  of  the 
natural  frequencies  of  water  in  the  reservoir.  In  particular,  the  response 
curves  have  a  double  resonant  peak  near  the  fundamental  frequencies  of  the 
dam  and  fluid  domain,  considered  separately.  In  comparison  with  the  response 
behavior  of  a  dam  without  water,  the  response  in  the  fundamental  mode  exhibits 
highly  resonant  behavior  when  the  'reservoir  is  full?  the  peak  response  at  the 
fundamental  frequency  is  especially  large  compared  to  that  at  higher  resonant 
frequencies.  Because  of  dam-water  interaction,  the  fundamental  resonant  fre¬ 
quency  of  the  dam  is  decreased  by  an  amount  depending  on  the  depth  of  water  — 
with  water  in  the  upper  parts  of  the  dam  having  the  most  influence  —  and 
modulus  of  elasticity  of  the  dam.  The  higher  resonant  frequencies  of  the  dam 
are  reduced  relatively  little  by  dam-water  interaction.  Qualitatively,  the 
effects  of  water  on  the  dam  response  are  generally  similar  whether  the  founda¬ 
tion  rock  is  rigid  or  flexible.  Hydrodynamic  effects  reduce  the  fundamental 
resonant  frequency  of  the  dam,  including  dam- foundation  interaction,  by  rough¬ 
ly  the  same  degree,  independent  of  the  foundation  material  properties.  How¬ 
ever,  the  apparent  damping  for  the  fundamental  vibration  mode  is  dominated  by 
effects  of  structure-foundation  interaction  and  varies  little  with  the  depth 
of  water. 

The  displacements  and  stresses  of  Pine  Plat  Dam  due  to  the  Taft  ground 
motion  are  increased  significantly  because  of  hydrodynamic  effects. 

The  fundamental  resonant  frequency  of  the  dam  is  reduced  by  dam-water 
interaction  and  by  dam- foundation  interaction.  The  two  reductions  are  addi¬ 
tive  with  the  influence  of  the  water  usually  being  larger.  This  trend  may 
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or  may  not  exist  at  higher  resonant  frequencies.  Similarly  there  are  no 
general  trends  regarding  the  comparative  effects  of  the  water  and  foundation 
on  the  resonant  responses  of  the  dam.  Which  effect  is  more  significant  de¬ 
pends  on  the  E-/E  ratio,  the  depth  of  water,  the  order  of  the  resonant  fre- 
r  s 

quency  (fundamental  or  higher),  and  the  ground  motion  component  (horizontal  or 
vertical) . 

Compared  to  the  response  of  Pine  Flat  Dam  including  only  hydrodynamic 
effects,  the  stresses  in  the  upper  parts  of  the  dam  are  significantly  increased 
due  to  structure- foundation  interaction.  Stresses  at  the  heel  of  the  dam 
are  increased  to  a  much  lesser  extent  because  of  the  stress  relaxation  due  to 
foundation  flexibility. 

The  response  of  the  dam,  without  water,  to  vertical  ground  motion  is  small 
relative  to  that  due  to  horizontal  ground  motion,  but  it  becomes  relatively 
signficant  when  hydrodynamic  effects  are  included.  However,  the  contribution 
of  the  vertical  component  of  ground  motion  to  the  total  response,  including 
hydrodynamic  effects,  depends  on  the  relative  phasing  of  the  responses  to 
horizontal  and  vertical  ground  motion,  which  in  turn  depends  on  the  phasing 
of  the  ground  motion  components  and  the  vibration  properties  of  the  dam. 
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APPENDIX  A  -  USERS  GUIDE  TO  COMPUTER  PROGRAM 


IDENTIFICATION 

EAGD  Earthquake  Analysis  of  Gravity  Dams 

Programmed:  P.  Chakrabarti ,  Sunil  Gupta ,  G.  Dasgupta 


PURPOSE 

This  computer  program  has  been  developed  to  determine  the  elastic  dynamic 
response  of  monoliths  of  concrete  gravity  dams  to  earthquake  ground  motion. 
The  effects  of  dam-water  interaction  and  dam-foundation  interaction  on  the 
response  of  a  dam  are  included.  The  computer  program  is  based  on  the  sub¬ 
structure  analysis  procedure  developed  in  Chapter  6,  in  which  the  dam  is 
treated  as  a  finite  element  system,  the  impounded  water  as  continum,  and  the 
foundation  rock  or  soil  region  as  a  finite  element  system  or  a  viscoelastic 
halfspace,  as  appropriate  for  the  site  conditions. 

The  computer  program  is  applicable  to  two-dimensional  structural  systems. 
Thus,  the  dam  monolith  and  foundation  must  be  idealized  as  a  system  in  plane 
stress  or  plane  strain.  Compressiblity  of  water  is  recognized  in  the  anal¬ 
ysis.  The  excitation  includes  the  transverse  (to  dam  axis)  and  vertical 
components  of  free-field  ground  motion,  assumed  to  be  the  same  across  the 
base  of  the  dam. 

In  order  to  define  the  input  to  the  computer  program,  a  cross-section  or 
monolith  of  the  dam  must  be  idealized  as  an  assemblage  of  planar,  quadri¬ 
lateral  finite  elements  as  shown  by  the  example  later  in  this  appendix. 
Elements  in  the  idealization  are  identified  by  a  sequence  of  numbers  start¬ 
ing  with  one.  All  nodal  points  are  identified  by  a  separate  numbering 
sequence  starting  with  one. 

Because  of  dam- foundation  interaction,  a  dynamic  stiffness  matrix  for  the 
foundation  appears  in  the  equation  of  motion.  This  matrix  depends  on  the 
excitation  frequency  and  is  defined  with  respect  to  the  degrees-of- freedom 
of  nodal  points  at  the  dam- foundation  interface.  This  matrix  must  be  deter¬ 
mined  by  a  separate  analysis  and  provided  to  the  computer  program  as  am 
input . 

The  computer  program  is  written  in  FORTRAN  IV  and  was  develped  on  the  CDC 
6400  Computer  at  the  University  of  California,  Berkeley.  It  can  be  executed 
on  other  computers  with  minimal  changes.  However,  the  Fast  Fourier  Transform 
package  included  is  partly  written  in  COMPASS,  developed  specifically  for  the 
CDC  6400  computer.  In  using  this  computer  prograun  on  other  coiqputer  systems, 
a  comparable  subroutine  package  for  FFT  computations  should  be  provided. 
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INPUT  DATA 


The  following  sequence  of  punched  cards  and  data  on  a  tape  numerically 
define  the  dam-water- foundation  system,  the  ground  motion  and  control  para¬ 
meters  for  the  anlaysis. 

A.  TITLE  CARD  (8A  W) 


Columns  1-80:  Contain  title  or  any  information  to  be  printed  with 
results. 


B. 


CONTROL 

CARD 

(515 

,  F10.0 

Columns 

1  - 

5 

NUMNP : 

6  - 

10 

NUMEL: 

11  - 

15 

NUMMAT: 

16  - 

20 

NBASE: 

21  - 

25 

NEV: 

26  - 

35 

WL: 

36  - 

40 

NPP : 

41  - 

45 

IGRAV: 

46  - 

50 

PSP: 

51  - 

55 

IRES: 

56  - 

60 

IOPR: 

61  - 

65 

I0PP: 

66  - 

70 

IRIG: 

215,  F5.0,  515) 

Number  of  nodal  points. 

Number  of  elements. 

Number  of  different  materials. 

Number  of  nodal  points  at  base  of  the  dam. 

Number  of  eigenvalues  (NEV  £  MBAND  +  1) . 

If  NEV  *  0  and  IGRAV  ^  0  only  static  analysis 
is  performed. 

Water  level  in  the  reservoir,  in  feet. 

Number  of  nodal  points  on  the  upstream  face 
of  the  dam  affected  by  the  water  pressure 
(see  G) . 

^  0,  to  perform  static  analysis. 

=  0.,  if  plane  stress  problem. 

=  1.,  if  plane  strain  problem. 

=  0,  to  perform  dynamic  response  analysis, 
otherwise,  only  static  analysis  and  mode 
shapes  are  computed. 

0,  to  skip  calculation  of  frequencies  and 
mode  shapes.  They  are  read  from  cards  in 
this  case. 

j  0,  to  punch  frequencies  and  mode  shapes  on 
cards . 

?  0,  if  the  foundation  is  rigid. 


C.  FOUNDATION  CARD  (3F10.0) 

Omit  this  card  if  IRIG  j*  0,  (i.e.  for  rigid  foundation). 
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Columns  1-10 


E:  Modulus  of  elasticity  of  foundation,  in  Ksf. 


11  -  20  RHO:  Mass  density  of  foundation,  in  K-sec^/ft  4. 

21  -  30  RBASE:  Spacing  between  equally-spaced  nodal  points 

at  base  of  the  dam. 

D.  DAM  MATERIAL  PROPERTY  CARD  (I5,3F10.0) 

The  following  cards  must  be  supplied  for  each  different  material 
(NUMMAT  cards). 

Columns  1  -  5  Material  identification  number. 

6-15  Modulus  of  elasticity,  in  Ksf. 

16  -  25  Posson's  ratio. 

2 

26  -  35  Mass  density  of  material,  in  K-sec  /ft  4. 

E.  NODAL  POINT  CARDS  (I5,F5.0,2F10.0,2I5,2F10.0) 

Columns  1  -  5  Nodal  Point  Number. 

6-10  Boundary  condition  code  "p". 

11  -  20  X-ordinate,  in  ft. 

21  -  30  Y-ordinate,  in  ft. 

31  -  60  Used  for  layer  generation,  otherwise  leave  blank. 

Specification  for  Code  "p" : 

p  «  0.  Both  displacements  unknown. 

p  ■  1.  Zero  displacement  in  the  X-direction. 

Unknown  displacement  in  the  Y-direction. 

p  *  2.  Unknown  displacement  in  the  X-direction. 

Zero  displacement  in  the  Y-direction. 

p  ■  3.  Zero  displacement  in  the  X-direction. 

Zero  displacement  in  the  Y-direction. 

Nodal  point  cards  must  be  in  numerical  sequence.  If  cards  are  omitted 
and  Cols.  31  -  60  are  left  blank,  the  omitted  nodal  points  are  generated 
along  a  straight  line  between  the  defined  nodal  points  (see  Note  1) . 

Or,  if  Cols.  31  -  60  are  used  they  are  generated  in  layers  (see  Note  2) . 
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Note  1:  Straight  line  generation. 

If  the  (L-I)  cards  for  nodal  points  N+l,  N+2  ....N+L-l  are  omitted  and 
Cols.  31-  60  of  the  card  for  nodal  point  N  are  left  blank,  the  omitted 
nodal  points  are  generated  at  equal  intervals  on  the  straight  line 
joining  nodes  N  and  (N+L) . 

Note  2:  Layer  generation. 

Layer  generation  may  be  used  after  two  rows  of  nodal  points  are  com  - 
pletely  defined.  If  on  the  card  for  node  N  the  following  data  is 
specified: 


Columns  31  -  35 

MOD: 

Module,  m  (>  0). 

36  -  40 

NLIM: 

Limit  of  generation  (>  N) . 

41  -  50 

FACX: 

Amplification  factor  fx. 

(if  left  blank,  assumed  to  be  1) 

51  -  60 

FACY: 

Amplification  factor  fy. 

(if  left  blank,  assumed  to  be  1) 


the  X-Y  coordinates  of  points  N+l ,  N+2  ....  NLIM  are  generated  by  the 
formulas 


X,_  +  f 

k-m  X 


«*-»  '  W 


Y,  =  Y,  +  f„  (Y,  -  Y,  „  ) 

k  k-m  Y  k-m  k-2m 


for  k  =  N+l,  ....  NLIM.  If  NLIM  =  NUMNP  no  more  nodal  cards  are  needed. 
If  NLIM  <  NUMNP,  the  card  for  point  (NLIM+1)  must  follow. 


The  boundary  condition  code  for  generated  nodal  points  is  set  equal  to 
zero. 


F.  ELEMENT  CARDS  (615) 


Columns  1  -  5 


6  - 

10 

11  - 

15 

16  - 

20 

21  - 

25 

26  - 

30 

Element  number 

The  maximum  difference  "b"  between 
these  numbers  is  an  indication  of 
the  band  width  of  the  Stiffness 
Matrix.  ”b"  may  be  minimized  by  a 
judicious  numbering  of  nodal  points. 

Material  identification 


Nodal  Point  I 
Nodal  Point  J 
Nodal  Point  K 
Nodal  Point  L 


For  a  right-hand  coordinate  system  the  nodal  point  number  I,  J,  K,  and  L 
must  be  in  sequence  in  a  counter-clockwise  direction  around  the  element 
(see  Fig.  A.l).  Element  cards  must  be  in  element  number  sequence.  If 
element  cards  are  omitted  the  program  automatically  generated  the 
omitted  information  by  incrementing  by  one  the  preceedir.g  I,  J,  K,  and 
L.  The  material  identification  for  the  generated  card  is  set  equal 
to  the  corresponding  value  on  the  last  card.  The  last  element  card 


must  always  be  supplied.  Triangular  elements  are  also  permissible; 
they  are  identified  by  repeating  the  last  nodal  number  (i.e.  I,  J,  K, 
K). 


G.  WATER  PRESSURE  CARDS  (1615) 

These  cards  are  necessary  only  if  NPP  >  0. 

Columns  1  -  5  =  0,  if  the  reservoir  is  to  the  left  of  the  dam. 

=  1,  if  the  reservoir  is  to  the  right  of  the  dam. 


6-10 
11  -  15 
16  -  20 
etc. 


H.  BASE  CONNECTION  NODE  CARDS  (1615) 

Omit  if  IRIG  ?  0  (i.e.  for  a  rigid  foundation). 

Nodal  point  numbers  at  the  base  of  the  dam  which 
are  connected  to  the  foundation  rock,  starting 
from  left  to  right;  a  total  of  NBASE  numbers,  as 
shown  in  Fig.  A. 3. 


I.  FREQUENCY  CARDS  (I12,F15.8) 

This  array  is  necessary  only  if  the  control  variable  (Card  B)  IOPR  ^  0. 
Columns  1-12  mode  number. 

13  -  22  frequency  in  rad/sec. 


J.  MODE  SHAPE  CARDS  (I12,2E15.6) 

This  array  is  also  necessary  only  if  the  control  variable  (Card  B) 

IOPR  jt  0.  The  mode  shapes  of  the  associated  dam- foundation  system  have 
to  be  normalized  in  the  sense  m  ^  *  I,  where  is  the  mode  shape 
matrix,  m  is  the  mass  matrix,  and"l"the~identity~matrix.  If  the  array 
is  the  output  of  a  previous  run  of”this  program  on  puched  cards,  then 
the  mode  shapes  are  already  normalized.  This  array  is  as  follows: 


Columns  1  -  5 

6-10 
11  -  15 


Nodal  point  numbers  affected  by  the 
water  starting  from  the  top  as  shown 
in  Fig.  A. 2. 


FIG.  A.  3  NODAL  POINTS  AT  DAM-  FIG.  A.  4  DEFINITION  OF  N 

FOUNDATION  INTERFACE  AND  "GRACE  BAND' 


One  card  for  each  nodal  point 


Columns  1-12 
13  -  27 
28  -  42 


Nodal  point  number 

X-ordinate 

X-ordinate 


one  set  for  each 
mode  shape. 


t 


K.  RESPONSE  CONTROL  CARD  (3I5,2F10.0) 


Columns  1  -  5  IHVs  Code  for  ground  motion  component 

*  0,  if  only  horizontal  ground  accelerations 
are  to  be  included. 

*  1,  if  only  vertical  ground  accelerations 
are  to  be  included. 

=  2,  if  both  the  horizontal  and  vertical 
ground  accelerations  are  to  be  included. 


6-10  NEXP:  Defines  the  number  of  values  of  excitation 
frequency  for  which  complex  frequency 
response  is  to  be  evaluated.  The  number  of 
excitation  frequency  values  is  2NEXP. 


11  -  15  IHYD:  Code  for  hydrodynamic  interaction.  IHYD  ?  0 
if  hydrodynamic  interaction  is  to  be 
neglected.  In  this  case  only  hydrostatic 
effects  are  included.  (If  NPP  «  0  and 
WL  *  0  in  card  B  then  this  must  be  non -zero) . 

16  -  25  DT:  Time  interval  at  which  dynamic  response  is 
generated.  The  frequency  increment  for 
computing  complex  frequency  response  is 
determined  from  DT  and  NEXP.  (See  note  3.) 


26  -  35  ALPHA:  Coefficient  of  reflection  between  water  and 
ground  rock  below  reservoir.  This  may  be 
computed  as  ALPHA  *  (k-l)/(k+l),  where 
k  *  Cr  wr/Cw  with  wr  and  w  neing  the  unit 
weights  of  rock  and  water  respectively,  Cr 
the  P-wave  velocity  in  rock  and  C  the 
velocity  of  sound  in  water  (4720  ft/sec) . 
ALPHA  <  1. 


Note  3:  Choice  of  DT  and  NEXP. 

For  using  the  FFT  algorithm  any  earthquake  ground  motion  input  data 
(see  Card  Group  M)  is  interpolated  at  equal  intervals  of  time  At  =  DT. 
The  total  number  of  ordinates  is  determined  by  NEXP,  and  is  given  by 
^  »  2 NEXP.  Thus,  the  total  duration  of  ground  excitation  used  in  the 
computations  is  j  ■  At  •  ^  (See  Fig.  A. 4.)  A  "grace  band"  of  ^ 
additional  zeros  are  appended  to  the  record  to  minimize  the  aliasing 
errors  in  the  FFT  computations.  The  uniform  frequency  increment  Af 
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and  the  frequency  range  2  ^or  which  the  complex  frequency  response 
function  has  to  be  computed  is  given  by  the  following  relationships 


N  *r 
2  *  Af 


After  the  Fourier  transform  computations,  which  include  the  forward 
transform  of  the  ground  motion  and  the  inverse  transforms  of  the  re¬ 
sponses  are  obtained  at  the  first  -=r  points  spaced  at  equal  intervals 
of  time  At.  By  specifying  TD  less  than  T  (see  Card  Group  M) ,  response 
may  be  obtained  for  a  shorter  duration.  As  is  shown  in  Fig.  A. 4  the 
second  part  of  the  response  contains  free  vibration  response  of  the 
structure  after  the  ground  excitation  has  ceased. 

It  has  been  stated  earlier  that  At(i.e.,  DT)  and  N  (or  NEXP)  determine 
the  values  of  T,  Af  and  F.  It  is  crucial  to  choose  these  appropriately 
so  that  maximum  accuracy  in  the  results  is  ensured.  The  following  con¬ 
siderations  govern  the  selection  of  the  time  step  At  and  total  number 
of  points 

1.  The  frequency  increment  is  small  enough  to  permit  an  accurate 
description,  especially  near  the  resonant  frequencies, .of  the 
frequency  responses  in  the  generalized  coordinates.  It  is 
recommended  that 


Af  <  st 

where  f^  is  the  fundamental  frequency  of  vibration  of  the  dam  in 
cps.  An  estimate  of  this  may  be  obtained  from  a  preliminary  mode 
shape  analysis  of  the  structure  using  this  computer  program. 

The  frequency  range  should  be  large  enough  to  include  the  signifi¬ 
cant  portion  of  the  frequency  response  function  for  generalized 
displacement  in  the  highest  mode  of  vibration  included  in  the 
analysis.  If  f£  is  the  frequency  of  vibration  of  the  highest  mode 
included  in  the  analysis,  it  is  recommended  that 

I  >  2  fl 

The  frequency  range  should  include  all  the  significant  frequencies 
contained  in  typical  earthquake  records.  Earthquake  data  processed 
by  modern  techniques  accurately  reproduces  frequencies  up  to  about 
25  cps.  Based  on  these  considerations,  in  general,  it  is  recom¬ 
mended  that 


—  >  25  cps. 
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M.  GROUND  ACCELERATION  INFORMATION  CARDS. 

1.  CONTROL  CARD  (2I5,F10.0) 

Columns  1  -  5  NXUGH:  Number  of  ordinates  describing  time- 

history  of  horizontal  ground  accelera¬ 
tion. 

6-10  NXUGVs  Number  of  ordinates  describing  time- 

history  of  vertical  ground  acceleration. 

11-20  TD:  Time  duration  for  which  response  is 

desired.  TD  <  DT.  (2NEXP-  1) . 

2.  HORIZONTAL  GROUND  ACCELERATION  CARDS  (6 (F6. 3 ,F6.4) ) 

These  cards  are  to  be  omitted  if  IHV  -  1.  NXUGH  time-acceleration 
pairs  describing  the  time-history  of  transverse  horizontal  ground 
acceleration  are  to  be  specified  on  these  cards,  with  six  pairs 
per  card.  Time  must  be  expressed  in  seconds  and  the  acceleration 
as  multiples  of  g,  the  acceleration  due  to  gravity. 
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3.  VERTICAL  GROUND  ACCELERATION  CARDS  (6(F6.3,F6.4) ) 


These  cards  are  to  be  omitted  if  IHV  *  0.  NXUGV  time  acceleration 
points  describing  the  time-history  of  vertical  ground  acceleration 
jure  to  be  specified  on  these  cards,  with  six  pairs  per  card.  Time 
must  be  expressed  in  seconds  and  the  acceleration  as  multiples  of 
g,  the  acceleration  due  to  gravity. 

N.  OUTPUT  INFORMATION  CARDS 

1.  OUTPUT  CONTROL  CARD  (515) 

Columns  1-5  NPRINTt  Print  interval. 

Nodal  point  displacements  and  element 
stresses  are  printed  every  NPRINT  time 
intervals. 

6-10  ICOMBs  =0,  if  dynamic  response  is  needed 
separately. 

=  1,  if  the  static  and  dynamic  responses 
are  to  be  combined. 

11  -  15  ISEL:  Selection  code 

r1  0,  if  output  is  desired  for  only 
selected  nodal  points  and  elements; 
otherwise  displacements  of  all  nodal 
points  and  stresses  in  all  elements 
are  printed. 

16  -  20  NNODE:  Total  number  of  nodes  for  which  dis¬ 
placements  are  to  be  printed,  >_  1. 

21  -  25  NNELM:  Total  number  of  elements  for  which 

stresses  are  to  be  printed,  >_  1. 

NNODE  and  NNEIM  may  be  left  blank  if  ISEL  *  0. 

2.  NODAL  POINT  SELECTION  CARDS  (1615) 

These  cards  are  to  be  omitted  if  ISEL  *  0  or  NNODE  *  0.  List  the 
NNODE  nodal  point  numbers  at  which  displacements  are  to  be  printed. 

3.  ELEMENT  SELECTION  CARDS  (1615) 

These  cards  are  to  be  omitted  if  ISEL  «  0  or  NNELM  *  0.  List  the 
NNELM  element  numbers  for  which  stresses  are  to  be  printed. 

O.  FOUNDATION  MATRIX  TAPE 

If  IRIG  ■  0,  the  dynamic  stiffness  matrix  for  the  foundation  region 
should  be  available  on  a  tape.  To  be  compatible  with  this  program, 
this  data  should  have  been  generated  by  Computer  Programs  described 
in  "Dynamic  Stiffness  Matrices  for  Homogenous  Viscoelastic  Half  Planes," 
by  G.  Dasgupta  and  A.  K.  Chopra,  Report  No.  UCB/EERC-77/26;  otherwise 
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the  input  statements  in  this  program  should  be  modified  by  the  user. 


The  data  read  from  the  tape  should  then  be  copied  onto  a  scratch  file 
called  TAPE90,  which  is  read  by  the  computer  program.  After  execution, 
the  data  is  destroyed. 


OUTPUT 


* 


The  following  is  printed  by  the  program.  (Note  that  some  of  these  may  be 
suppressed  according  to  the  options  provided  in  card  B.) 

1.  First  set  of  input  data:  structural  and  material  properties, 
foundation  properties,  options,  etc. 

2.  Hydrostatic  Loads:  i.e.  equivalent  nodal  point  loads  due  to  hydro¬ 
static  pressure  of  water  in  the  reservoir. 

Nodal  point  displacements  and  element  stresses  for  static  loads. 

3.  Frequencies  and  Mode  Shapes. 

4.  Second  set  of  input:  response  data  including  damping,  time  step, 
etc. 

5.  Absolute  value  of  complex  frequency  responses  for  generalized 
accelerations  of  the  dam  due  to  horizontal  and  vertical  ground 
motions  separately. 

6.  Third  set  of  input:  earthquake  acceleration  data,  etc. 

7.  Displacements  of  selected  nodal  points  (according  to  ISEL  and 
NNODE  in  Card  Group  N)  and  stresses  in  selected  elements  (accord¬ 
ing  to  ISEL  and  NNELM  in  Card  Group  N)  at  instants  of  time  deter¬ 
mined  by  the  print  interval  (NPRINT  in  Card  Group  N) . 

8.  The  peak  values  of  major  and  minor  principal  stresses  in  each 
element  during  the  earthquake  and  the  times  at  which  they  occur. 

9.  The  following  quantities  are  written  on  tape  unit  3. 

Logical  Record  1:  NUMNP ,  NUMEL,  NEV,  ND,  DT.ND  is  the  total  number 
of  time  intervals  of  response  determined  from  TD  and  DT  using 
integer  arithmetic. 

Starting  with  Record  2,  for  each  of  the  ND  time  intervals,  two 
records  are  written  using  the  following  two  statements: 

WRITE  (3)  X 

WRITE  (3)  STRES 

where  X  and  STRES  are  one-dimensional  arrays,  dimensioned  properly, 
so  that  X (2  *  I  -  1)  and  X(2  *  I)  are  the  X-component  and  the 
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Y-component  of  the  displacement  at  node  number  I,  I  -  1,2  ....  NUMNP ,  Sim¬ 
ilarly  STRES  (3  *  N  -  2) ,  STRES  (3  *  N  -  IK  and  STRES  (3  *  N)  are  the  three 
components  of  stress,  i.e.  crx,  ay,  and  CT^y  in  the  element  number  N,  N*l,2 
....  NUMEL.  A  physical  tape  has  to  be  requested,  with  the  file  name  TAPE3,  in 
the  program  control  cards  if  these  results  are  to  be  saved  for  subsequent  use. 
Plotting  programs  are  available  to  read  from  this  tape  and  plot  time-history 
of  nodal  point  displacements  and  element  stresses. 


EXAMPLE 

The  preparation  of  input  data  for  this  program  is  illustrated  by  an  example. 
A  very  coarse  finite  element  idealization  of  a  cross-section  of  a  concrete 
gravity  dam  is  shown  in  Fig.  A. 5,  followed  by  the  input  data  required  for 
the  analysis  of  this  structure  subjected  to  horizontal  earthquake  motion. 


COMPUTER  PROGRAM  INFORMATION,  LIMITATIONS,  AND  TIMING 

The  program  is  written  in  FORTRAN  IV  and  was  developed  on  a  CDC6400  com¬ 
puter  at  the  University  of  California,  Berkeley.  The  core  storage  require¬ 
ments  of  the  program  are  separated  into  fixed  and  variable  parts  with  the 
fixed  part  consisting  of  instructions,  non-subscripted  variables,  and  those 
arrays  which  do  not  depend  on  the  size  of  the  individual  problem.  The  vari¬ 
able  part  is  stored  in  array  A  which  appears  in  the  blank  COMMON  statement. 

The  high  speed  storage  requirements  of  the  program  can  be  changed  depending 
on  the  size  of  the  problem  to  be  solved.  This  is  done  by  changing  two 
FORTRAN  statements  at  the  start  of  the  program,  i.e. 

COMMON  A  (N) 

MS TOR  »  N 

The  value  specified  for  N  must  exceed  each  of  the  following. 

1.  N0  +  2  *NUMNP*  MBAND+6  *NB  2 

2.  N0  +  41*NUMEL 

3.  Nq  +  6  *NUMNP+NEV*  NEV 

4.  N0  +  3*NNB+2*NB2+NPP+4*NTERM+NEV*NEV* (4+NTERM) 

+NEV* (2*NUMNP+3*NCOMP+2*NWR*NCOMP+2*NTERM+2) 

5.  NDATA* (4+2*NEV)+NEV* (1+2*NWR)+NXUG 

6.  NUMNP*(5+NEV)+NEV* (2*NWR+1)+42*NUMEL 
where, 


N0  ■  6*NUMNP+NBC+NPP+NEV+NNB+NB2 


NODAL  POINT  NUMBER 


FIG.  A. 5  EXAMPLE  FINITE  ELEMENT  MESH  FOR  A  CONCRETE  GRAVITY  DAM  TO 
ILLUSTRATE  PREPARATION  OF  INPUT  DATA  FOR  COMPUTER  PROGRAM 
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and 


NUMNP  =  Number  of  nodal  points 

NUMEL  =  Number  of  elements 

NBC  *  Number  of  displacement  constratints 

NPP  =  Number  of  nodal  points  affected  by  water  pressure 

NCOMP  =  Number  of  components  of  earthquake  included  in  the  analysis. 

NDATA  -  2  **  NEXP 

NEV  *  Number  of  generalized  coordinates  included. 

NEXP  *  Exponent  for  the  number  of  ordinates  to  define  complex 
requency  responses. 

MBAND  *  2  *  (MB  +  1) 

MB  **  max  MB^ ,  i  *  1,  NUMEL 

MBi  *  difference  between  the  largest  and  smallest  nodal  point 
numbers  for  element  i. 

NTERK  *  50 

NXUG  -  max  (fclXUGH,  NXUGV) 

NXUGH  *  Number  of  horizontal  ground  acceleration  ordinates 
NXUGV  *=  Number  of  vertical  ground  acceleration  ordinates 
NBASE  *  Number  of  nodal  points  at  base  of  the  dam. 

NB2  “  NBASE *2 

NNB  =  NB2*NB2*2 

NWR  -  max(100,Nt>ATA+l) 

If  only  frequencies  and  mode  shapes  are  desired  along  with  static  analysis 
it  suffices  to  check  (1)  and  (2)  above. 

The  computer  time  required  for  solution  depends  on  a  number  of  factors. 

The  more  important  ones  are  the  number  of  time  ordinates  (which  is  deter¬ 
mined  by  NEXP)  the  number  of  nodal  points,  the  bandwidth  of  the  stiffness 
matrix  (the  nodal  points  should  be  numbered  in  a  manner  which  minimizes 
the  bandwidth),  the  number  of  elements,  the  number  of  generalized  coordin¬ 
ates  to  be  included,  and  the  interval  for  printing  and  writing  the  output. 

Some  representative  execution  times  on  a  CDC  6400  computer  are  given  for 
four  cases  of  analysis  (Tables  A.l  and  A. 2)  for  Pine  Flat  dam  (Chapter  8), 
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with  NUMNP  =  162,  NUMEL  =  136,  MBAND  =  22,  and  NEXP  =10,  NEV  =  5 
for  cases  1  and  2,  10  for  cases  3  and  4,  NBASE  =  9.  Table  A.l  pertains  to 
analysis  for  transverse  ground  motion  only.  Table  A. 2  is  for  analyses  con¬ 
sidering  transverse  and  vertical  ground  motions,  simultaneously.  The  results 

include  the  complete  history  of  displacements  and  stresses  during  the  earth¬ 
quake  . 

Tables  A.l  and  A. 2  provide  some  indication  of  the  increase  in  computational 
effort  required  to  include  dam-water  and  dam.- foundation  interaction  effects 
in  the  analysis.  Considering  that  these  interaction  effects  complicate  the 
analysis  greatly,  the  additional  computer  time  required  is  modest.  Compar¬ 
atively,  more  computational  effort  is  required  to  include  dam-foundation 
interaction  effects,  primarily  because  a  larger  number  of  generalized 
coordinates  need  to  be  included  in  the  analysis.  The  efficiency  of  the 
analysis  lies  in  use  of  the  substructure  method  along  with  transformation 
to  generalized  coordinates. 


TABLE  A. 1  COMPUTATION  TIMES:  PINE  FLAT  DAM , HORI ZONTAL  GROUND  MOTION  ONLY 


CASE 

FOUNDATION 

HYDRODYNAMIC 

EFFECTS 

NO.  OF 

GENERALIZED 

COORDINATES 

CENTRAL 
PROCESSOR 
TIME  (SEC.) 

REMARKS 

1 

RIGID 

EXCLUDED 

5 

314 

Eigenvalues , 
Eigenvectors, 
and 

Dynamic 
Foundation 
Stiffness 
Matrix 
are  read 

as 

input 

2 

RIGID 

INCLUDED 

5 

474 

3 

FLEXIBLE 

EXCLUDED 

10 

810 

4 

FLEXIBLE 

INCLUDED 

10 

m3  J 

TABLE  A. 2  COMPUTATION  TIMES:  PINE  FLAT  DAM,  HORIZONTAL  AND  VERTICAL 

GROUND  MOTION 


CASE 

FOUNDATION 

HYDRODYNAMIC 

EFFECTS 

NO.  OF 

GENERALIZED 

COORDINATES 

CENTRAL 
PROCESSOR 
TIME  (SEC.) 

REMARKS 

■ 

RIGID 

EXCLUDED 

5 

325 

Eigenvalues, 

Eigenvectors 

and 

Dynamic 
Foundation 
Stiffness 
Matrix 
are  read 

as 

input 

2 

RIGID 

INCLUDED 

5 

500 

3 

FLEXIBLE 

EXCLUDED 

10 

832 

fl 

FLEXIBLE 

INCLUDED 

10 

1356 

APPENDIX  B 

LISTING  OF  EAGD  COMPUTER  PROGRAM 
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